Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

Mechanistic Models of Anti-HIV Microbicide Drug Delivery

A new modality for preventing HIV transmission is emerging in the form of topical microbicides. Some clinical trials have shown some promising results of these methods of protection while other trials have failed to show efficacy. Due to the relatively novel nature of microbicide drug transport, a rigorous, deterministic analysis of that transport can help improve the design of microbicide vehicles and understand results from clinical trials. This type of analysis can aid microbicide product design by helping understand and organize the determinants of drug transport and the potential efficacies of candidate microbicide products.

Microbicide drug transport is modeled as a diffusion process with convection and reaction effects in appropriate compartments. This is applied here to vaginal gels and rings and a rectal enema, all delivering the microbicide drug Tenofovir. Although the focus here is on Tenofovir, the methods established in this dissertation can readily be adapted to other drugs, given knowledge of their physical and chemical properties, such as the diffusion coefficient, partition coefficient, and reaction kinetics. Other dosage forms such as tablets and fiber meshes can also be modeled using the perspective and methods developed here.

The analyses here include convective details of intravaginal flows by both ambient fluid and spreading gels with different rheological properties and applied volumes. These are input to the overall conservation equations for drug mass transport in different compartments. The results are Tenofovir concentration distributions in time and space for a variety of microbicide products and conditions. The Tenofovir concentrations in the vaginal and rectal mucosal stroma are converted, via a coupled reaction equation, to concentrations of Tenofovir diphosphate, which is the active form of the drug that functions as a reverse transcriptase inhibitor against HIV. Key model outputs are related to concentrations measured in experimental pharmacokinetic (PK) studies, e.g. concentrations in biopsies and blood. A new measure of microbicide prophylactic functionality, the Percent Protected, is calculated. This is the time dependent volume of the entire stroma (and thus fraction of host cells therein) in which Tenofovir diphosphate concentrations equal or exceed a target prophylactic value, e.g. an EC50.

Results show the prophylactic potentials of the studied microbicide vehicles against HIV infections. Key design parameters for each are addressed in application of the models. For a vaginal gel, fast spreading at small volume is more effective than slower spreading at high volume. Vaginal rings are shown to be most effective if inserted and retained as close to the fornix as possible. Because of the long half-life of Tenofovir diphosphate, temporary removal of the vaginal ring (after achieving steady state) for up to 24h does not appreciably diminish Percent Protected. However, full steady state (for the entire stromal volume) is not achieved until several days after ring insertion. Delivery of Tenofovir to the rectal mucosa by an enema is dominated by surface area of coated mucosa and whether the interiors of rectal crypts are filled with the enema fluid. For the enema 100% Percent Protected is achieved much more rapidly than for vaginal products, primarily because of the much thinner epithelial layer of the mucosa. For example, 100% Percent Protected can be achieved with a one minute enema application, and 15 minute wait time.

Results of these models have good agreement with experimental pharmacokinetic data, in animals and clinical trials. They also improve upon traditional, empirical PK modeling, and this is illustrated here. Our deterministic approach can inform design of sampling in clinical trials by indicating time periods during which significant changes in drug concentrations occur in different compartments. More fundamentally, the work here helps delineate the determinants of microbicide drug delivery. This information can be the key to improved, rational design of microbicide products and their dosage regimens.

Orbitally tuned time scale and astronomical forcing in the middle Eocene to early Oligocene determined on ODP and IODP sediment cores

