Patterns and processes in the California Current System

The California Current System (CCS) is forced by the distribution of atmospheric pressure and associated winds in relation to the west coast of North America. In this paper, we begin with a simplified case of winds and a linear coast, then consider variability characteristic of the CCS, and conclude by considering future change. The CCS extends from the North Pacific Current (～50°N) to off Baja California, Mexico (～15–25°N) with a major discontinuity at Point Conception (34.5°N). Variation in atmospheric pressure affects winds and thus upwelling. Coastal, wind-driven upwelling results in nutrification and biological production and a southward coastal jet. Offshore, curl-driven upwelling results in a spatially large, productive habitat. The California Current flows equatorward and derives from the North Pacific Current and the coastal jet. Dominant modes of spatial and temporal variability in physical processes and biological responses are discussed. High surface production results in deep and bottom waters depleted in oxygen and enriched in carbon dioxide. Fishing has depleted demersal stocks more than pelagic stocks, and marine mammals, including whales, are recovering. Krill, squid, and micronekton are poorly known and merit study. Future climate change will differ from past change and thus prediction of the CCS requires an understanding of its dynamics. Of particular concern are changes in winds, stratification, and ocean chemistry.     

Comparing spatially varying coefficient models: a case study examining violent crime rates and their relationships to alcohol outlets and illegal drug arrests

In this paper, we compare and contrast a Bayesian spatially varying coefficient process (SVCP) model with a geographically weighted regression (GWR) model for the estimation of the potentially spatially varying regression effects of alcohol outlets and illegal drug activity on violent crime in Houston, Texas. In addition, we focus on the inherent coefficient shrinkage properties of the Bayesian SVCP model as a way to address increased coefficient variance that follows from collinearity in GWR models. We outline the advantages of the Bayesian model in terms of reducing inflated coefficient variance, enhanced model flexibility, and more formal measuring of model uncertainty for prediction. We find spatially varying effects for alcohol outlets and drug violations, but the amount of variation depends on the type of model used. For the Bayesian model, this variation is controllable through the amount of prior influence placed on the variance of the coefficients. For example, the spatial pattern of coefficients is similar for the GWR and Bayesian models when a relatively large prior variance is used in the Bayesian model.      

Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

The Use of LiDAR in Digital Surface Modelling: Issues and Errors

Modelling and extracting 3D geographical data presents numerous challenges that require continual research to attempt to evolve an efficient, reliable and accurate solution. LiDAR data capture and analysis has become a preferred acquisition choice for elevation data because the resulting quality and level of detail far exceeds traditional methods for large survey areas. As with any data collection system, LiDAR is prone to errors. Analysing these errors, ascertaining causes and producing error correction strategies is vital if accurate and confident results are to be obtained. Eight years of LiDAR datasets (from 1998 to 2005) have been closely analysed for a large coastal area of South Wales. This article provides a detailed and accurate summary of the identified LiDAR data issues and subsequent errors which affect the accuracy of end products such as Digital Surface Models (DSMs).     

Late glacial fans in the eastern Skagerrak; depositional environment interpreted from swath bathymetry and seismostratigraphy

The origin of acoustically transparent fan deposits overlying glacial till and ice-proximal sediments on the southern margin of the Norwegian Channel has been studied using high-resolution seismic-reflection profiles and multibeam bathymetry. The first deposits overlying glacigenic sediments are a series of stacked, acoustically transparent submarine fans. The lack of glaciomarine sediments below and between individual fans indicates that deposition was rapid and immediately followed the break up of the Late Weichselian ice cover. The fans are overlain by stratified glaciomarine sediments and Holocene mud. Because of the uniformity of this drape, the upper surface of the fan deposits is mimicked at the present seafloor, and the bathymetric images clearly show the spatial relationship of the fans to bedrock ridges and the presence of braided channel-levee systems on the surface of the youngest fans. The acoustically transparent character of the fan deposits indicates that they comprise silt and clay, and their lobate form and lack of internal stratification indicates that they were deposited by debris flows. The channel-levee morphology indicates deposition from more watery hyperconcentrated fluid flows. The fan sediments were either derived from 1) erosion of Mid Weichselian lake deposits in southern Skagerrak or 2) from Late glacial ice-margin lake deposits, ponded against the Norwegian Channel ice stream, which collapsed catastrophically when the lateral support was removed as the ice disintegrated. Fans composed almost exclusively of fine-grained sediment need not, therefore, rule out an origin in a deglacial setting relatively close to the former margins of glaciers and ice sheets.     

