Development of a Conceptual Site Model for Evaluating Seawater Intrusion under Sea Level Rise Scenarios Using Analytical Methods

Island County is located in the Puget Sound of Washington State and includes several islands, the largest of which is Whidbey Island. Central Whidbey Island was chosen as the project site, as residents use groundwater for their water supply and seawater intrusion near the coast is known to contaminate this resource. In 1989, Island County adopted a Saltwater Intrusion Policy and used chloride concentrations in existing wells in order to define and map “risk zones.” In 2005, this method of defining vulnerability was updated with the use of water level elevations in conjunction with chloride concentrations. The result of this work was a revised map of seawater intrusion vulnerability that is currently in use by Island County. This groundwater management strategy is defined as trigger-level management and is largely a reactive tool. In order to evaluate trends in the hydrogeologic processes at the site, including seawater intrusion under sea level rise scenarios, this report presents a workflow where groundwater flow and discharge to the sea are quantified using a revised conceptual site model. The revised conceptual site model used several simplifying assumptions that allow for first-order quantitative predictions of seawater intrusion using analytical methods. Data from water well reports included lithologic and well construction information, static water levels, and aquifer tests for specific capacity. Results from specific capacity tests define the relationship between discharge and drawdown and were input for a modified Theis equation to solve for transmissivity (Arihood, 2009). Components of the conceptual site model were created in ArcGIS and included interpolation of water level elevation, creation of groundwater basins, and the calculation of net recharge and groundwater discharge for each basin. The revised conceptual site model was then used to hypothesize regarding hydrogeologic processes based on observed trends in groundwater flow. Hypotheses used to explain a reduction in aquifer thickness and hydraulic gradient were: (1) A large increase in transmissivity occurring near the coast. (2) The reduced aquifer thickness and hydraulic gradient were the result of seawater intrusion. (3) Data used to create the conceptual site model were insufficient to resolve trends in groundwater flow. For Hypothesis 2, analytical solutions for groundwater flow under Dupuit assumptions were applied in order to evaluate seawater intrusion under projected sea level rise scenarios. Results indicated that a rise in sea level has little impact on the position of a saltwater wedge; however, a reduction in recharge has significant consequences. Future work should evaluate groundwater flow using an expanded monitoring well network and aquifer recharge should be promoted by reducing surface water runoff.

Comparisons of physical separation methods of Kunjin virus-induced membranes

The two sets of connected membranes induced in Kunjin virus-infected cells are characterized by the presence of NS3 helicase/protease in both, and by RNA-dependent RNA polymerase (RdRp) activity plus the associated double-stranded RNA (dsRNA) template in vesicle packets (VP), or by the absence of both the VP-specific markers in the convoluted membranes/paracrystalline arrays (CM/PC). Attempts were made to separate flavivirus-induced membranes by sedimentation or flotation analyses in density gradients of sucrose or iodixanol, respectively, after treatment of cell lysates by sonication, osmotic shock, or tryptic digestion. Only osmotic shock treatment provided suggestive evidence of separation. This was explored by flow cytometry analysis (FCA) of RdRp active membrane fractions from a sucrose gradient, using dual fluorescent labelling via antibodies to NS3 and dsRNA. FCA revealed the presence of a dual labelled membrane population indicative of VP, and in a faster sedimenting fraction a membrane population able to be labelled only in NS3, representative of CM/PC and associated (R)ER. It was postulated that osmotic shock ruptured the bounding membrane of the VP, releasing the enclosed small vesicles associated with the Kunjin virus replication complex characterized previously. Notably, the presence of the full spectrum of nonstructural proteins in some membrane fractions was not a reliable marker for RdRp activity. These experiments may provide the opportunity for isolation of relatively pure flavivirus replication complexes in their native membrane-associated state by fluorescence-activated cell sorting. (C) 2004 Elsevier B.V. All rights reserved.

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

The design of planar gradient coils. Part ll: A weighted superposition method

In this work a superposition technique for designing gradient coils for the purpose of magnetic resonance imaging is outlined, which uses an optimized weight function superimposed upon an initial winding similar to that obtained from the target field method to generate the final wire winding. This work builds on the preliminary work performed in Part I on designing planar insertable gradient coils for high resolution imaging. The proposed superposition method for designing gradient coils results in coil patterns with relatively low inductances and the gradient coils can be used as inserts into existing magnetic resonance imaging hardware. The new scheme has the capacity to obtain images faster with more detail due to the deliver of greater magnetic held gradients. The proposed method for designing gradient coils is compared with a variant of the state-of-the-art target field method for planar gradient coils designs, and it is shown that the weighted superposition approach outperforms the well-known the classical method.

A gradient descent algorithm for minimizing amino acid coupling reactions when synthesizing cyclic-peptide libraries