Deciphering the driving mechanisms of Earth system processes, including the climate dynamics expressed as paleoceanographic events, requires a complete, continuous, and high-resolution stratigraphy that is very accurately dated. In this study, we construct a robust astronomically calibrated age model for the middle Eocene to early Oligocene interval (31-43 Ma) in order to permit more detailed study of the exceptional climatic events that occurred during this time, including the Middle Eocene Climate Optimum and the Eocene/Oligocene transition. A goal of this effort is to accurately date the middle Eocene to early Oligocene composite section cored during the Pacific Equatorial Age Transect (PEAT, IODP Exp. 320/321). The stratigraphic framework for the new time scale is based on the identification of the stable long eccentricity cycle in published and new high-resolution records encompassing bulk and benthic stable isotope, calibrated XRF core scanning, and magnetostratigraphic data from ODP Sites 171B-1052, 189-1172, 199-1218, and 207-1260 as well as IODP Sites 320-U1333, and -U1334 spanning magnetic polarity Chrons C12n to C20n. Subsequently we applied orbital tuning of the records to the La2011 orbital solution. The resulting new time scale revises and refines the existing orbitally tuned age model and the Geomagnetic Polarity Time Scale from 31 to 43 Ma. Our newly defined absolute age for the Eocene/Oligocene boundary validates the astronomical tuned age of 33.89 Ma identified at the Massignano (Italy) global stratotype section and point. Our compilation of geochemical records of climate-controlled variability in sedimentation through the middle-to-late Eocene and early Oligocene demonstrates strong power in the eccentricity band that is readily tuned to the latest astronomical solution. Obliquity driven cyclicity is only apparent during very long eccentricity cycle minima around 35.5 Ma, 38.3 Ma and 40.1 Ma.

(Table 1) Mn, Fe, TOC and CaCO3 concentrations in bottom sediments and silt waters of Kandalaksha Bay, White Sea

The redox stratification of bottom sediments in Kandalaksha Bay, White Sea, is characterized by elevated concentrations of Mn (3-5%) and Fe (7.5%) in the uppermost layer, which is two orders of magnitude and one and a half times, respectively, higher than the average concentrations of these elements in the Earth's crust. The high concentrations of organic matter (Corg = 1-2%) in these sediments cannot maintain (because of its low reaction activity) the sulfate-reducing process (the concentration of sulfide Fe is no higher than 0.6%). The clearest manifestation of diagenesis is the extremely high Mn2+ concentration in the silt water (>500 µM), which causes its flux into the bottom water, oxidation in contact with oxygen, and the synthesis of MnO2 oxy-hydroxide enriching the surface layer of the sediments. Such migrations are much less typical of Fe. Upon oxygen exhaustion in the uppermost layer of the sediments, the synthesized oxyhydroxides (MnO2 and FeOOH) serve as oxidizers of organic matter during anaerobic diagenesis. The calculated diffusion-driven Mn flux from the sediments (280 µM/m**2 day) and corresponding amount of forming Mn oxyhydrate as compared to opposite oxygen flux to sediments (1-10 mM/m**2 day) indicates that >10% organic matter in the surface layer of the sediments can be oxidized with the participation of MnO2. The roles of other oxidizers of organic matter (FeOOH and SO4**2-) becomes discernible at deeper levels of the sediments. The detailed calculation of the balance of reducing processes testifies to the higher consumption of organic matter during the diagenesis of surface sediments than it follows from the direct determination of Corg. The most active diagenetic redox processes terminate at depths of 25-50 cm. Layers enriched in Mn at deeper levels are metastable relicts of its surface accumulation and are prone to gradual dissemination.

The Crisis in the UK's High Streets: can the evening and night-time economy help?

Changes in retailing over the last half century have had a detrimental effect on the UK’s local high streets. The recession of 2008-2012 exacerbated these trends leading to a high number of vacancies and neglected properties. The impact was sufficiently severe for the term ‘crisis’ to be used in connection with the British high street. In the academic and commercial reports generated by the recognition that the high street needed to adapt to changing circumstances, a view emerged that the leisure component of high street activity would gain in importance. This article reviews the relationship of the evening and nighttime economy to the high street and considers its potential in reinventing the vitality that is normally associated with these mixed-use urban corridors. The article argues that there is hope in the high street offering a different type of experience to the mainstream forms of entertainment that are consolidating in major town and city centres. It concludes by suggesting that for this to be successful, some public support is necessary.

Development of computer application for determination of chemical diffusion properties in biological tissues

