The Use of Geographically Weighted Regression for Spatial Prediction: An Evaluation of Models Using Simulated Data Sets

Increasingly, the geographically weighted regression (GWR) model is being used for spatial prediction rather than for inference. Our study compares GWR as a predictor to (a) its global counterpart of multiple linear regression (MLR); (b) traditional geostatistical models such as ordinary kriging (OK) and universal kriging (UK), with MLR as a mean component; and (c) hybrids, where kriging models are specified with GWR as a mean component. For this purpose, we test the performance of each model on data simulated with differing levels of spatial heterogeneity (with respect to data relationships in the mean process) and spatial autocorrelation (in the residual process). Our results demonstrate that kriging (in a UK form) should be the preferred predictor, reflecting its optimal statistical properties. However the GWR-kriging hybrids perform with merit and, as such, a predictor of this form may provide a worthy alternative to UK for particular (non-stationary relationship) situations when UK models cannot be reliably calibrated. GWR predictors tend to perform more poorly than their more complex GWR-kriging counterparts, but both GWR-based models are useful in that they provide extra information on the spatial processes generating the data that are being predicted.     

The ebb and flow of capital in Indonesian coastal production systems

The global fisheries sector has undergone both rapid industrialization and considerable resource depletion. Unlike fisheries in the Northern Hemisphere, the Indonesian (and indeed Southeast Asian) sector is still largely dominated by small-scale producers, who are partially embedded within a subsistence economy. Changes in the nature of production and livelihoods in the fisheries sector appear similar to those in land-based agriculture but have received far less attention in the literature and demand further analysis given the distinct characteristics of the natural resource base. Using national datasets complemented by insights from a two-month period of fieldwork in South Sulawesi, this paper presents the process of capital intensification underpinning national fisheries growth and how it is transforming small-scale production systems. Despite increasing market integration, we found that smallholders have persisted across coastal production systems to an even stronger degree than land-based agriculture. We suggest some reasons why this is so. However, we also observed evidence of internal class differentiation within coastal communities. Such differentiation, combined with resource degradation and depletion, exposes the poorest in the community to enhanced livelihood vulnerability.     

Transient climate changes in a perturbed parameter ensemble of emissions-driven earth system model simulations

We describe results from a 57-member ensemble of transient climate change simulations, featuring simultaneous perturbations to 54 parameters in the atmosphere, ocean, sulphur cycle and terrestrial ecosystem components of an earth system model (ESM). These emissions-driven simulations are compared against the CMIP3 multi-model ensemble of physical climate system models, used extensively to inform previous assessments of regional climate change, and also against emissions-driven simulations from ESMs contributed to the CMIP5 archive. Members of our earth system perturbed parameter ensemble (ESPPE) are competitive with CMIP3 and CMIP5 models in their simulations of historical climate. In particular, they perform reasonably well in comparison with HadGEM2-ES, a more sophisticated and expensive earth system model contributed to CMIP5. The ESPPE therefore provides a computationally cost-effective tool to explore interactions between earth system processes. In response to a non-intervention emissions scenario, the ESPPE simulates distributions of future regional temperature change characterised by wide ranges, and warm shifts, compared to those of CMIP3 models. These differences partly reflect the uncertain influence of global carbon cycle feedbacks in the ESPPE. In addition, the regional effects of interactions between different earth system feedbacks, particularly involving physical and ecosystem processes, shift and widen the ESPPE spread in normalised patterns of surface temperature and precipitation change in many regions. Significant differences from CMIP3 also arise from the use of parametric perturbations (rather than a multimodel ensemble) to represent model uncertainties, and this is also the case when ESPPE results are compared against parallel emissions-driven simulations from CMIP5 ESMs. When driven by an aggressive mitigation scenario, the ESPPE and HadGEM2-ES reveal significant but uncertain impacts in limiting temperature increases during the second half of the twenty-first century. Emissions-driven simulations create scope for development of errors in properties that were previously prescribed in coupled ocean–atmosphere models, such as historical CO₂ concentrations and vegetation distributions. In this context, historical intra-ensemble variations in the airborne fraction of CO₂ emissions, and in summer soil moisture in northern hemisphere continental regions, are shown to be potentially useful constraints, subject to uncertainties in the relevant observations. Our results suggest that future climate-related risks can be assessed more comprehensively by updating projection methodologies to support formal combination of emissions-driven perturbed parameter and multi-model earth system model simulations with suitable observational constraints. This would provide scenarios underpinned by a more complete representation of the chain of uncertainties from anthropogenic emissions to future climate outcomes.     

