A systematic review of strong gravitational lens modeling software

Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This review of 13 software packages was undertaken to provide a single source of information. Gravitational lens models are classified as parametric models or non-parametric models, and can be further divided into research and educational software. Software used in research includes the GRAVLENS package (with both gravlens and lensmodel), Lenstool, LensPerfect, glafic, PixeLens, SimpLens, Lensview, and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing, lens and MOWGLI are categorized as educational programs that are useful for demonstrating various aspects of lensing. Each of the 13 software packages is reviewed with regard to software features (installation, documentation, files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a brief review of studies where they have been used. Recent studies have demonstrated the utility of strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of strong gravitational lens systems are needed. This is the first systematic review of strong gravitational lens modeling software, providing investigators with a starting point for future software development to further advance gravitational lens modeling research.     

Stratification and mixing of a post-glacial Neoproterozoic ocean: Evidence from carbon and sulfur isotopes in a cap dolostone from northwest China

To improve our knowledge about the geochemical and environmental aftermath of Neoproterozoic global glaciations, we analyzed stable isotopes (δ¹³C, δ¹⁸O, δ³⁴S) and elemental concentrations (Ca, Mg, S, Sr, Fe, and Mn) of the ～ 10-m-thick Zhamoketi cap dolostone atop the Tereeken diamictite in the Quruqtagh area, eastern Chinese Tianshan. Available chemostratigraphic data suggest that the Tereeken diamictite is probably equivalent to the Marinoan glaciation. Our new data indicate that organic and carbonate carbon isotopes of the Zhamoketi cap dolostone show little stratigraphic variations, averaging ? 28.2‰ and ? 4.6‰, respectively. In contrast, sulfur isotopes show significant stratigraphic variations. Carbonate associated sulfate (CAS) abundance decreases rapidly in the basal cap dolostone and δ³⁴S_CAS composition varies between + 9‰ and + 15‰ in the lower 2.5 m. In the overlying interval, CAS abundance remains low while δ³⁴S_CAS rises ～ 5‰ and varies more widely between + 10‰ and + 21‰. The range of δ³⁴S_py of the cap dolostone overlaps with that of δ³⁴S_CAS, but direct comparison shows that δ³⁴S_py is typically greater than δ³⁴S_CAS measured from the same samples. Hypotheses to explain the observations must account for both the remarkable sulfur isotope enrichment of pyrites and the inverse fractionation. We propose that CAS and pyrite were derived from two isotopically distinct reservoirs in a chemically stratified basin or a basin with a sulfate minimum zone. In this model, CAS was derived from shallow, oxic surface waters with moderate sulfate concentration and depleted in ³⁴S due to the post-glacial influx of sulfur from continental weathering. In contrast, pyrite was derived from anoxic bottom waters (or a sulfate minimum zone) with low sulfate concentration and ³⁴S enrichment due to long-term syn-glacial sulfate reduction. The rapid shift in CAS abundance and sulfur isotope composition within the cap dolostone is interpreted to reflect the mixing of the two reservoirs after initial deglaciation. Comparison with other post-Marinoan cap carbonates shows significant spatial heterogeneity in δ³⁴S_CAS, which together with strong temporal variation in δ³⁴S_CAS, points to generally low sulfate concentrations in post-Marinoan oceans.     

Groundwater discharge along a channelized Coastal Plain stream

In the Coastal Plain of the southeastern USA, streams have commonly been artificially channelized for flood control and agricultural drainage. However, groundwater discharge along such streams has received relatively little attention. Using a combination of stream- and spring-flow measurements, spring temperature measurements, temperature profiling along the stream-bed, and geologic mapping, we delineated zones of diffuse and focused discharge along Little Bayou Creek, a channelized, first-order perennial stream in western Kentucky. Seasonal variability in groundwater discharge mimics hydraulic-head fluctuations in a nearby monitoring well and spring-discharge fluctuations elsewhere in the region, and is likely to reflect seasonal variability in recharge. Diffuse discharge occurs where the stream is incised into the semi-confined regional gravel aquifer, which is comprised of the Mounds Gravel. Focused discharge occurs upstream where the channel appears to have intersected preferential pathways within the confining unit. Seasonal fluctuations in discharge from individual springs are repressed where piping results in bank collapse. Thereby, focused discharge can contribute to the morphological evolution of the stream channel.     

Large-scale cold water dispersant effectiveness experiments with Alaskan crude oils and Corexit 9500 and 9527 dispersants

There continues to be reluctance in some jurisdictions to use chemical dispersants as a viable countermeasure for accidental oil spills. One argument used by some opponents to dispersant use is that “chemical dispersants do not work effectively in cold water”. To address this issue, the U.S. Minerals Management Service (MMS) funded and conducted two series of large-scale dispersant experiments in very cold water at Ohmsett - The National Oil Spill Response Test Facility, located in Leonardo, New Jersey in February-March 2006 and January-March 2007. Alaska North Slope, Endicott, Northstar and Pt. McIntyre crude oils and Corexit 9500 and Corexit 9527 dispersants were used in the two test series. The crude oils were tested both when fresh and after weathering. Results demonstrated that both Corexit 9500 and Corexit 9527 dispersants were 85-99% effective in dispersing the fresh and weathered crude oils tested at cold temperatures. The MMS expects that results from these test series will assist government regulators and responders in making science based decisions on the use of dispersants as a response tool for oil spills in the Arctic.     

The use of monitoring data for identifying factors influencing phytoplankton bloom dynamics in the eutrophic Taw Estuary, SW England

Using the Taw Estuary as an example, data routinely collected by the Environment Agency for England and Wales over the period 1990-2004 were interrogated to identify the drivers of excessive algal growth. The estuary was highly productive with chlorophyll concentrations regularly exceeding 100 μg L⁻¹, mostly during periods of low freshwater input from the River Taw when estuarine water residence times were longest. However, algal growth in mid estuary was often inhibited by ammonia inputs from the adjacent sewage treatment works. The reported approach demonstrates the value of applying conventional statistical analyses in a structured way to existing monitoring data and is recommended as a useful tool for the rapid assessment of eutrophication. However, future estuarine monitoring should include the collection of dissolved organic nutrient data and targeted high temporal resolution data because the drivers of eutrophication are complex and often very specific to a particular estuary.     

Effect of coherent noise on single-station direction of arrival estimation

Polarization analysis of multi-component seismic data is used in both exploration seismology and earthquake seismology. In single-station polarization processing, it is generally assumed that any noise present in the window of analysis is incoherent, i.e., does not correlate between components. This assumption is often violated in practice because several overlapping seismic events may be present in the data. The additional arrival(s) to that of interest can be viewed as coherent noise. This paper quantifies the error because of coherent noise interference. We first give a general theoretical analysis of the problem. A simple mathematical wavelet is then used to obtain a closed-form solution to the principal direction estimated for a transient incident signal superposed with a time-shifted, unequal amplitude version of itself, arriving at an arbitrary angle to the first wavelet. The effects of relative amplitude, arrival angle, and the time delay of the two wavelets on directional estimates are investigated. Even for small differences in angle of arrival, there may be significant error (>10°) in the azimuth estimate.     

Quantifying and predicting the accuracy of radar-based quantitative precipitation forecasts

On the basis that hydrological users need to know the forecast uncertainty at the time that the forecast is issued, we computed distributions of radar rainfall forecast uncertainty as a function of forecast lead time, basin size, and forecasted rainfall intensity using data from the US 3-D National Mosaic of radar data. We document how exceptional forecasts such as those of heavy rainfall are generally biased. Since forecast uncertainty is also weather dependent, we tried to find good predictors to help either reduce the forecast uncertainty or better define it. These predictors were based either on characteristics of the current precipitation field or on the performance of the nowcast in the immediate past. The value of some predictors, especially those based on the properties of large-scale rainfall patterns, was significant though modest, the predictors being generally more skillful at characterizing forecast uncertainty than at improving forecast accuracy.     

A parallel mesh-free contaminant transport model based on the Analytic Element and Streamline Methods

This paper introduces a new method for simulating large-scale subsurface contaminant transport that combines an Analytic Element Method (AEM) groundwater flow solution with a split-operator Streamline Method for modeling reactive transport. The key feature of the method is the manner in which the vertically integrated AEM flow solution is used to construct three-dimensional particle tracks that define the geometry of the Streamline Method. The inherently parallel nature of the algorithm supports the development of reactive transport models for spatial domains much larger than current grid-based methods. The applicability of the new approach is verified for cases with negligible transverse dispersion through comparisons to analytic solutions and existing numerical solutions, and parallel performance is demonstrated through a realistic test problem based on the regional-scale transport of agricultural contaminants from spatially distributed sources.