High resolution micromill sampling for analysis of fish otoliths by ICP-MS: effects of sampling and specimen preparation on trace element fingerprints

Otoliths are calcified structures in the head of fish that record environmental information about fish's life history. Gathering the elemental information from the core of an otolith corresponding to the juvenile period of fish's life is critical to discriminate the adult fish to their natal habitats reliably. A high resolution micromill has been used to isolate the otolith core from a whole otolith for elemental analysis. The effects of micromilling procedures (e.g., sectioning, embedding and drilling) on contamination to otolith trace element levels were examined using paired blackfin tuna (Thunnus atlanticus) otoliths. Otoliths were decontaminated by dilute hydrogen peroxide and nitric acid throughout to remove surface contamination. A preconcentration procedure was used to determine the trace elements from the small core material by ICP-MS. It was found that micromilling procedures introduce significant contamination to otoliths, especially for Al, Cu, Pb and Zn. The sectioning procedure caused significant contamination for Co and Cu, while the embedding procedure resulted in contamination for nearly all trace elements (Al, Cd, Co, Cu, Ga, Mn, Ni, Pb, V and Zn). The combined sectioning, embedding and drilling procedure also resulted in contamination for most trace elements. Despite the contamination across all procedural steps, the decontamination procedure effectively removed the surface contamination with the exception of Pb and Zn. Bias (e.g., residual contamination) on Pb was small in comparison to overall concentration of Pb expected to occur in fish otoliths, therefore, its effect may be minor in discriminating individuals. Bias on Zn was larger that could limit its application in discriminating individuals.     

Measuring species richness on sandy beach transects: extrapolative estimators and their implications for sampling effort

Species richness is a measure that is fundamental to many studies in ecology, and it is particularly important on sandy beaches, where it underlies patterns described by the broadly accepted swash exclusion hypothesis. However, its estimation in practice is problematic. This has led ecologists in other fields to adopt extrapolative estimators of species richness, which project the total number of species present in a habitat by adjusting upward the number of species observed by an amount related to the number of rare species encountered in the samples. In so doing, the species richness can be estimated, with confidence intervals, at any level of sampling effort. Despite the availability and advantages of these methods, beach ecologists have continued to use the observed species richness as a point estimate of biodiversity for beaches. Here, we employ a Monte Carlo resampling approach over a range of routine transect designs used to sample sandy beaches, and evaluate the performance of seven non-parametric extrapolative estimators for species richness relative to that of the more conventionally used observed species richness. We find that the first-order Jackknife estimator (Jack 1) is the least biased, most accurate and most consistent across sites. Employing this estimator would allow accurate estimation of species richness on short (tens of metres) stretches of beach without exceeding the acceptable levels of sampling effort (4–5 m²). Spreading this effort evenly over three across-shore transects, each with a minimum of 13 equally spaced levels seems appropriately efficient. Although a greater number of research studies is required to ascertain the generality of these results beyond the beaches we sampled, we tentatively recommend the application of our results in biodiversity surveys on sandy beaches.     

Site fidelity and movements of gilthead sea bream (Sparus aurata) in a coastal lagoon (Ria Formosa, Portugal)

Short-time site fidelity and movements of gilthead sea bream (Sparus aurata) in a coastal lagoon were determined using passive acoustic telemetry. Nine fish, ranging from 20.1 to 32.5 cm total length, were surgically implanted with acoustic transmitters and monitored for up to 179 days. Minimum convex polygon areas ranged from 18,698.6 m² to 352,711.9 m². Home range sizes were small, with individuals using core areas on a daily basis. However, these core areas shifted within the study site over time towards the opening to the sea. Two different diel behaviors were recorded, with some individuals more active at night and others during day time. Some individuals also demonstrated homing abilities, returning to the capture site after being released more than 4 km away.