Changes to the drivers of fire weather with a warming climate – a case study of southeast Tasmania

Projected changes to the global climate system have great implications for the incidence of large infrequent fires in many regions. Here we examine the synoptic-scale and local-scale influences on the incidence of extreme fire weather days and consider projections of the large-scale mean climate to explore future fire weather projections. We focus on a case study region with periodic extreme fire dangers; southeast Tasmania, Australia. We compare the performance of a dynamically downscaled regional climate model with Global Climate Model outputs as a tool for examining the local-scale influences while accounting for high regional variability. Many of the worst fires in Tasmania and the southeast Australian region are associated with deep cold fronts and strong prefrontal winds. The downscaled simulations reproduce this synoptic type with greater fidelity than a typical global climate model. The incidence of systems in this category is projected to increase through the century under a high emission scenario, driven mainly by an increase in the temperature of air masses, with little change in the strength of the systems. The regional climate model projected increase in frequency is smaller than for the global climate models used as input, with a large model range and natural variability. We also demonstrate how a blocking Foehn effect and topographic channelling contributed to the extreme conditions during an extreme fire weather day in Tasmania in January 2013. Effects such as these are likely to contribute to high fire danger throughout the century. Regional climate models are useful tools that enable various meteorological drivers of fire danger to be considered in projections of future fire danger.     

Trends in Abundance Indices of Fishes in Maryland’s Coastal Bays During 1972–2009

Maryland’s coastal bays provide habitat for juveniles of many commercially and recreationally important species of shellfish and finfish. Since 1972, the Maryland Department of Natural Resources has conducted the Maryland Coastal Bays Trawl and Seine Survey to monitor the populations of key species. The survey has undergone substantial spatial and methodological changes affecting the interpretation of simple indices of abundance. We developed generalized linear models to standardize the indices of abundance of five commonly caught fish species: Atlantic menhaden Brevoortia tyrannus, weakfish Cynoscion regalis, spot Leiostomus xanthurus, bay anchovy Anchoa mitchilli, and summer flounder Paralichthys dentatus. Density declined significantly since 1972 for menhaden, bay anchovy, and spot in at least one region within the coastal bays. The northern bays had significantly higher densities than the southern bays for all species. Changes in abundance indices of the five species examined were not related to sea grass coverage, temperature, salinity, nitrogen-to-phosphorus ratios, and other habitat variables but were likely a result of stock-wide recruitment processes.     

Discovery of possible mega-thrust earthquake along the Seram Trough from records of 1629 tsunami in eastern Indonesian region