The main objective of this work was to develop an application capable of determining the diffusion times and diffusion coefficients of optical clearing agents and water inside a known type of muscle. Different types of chemical agents can also be used with the method implemented, such as medications or metabolic products. Since the diffusion times can be calculated, it is possible to describe the dehydration mechanism that occurs in the muscle. The calculation of the diffusion time of an optical clearing agent allows to characterize the refractive index matching mechanism of optical clearing. By using both the diffusion times and diffusion of water and clearing agents not only the optical clearing mechanisms are characterized, but also information about optical clearing effect duration and magnitude is obtained. Such information is crucial to plan a clinical intervention in cooperation with optical clearing. The experimental method and equations implemented in the developed application are described in throughout this document, demonstrating its effectiveness. The application was developed in MATLAB code, but the method was personalized so it better fits the application needs. This process significantly improved the processing efficiency, reduced the time to obtain he results, multiple validations prevents common errors and some extra functionalities were added such as saving application progress or export information in different formats. Tests were made using glucose measurements in muscle. Some of the data, for testing purposes, was also intentionally changed in order to obtain different simulations and results from the application. The entire project was validated by comparing the calculated results with the ones found in literature, which are also described in this document.

In situ differentiation of labile/inert metal species in Brazilian tropical rivers by means of a time-controlled batch-procedure based on TEPHA resin

The present study deals with a new analytical procedure based on a cellulose diffusion membrane and immobilised tetraethylene-pentamine-hexaacetate chelator (DM-TEPHA) for an in situ differentiation of labile and inert metal species in aquatic systems. The DM-TEPHA system was prepared by placing TEPHA chelator in pre-purified cellulose bags and in situ applied immersing the system in two Brazilian rivers to study the relative lability of metal species (Cu, Pb, Fe, Mn and Ni) as a function of the time and the quantity of exchanger, respectively. The procedure is simple and enables a new perspective for understanding the complexation, transport, stability and lability of metal species in aquatic systems rich in organic matter.

The Predation of Artemia Nauplii by the Larvae of the Amazon River Prawn, Macrobrachium amazonicum (Heller, 1862), is Affected by Prey Density, Time of Day, and Ontogenetic Development

The present study analyzed the effects of prey density, the time of day, and ontogenetic development on the predation of Artemia nauplii by the larvae of the Amazon river prawn, Macrobrachium amazonicum, as well as possible synergy among these factors. Larvae were raised in 120-L tanks with biological filter systems, and fed on recently hatched Artemia nauplii, using two feeding management protocols: (a) fed once per day at 2000 h (high density HD) and (b) half of the ration provided at 2000 h, complemented at 0800 h the following day by a replacement of the nauplii consumed up to a maximum of the full ration (low density with replacement LDWR). Each treatment consisted of six replicates. The consumption of nauplii was estimated prior to the feeding times. Consumption varied according to time of day, ontogenetic development, and feeding protocol. The larvae ingested more nauplii during the daytime at most developmental stages. Ingestion rates were similar during the day under both treatments, but at night the higher density of prey in the HD treatment caused a higher encounter rate and increased ingestion of nauplii by the larvae. Among the performance indicators only survival was greater in HD in comparison with LDWR; productivity and dry weight were similar. The results indicate a circadian trophic rhythm in M. amazonicum, with the encounter rate being an important mechanism for the capture of prey during the night. A second mechanism probably the visual system aids the perception of prey during the daytime. Based on these results, we suggest that feeding captive Amazon river prawn larvae only once a day would be appropriate and economically beneficial. Further work is necessary to determine the most effective time that this single feed should be applied.

Nano-Engineering of Densely Packed Electrochemical Energy Storage Architectures by Atomic Layer Deposition