Combinatorial chemistry has become an invaluable tool in medicinal chemistry for the identification of new drug leads. For example, libraries of predetermined sequences and head-to-tail cyclized peptides are routinely synthesized in our laboratory using the IRORI approach. Such libraries are used as molecular toolkits that enable the development of pharmacophores that define activity and specificity at receptor targets. These libraries can be quite large and difficult to handle, due to physical and chemical constraints imposed by their size. Therefore, smaller sub-libraries are often targeted for synthesis. The number of coupling reactions required can be greatly reduced if the peptides having common amino acids are grouped into the same sub-library (batching). This paper describes a schedule optimizer to minimize the number of coupling reactions by rotating and aligning sequences while simultaneously batching. The gradient descent method thereby reduces the number of coupling reactions required for synthesizing cyclic peptide libraries. We show that the algorithm results in a 75% reduction in the number of coupling reactions for a typical cyclic peptide library.

A wave equation technique for designing compact gradient coils

A primary purpose of this research is to design a gradient coil that is planar in construction and can be inserted within existing infrastructure. The proposed wave equation method for the design of gradient coils is novel within the field. it is comprehensively shown how this method can be used to design the planar x-, y-, and z-gradient wire windings to produce the required magnetic fields within a certain domain. The solution for the cylindrical gradient coil set is also elucidated. The wave equation technique is compared with the well-known target held method to gauge the quality of resultant design. In the case of the planar gradient coil design, it is shown that using the new method, a set of compact gradient coils with large field of view can be produced. The final design is considerably smaller in dimension when compared with the design obtained using the target field method, and therefore the manufacturing costs and materials required are somewhat reduced.

Convergence analysis of a discrete-time single-unit gradient ICA algorithm

We revisit the one-unit gradient ICA algorithm derived from the kurtosis function. By carefully studying properties of the stationary points of the discrete-time one-unit gradient ICA algorithm, with suitable condition on the learning rate, convergence can be proved. The condition on the learning rate helps alleviate the guesswork that accompanies the problem of choosing suitable learning rate in practical computation. These results may be useful to extract independent source signals on-line.

Particle tracking in a salinity gradient: A method for measuring sinking rate of individual phytoplankton in the laboratory

This paper presents a new method to measure the sinking rates of individual phytoplankton “particles” (cells, chains, colonies, and aggregates) in the laboratory. Conventional particle tracking and high resolution video imaging were used to measure particle sinking rates and particle size. The stabilizing force of a very mild linear salinity gradient (1 ppt over 15 cm) prevented the formation of convection currents in the laboratory settling chamber. Whereas bulk settling methods such as SETCOL provide a single value of sinking rate for a population, this method allows the measurement of sinking rate and particle size for a large number of individual particles or phytoplankton within a population. The method has applications where sinking rates vary within a population, or where sinking rate-size relationships are important. Preliminary data from experiments with both laboratory and field samples of marine phytoplankton are presented here to illustrate the use of the technique, its applications, and limitations. Whereas this paper deals only with sinking phytoplankton, the method is equally valid for positively buoyant species, as well as nonbiological particles.

Reconstruction of a stationary flow from boundary data using iterative methods

In this paper, three iterative procedures (Landweber-Fridman, conjugate gradient and minimal error methods) for obtaining a stable solution to the Cauchy problem in slow viscous flows are presented and compared. A section is devoted to the numerical investigations of these algorithms. There, we use the boundary element method together with efficient stopping criteria for ceasing the iteration process in order to obtain stable solutions.

Microarray-formatted clinical biomarker assay development using peptide aptamers to anterior gradient-2

Anterior gradient-2 protein was identified using proteomic technologies as a p53 inhibitor which is overexpressed in human cancers, and this protein presents a novel pro-oncogenic target with which to develop diagnostic assays for biomarker detection in clinical tissue. Combinatorial phage-peptide libraries were used to select 12 amino acid polypeptide aptamers toward anterior gradient-2 to determine whether methods can be developed to affinity purify the protein from clinical biopsies. Selecting phage aptamers through four rounds of screening on recombinant human anterior gradient-2 protein identified two classes of peptide ligand that bind to distinct epitopes on anterior gradient-2 protein in an immunoblot. Synthetic biotinylated peptide aptamers bound in an ELISA format to anterior gradient-2, and substitution mutagenesis further minimized one polypeptide aptamer to a hexapeptide core. Aptamers containing this latter consensus sequence could be used to affinity purify to homogeneity human anterior gradient-2 protein from a single clinical biopsy. The spotting of a panel of peptide aptamers onto a protein microarray matrix could be used to quantify anterior gradient-2 protein from crude clinical biopsy lysates, providing a format for quantitative screening. These data highlight the utility of peptide combinatorial libraries to acquire rapidly a high-affinity ligand that can selectively bind a target protein from a clinical biopsy and provide a technological approach for clinical biomarker assay development in an aptamer microarray format.