Arthur Wichmann’s “Earthquakes of the Indian Archipelago” documents several large earthquakes and tsunami throughout the Banda Arc region that can be interpreted as mega-thrust events. However, the source regions of these events are not known. One of the largest and well-documented events in the catalog is the great earthquake and tsunami affecting the Banda Islands on August 1, 1629. It caused severe damage from a 15-m tsunami that arrived at the Banda Islands about a half hour after violent shaking stopped. The earthquake was also recorded 230 km away in Ambon, but no tsunami is mentioned. This event was followed by at least 9 years of uncommonly frequent seismic activity in the region that tapered off with time, which can be interpreted as aftershocks. The combination of these observations indicates that the earthquake was most likely a mega-thrust event. We use an inverse modeling approach to numerically reconstruct the tsunami, which constrains the likely location and magnitude of the 1629 earthquake. Only, linear numerical models are applied due to the low resolution of bathymetry in the Banda Islands and Ambon. Therefore, we apply various wave amplification factors (1.5–4) derived from simulations of recent, well-constrained tsunami to bracket the upper and lower limits of earthquake moment magnitudes for the event. The closest major earthquake sources to the Banda Islands are the Tanimbar and Seram Troughs of the Banda subduction/collision zone. Other source regions are too far away for such a short arrival time of the tsunami after shaking. Moment magnitudes predicted by the models in order to produce a 15-m tsunami are M_w of 9.8–9.2 on the Tanimbar Trough and M_w 8.8–8.2 on the Seram Trough. The arrival times of these waves are 58 min for Tanimbar Trough and 30 min for Seram Trough. The model also predicts 5-m run-up for Ambon from a Tanimbar Trough source, which is inconsistent with the historical records. Ambon is mostly shielded from a wave generated by a Seram Trough source. We conclude that the most likely source of the 1629 mega-thrust earthquake is the Seram Trough. Only one earthquake >M_w 8.0 is recorded instrumentally from the eastern Indonesia region although high rates of strain (50–80 mm/a) are measured across the Seram section of the Banda subduction zone. Enough strain has already accumulated since the last major historical event to produce an earthquake of similar size to the 1629 event. Due to the rapid population growth in coastal areas in this region, it is imperative that the most vulnerable coastal areas prepare accordingly.     

A shallow seismic refraction study of the Versailles cryptoexplosion structure,central Kentucky

Abstract— The Versailles cryptoexplosion structure, located in central Kentucky, contains a partially brecciated central uplift and a boundary marked by arcuate faulting. Seismic refraction data were used to divide the structure into distinct units on the basis of velocity. A zone of brecciated material shows a decrease in velocity compared with surrounding rocks. Locally thickened strata are present on the flanks of the central uplift. An increase of the bedrock velocity in the area surrounding the structure suggests localized dolomization of fractured bedrock, and a plug of fractured strata near the base of the brecciation appears to be uplifted. Based on the geometry of deformation seen in known impact structures, and the integration of this study with previous geological and geophysical work in the area, it is suggested that our results support the interpretation of the Versailles structure as an eroded astrobleme (Black, 1964b; Seeger, 1972).     

underwThe Hawaii-2 Observatory seismic system

The Hawaii-2 Observatory seismic system is currently transmitting high-quality seismic data from the ocean floor in the central NE Pacific Ocean through Hawaii to the IRIS Data Management Center. The system includes broad-band seismic, geophone, acoustic, and ocean current sensors. The seismic sensors are buried about 0.4 m below the ocean floor to improve coupling to the ocean bottom and to reduce noise levels. The system can be remotely calibrated, leveled and locked, and gains can be changed on command from shore. Data are temporarily stored in the seismic package for retransmission as needed to correct for transmission problems and to prevent loss of data. Data generated are valuable for studies of the Earth's structure and the dynamics of earthquakes     

Shaking table tests of a resilient bridge system with precast reinforced concrete columns equipped with springs

This paper presents the shake table test results of a novel system for the design of precast reinforced concrete bridges. The specimen comprises a slab and four precast columns. The connections are dry and the columns are connected to the slab by an ungrouted tendon. One of the tendon ends is anchored above the slab, in series with a stack of washer springs, while the other end is anchored at the bottom of the column. The addition of such a flexible restraining system increases the stability of the system, while keeping it relatively flexible allowing it to experience negative post-uplift stiffness. It is a form of seismic isolation. Anchoring the tendon within the column, caps the design moment of the foundation, and reduces its size. One hundred and eighty-one shake table tests were performed. The first 180 caused negligible damage to the specimen, mainly abrasion at the perimeter of the column top ends. Hence, the system proved resilient. The 181^st excitation caused collapse, because the tendons unexpectedly failed at a load less than 50% of their capacity (provided by the manufacturer), due to the failure of their end socket. This highlights the importance of properly designing the tendons. The tests were used to statistically validate a rigid body model. The model performed reasonably well never underestimating the median displacement response of the center of mass of the slab by more than 30%. However, the model cannot predict the torsion rotation of the slab that was observed in the tests and is due to imperfections.