Nanostructures are highly attractive for future electrical energy storage devices because they enable large surface area and short ion transport time through thin electrode layers for high power devices. Significant enhancement in power density of batteries has been achieved by nano-engineered structures, particularly anode and cathode nanostructures spatially separated far apart by a porous membrane and/or a defined electrolyte region. A self-aligned nanostructured battery fully confined within a single nanopore presents a powerful platform to determine the rate performance and cyclability limits of nanostructured storage devices. Atomic layer deposition (ALD) has enabled us to create and evaluate such structures, comprised of nanotubular electrodes and electrolyte confined within anodic aluminum oxide (AAO) nanopores. The V2O5- V2O5 symmetric nanopore battery displays exceptional power-energy performance and cyclability when tested as a massively parallel device (~2billion/cm2), each with ~1m3 volume (~1fL). Cycled between 0.2V and 1.8V, this full cell has capacity retention of 95% at 5C rate and 46% at 150C, with more than 1000 charge/discharge cycles. These results demonstrate the promise of ultrasmall, self-aligned/regular, densely packed nanobattery structures as a testbed to study ionics and electrodics at the nanoscale with various geometrical modifications and as a building block for high performance energy storage systems[1, 2]. Further increase of full cell output potential is also demonstrated in asymmetric full cell configurations with various low voltage anode materials. The asymmetric full cell nanopore batteries, comprised of V2O5 as cathode and prelithiated SnO2 or anatase phase TiO2 as anode, with integrated nanotubular metal current collectors underneath each nanotubular storage electrode, also enabled by ALD. By controlling the amount of lithium ion prelithiated into SnO2 anode, we can tune full cell output voltage in the range of 0.3V and 3V. This asymmetric nanopore battery array displays exceptional rate performance and cyclability. When cycled between 1V and 3V, it has capacity retention of approximately 73% at 200C rate compared to 1C, with only 2% capacity loss after more than 500 charge/discharge cycles. With increased full cell output potential, the asymmetric V2O5-SnO2 nanopore battery shows significantly improved energy and power density. This configuration presents a more realistic test - through its asymmetric (vs symmetric) configuration – of performance and cyclability in nanoconfined environment. This dissertation covers (1) Ultra small electrochemical storage platform design and fabrication, (2) Electron and ion transport in nanostructured electrodes inside a half cell configuration, (3) Ion transport between anode and cathode in confined nanochannels in symmetric full cells, (4) Scale up energy and power density with geometry optimization and low voltage anode materials in asymmetric full cell configurations. As a supplement, selective growth of ALD to improve graphene conductance will also be discussed[3]. References: 1. Liu, C., et al., (Invited) A Rational Design for Batteries at Nanoscale by Atomic Layer Deposition. ECS Transactions, 2015. 69(7): p. 23-30. 2. Liu, C.Y., et al., An all-in-one nanopore battery array. Nature Nanotechnology, 2014. 9(12): p. 1031-1039. 3. Liu, C., et al., Improving Graphene Conductivity through Selective Atomic Layer Deposition. ECS Transactions, 2015. 69(7): p. 133-138.

ANALYSIS OF MACROINVERTEBRATE COMMUNITIES IN SEASONAL WETLANDS THROUGH TIME, ACROSS SPACE, AND USING SPECIES TRAITS

Restoration of natural wetlands may be informed by macroinvertebrate community composition. Macroinvertebrate communities of wetlands are influenced by environmental characteristics such as vegetation, soil, hydrology, land use, and isolation. This dissertation explores multiple approaches to the assessment of wetland macroinvertebrate community composition, and demonstrates how these approaches can provide complementary insights into the community ecology of aquatic macroinvertebrates. Specifically, this work focuses on macroinvertebrates of Delmarva Bays, isolated seasonal wetlands found on Maryland’s eastern shore. A comparison of macroinvertebrate community change over a nine years in a restored wetland complex indicated that the macroinvertebrate community of a rehabilitated wetlands more rapidly approximated the community of a reference site than did a newly created wetland. The recovery of a natural macroinvertebrate community in the rehabilitated wetland indicated that wetland rehabilitation should be prioritized over wetland creation and long-term monitoring may be needed to evaluate restoration success. This study also indicated that characteristics of wetland vegetation reflected community composition. The connection between wetland vegetation and macroinvertebrate community composition led to a regional assessment of predaceous diving beetle (Coleoptera: Dytiscidae) community composition in 20 seasonal wetlands, half with and half without sphagnum moss (Sphagnum spp.). Species-level identifications indicated that wetlands with sphagnum support unique and diverse assemblages of beetles. These patterns suggest that sphagnum wetlands provide habitat that supports biodiversity on the Delmarva Peninsula. To compare traits of co-occurring beetles, mandible morphology and temporal and spatial variation were measured between three species of predaceous diving beetles. Based on mandible architecture, all species may consume similarly sized prey, but prey characteristics likely differ in terms of piercing force required for successful capture and consumption. Therefore, different assemblages of aquatic beetles may have different effects on macroinvertebrate community structure. Integrating community-level and species-level data strengthens the association between individual organisms and their ecological role. Effective restoration of imperiled wetlands benefits from this integration, as it informs the management practices that both preserve biodiversity and promote ecosystem services.