Survival and growth responses of eight Everglades tree species along an experimental hydrological gradient on two tree island types

Questions: How are the early survival and growth of seedlings of Everglades tree species planted in an experimental setting on artificial tree islands affected by hydrology and substrate type? What are the implications of these responses for broader tree island restoration efforts? Location: Loxahatchee Impoundment Landscape Assessment (LILA), Boynton Beach, Florida, USA. Methods: An experiment was designed to test hydrological and substrate effects on seedling growth and survivorship. Two islands – a peat and a limestone-core island representing two major types found in the Everglades – were constructed in four macrocosms. A mixture of eight tree species was planted on each island in March of 2006 and 2007. Survival and height growth of seedlings planted in 2006 were assessed periodically during the next two and a half years. Results: Survival and growth improved with increasing elevation on both tree island substrate types. Seedlings' survival and growth responses along a moisture gradient matched species distributions along natural hydrological gradients in the Everglades. The effect of substrate on seedling performance showed higher survival of most species on the limestone tree islands, and faster growth on their peat-based counterparts. Conclusions: The present results could have profound implications for restoration of forests on existing landforms and artificial creation of tree islands. Knowledge of species tolerance to flooding and responses to different edaphic conditions present in wetlands is important in selecting suitable species to plant on restored tree islands

Changing Bacterial Growth Efficiencies across a Natural Nutrient Gradient in an Oligotrophic Estuary

Recent studies have characterized coastal estuarine systems as important components of the global carbon cycle. This study investigated carbon cycling through the microbial loop of Florida Bay by use of bacterial growth efficiency calculations. Bacterial production, bacterial respiration, and other environmental parameters were measured at three sites located along a historic phosphorus-limitation gradient in Florida Bay and compared to a relatively nutrient enriched site in Biscayne Bay. A new method for measuring bacterial respiration in oligotrophic waters involving tracing respiration of 13C-glucose was developed. The results of the study indicate that 13C tracer assays may provide a better means of measuring bacterial respiration in low nutrient environments than traditional dissolved oxygen consumption-based methods due to strong correlations between incubation length and δ13C values. Results also suggest that overall bacterial growth efficiency may be lower at the most nutrient limited sites.

Application of Numerical Methods to Study Arrangement and Fracture of Lithium-Ion Microstructure

The focus of this work is to develop and employ numerical methods that provide characterization of granular microstructures, dynamic fragmentation of brittle materials, and dynamic fracture of three-dimensional bodies.

We first propose the fabric tensor formalism to describe the structure and evolution of lithium-ion electrode microstructure during the calendaring process. Fabric tensors are directional measures of particulate assemblies based on inter-particle connectivity, relating to the structural and transport properties of the electrode. Applying this technique to X-ray computed tomography of cathode microstructure, we show that fabric tensors capture the evolution of the inter-particle contact distribution and are therefore good measures for the internal state of and electronic transport within the electrode.

We then shift focus to the development and analysis of fracture models within finite element simulations. A difficult problem to characterize in the realm of fracture modeling is that of fragmentation, wherein brittle materials subjected to a uniform tensile loading break apart into a large number of smaller pieces. We explore the effect of numerical precision in the results of dynamic fragmentation simulations using the cohesive element approach on a one-dimensional domain. By introducing random and non-random field variations, we discern that round-off error plays a significant role in establishing a mesh-convergent solution for uniform fragmentation problems. Further, by using differing magnitudes of randomized material properties and mesh discretizations, we find that employing randomness can improve convergence behavior and provide a computational savings.

The Thick Level-Set model is implemented to describe brittle media undergoing dynamic fragmentation as an alternative to the cohesive element approach. This non-local damage model features a level-set function that defines the extent and severity of degradation and uses a length scale to limit the damage gradient. In terms of energy dissipated by fracture and mean fragment size, we find that the proposed model reproduces the rate-dependent observations of analytical approaches, cohesive element simulations, and experimental studies.

Lastly, the Thick Level-Set model is implemented in three dimensions to describe the dynamic failure of brittle media, such as the active material particles in the battery cathode during manufacturing. The proposed model matches expected behavior from physical experiments, analytical approaches, and numerical models, and mesh convergence is established. We find that the use of an asymmetrical damage model to represent tensile damage is important to producing the expected results for brittle fracture problems.

The impact of this work is that designers of lithium-ion battery components can employ the numerical methods presented herein to analyze the evolving electrode microstructure during manufacturing, operational, and extraordinary loadings. This allows for enhanced designs and manufacturing methods that advance the state of battery technology. Further, these numerical tools have applicability in a broad range of fields, from geotechnical analysis to ice-sheet modeling to armor design to hydraulic fracturing.

Evaluating statistical and machine learning methods to predict risk of in-hospital child mortality in Uganda

Thesis (Master's)--University of Washington, 2016-08