Benchmarking spatial joins a la carte

Spatial joins are join operations that involve spatial data types and operators. Spatial access methods are often used to speed up the computation of spatial joins. This paper addresses the issue of benchmarking spatial join operations. For this purpose, we first present a WWW-based benchmark generator to produce sets of rectangles. Using a Web browser, experimenters can specify the number of rectangles in a sample, as well as the statistical distributions of their sizes, shapes, and locations. Second, using the generator and a well-defined set of statistical models we define several tests to compare the performance of three spatial join algorithms: nested loop, scan-and-index, and synchronized tree traversal. We also added two real-life data sets from the Sequoia 2000 storage benchmark. Our results show that the relative performance of the different techniques mainly depends on the selectivity factor of the join predicate. All of the statistical models and algorithms are available on the Web, which allows for easy verification and modification of our experiments.     

Generalized LSG models with spatially varying Coriolis parameter

In this paper, we derive and study approximate balance models for nearly geostrophic shallow water flow where the Coriolis parameter is permitted to vary across the domain as long as it remains nondegenerate. This situation includes, for example, the β-plane approximation to the shallow water equations at mid-latitudes. Our approach is based on changing configuration space coordinates in the underlying variational principle in such a way that consistent asymptotics in the transformed Lagrangian leads to a degenerate Lagrangian structure. In this paper, we restrict our attention to first-order models. We show that the resulting models can be formulated in terms of an advected potential vorticity with a nonlinear vorticity inversion relation. We study the associated solvability conditions and identify a subfamily of models for which these conditions are satisfied without additional restrictions on the data. Finally, we provide the link between our framework and the theory of constrained Hamiltonian systems.     

Effective sea ice area based on a thickness threshold

As Arctic sea ice declines in response to climate change, a shift from thick multiyear ice to a thinner ice cover is occurring. With this transition, ice thicknesses approach a threshold below which ice no longer insulates the atmosphere from oceanic surface fluxes. While this is well known, there are no estimates of the magnitude of this threshold, nor of the proportion of sea ice area that is below this threshold as ice thins. We determine this threshold by simulating the atmospheric response to varying thicknesses, ranging from 0.0 to 2.0 m and determine that threshold to be 0.40–0.50 m. The resulting “effective” ice area is 4–14% lower than reported total ice area, as 0.39–0.97 × 10⁶ km² of the total ice area falls below the threshold throughout the twentieth century, including during notable ice minima. The atmosphere above large non-insulating ice-covered regions is susceptible to more than 2 °C of warming despite ice presence. Observed mean Arctic Ocean ice thickness is projected to fall below this threshold as early as the mid-2020s. Studies on ocean–atmosphere interactions in relation to sea ice area should focus on this insulating sea ice area, where ice is at least 0.40–0.50 m thick, and treat ice regions below 0.40–0.50 m thickness with caution.

     

Aerobic degradation of organic carbon inferred from dinoflagellate cyst decomposition in Southern Ocean sediments

Organic carbon (OC) burial is an important process influencing atmospheric CO₂ concentration and global climate change; therefore it is essential to obtain information on the factors determining its preservation. The Southern Ocean (SO) is believed to play an important role in sequestering CO₂ from the atmosphere via burial of OC. Here we investigate the degradation of organic-walled dinoflagellate cysts (dinocysts) in two short cores from the SO to obtain information on the factors influencing OC preservation. On the basis of the calculated degradation index kt, we conclude that both cores are affected by species-selective aerobic degradation of dinocysts. Further, we calculate a degradation constant k using oxygen exposure time derived from the ages of our cores. The constant k displays a strong relationship with pore-water O₂, suggesting that decomposition of OC is dependent on both the bottom- and pore-water O₂ concentrations.     

Hydrothermal models of the Perth metropolitan area, Western Australia: implications for geothermal energy

Hydrothermal simulations are used to provide insight into the subsurface thermal regime of the Perth metropolitan area (PMA) in Western Australia. High average permeabilities and estimated fluid flow rates in shallow aquifers of the PMA suggest that advection and convection may occur in these aquifers. These processes are simulated, using a new geological model of the PMA to constrain the geometry of aquifers, aquitards and faults. The results show that advection has a strong influence on subsurface temperature, especially in the north of the PMA, where aquifer recharge creates an area of anomalously low temperature. Convection may be important, depending on the permeability of the Yarragadee Aquifer. If convection occurs, it creates thermal highs and lows with a spacing of approximately 5 km. Some of these thermal anomalies migrate over geological time due to coupling between advection and convection, but they are stationary on human timescales. Fault permeability influences the pattern of convection. Advection and convection cause variations in the geothermal gradient which cannot be predicted by conductive models; therefore, these processes should be considered in any model that is used for assessment of geothermal resources in the PMA.     

Flux and dispersion of gases from the “Drachenschlund” hydrothermal vent at 8°18' S, 13°30' W on the Mid-Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH₄ to ³He in the neutrally buoyant plume is quite high, 4 × 10⁸. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol s^− 1 in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH₄ calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of ³He is estimated to be between 0.14 and 1.2 nmol s^− 1. In either case, the ³He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23 yr^− 1, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m²m² s^− 1. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.     

Interactive Visualization to Advance Earthquake Simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth’s surface and interior. Virtual mapping tools allow virtual “field studies” in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method’s strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations.     

Rural Social Networks along Amazonian Rivers: Seeds,Labor and Soccer among Communities on the Napo River,Peru

Research on Amazonian communities has focussed more often on rural‐urban linkages than on links among rural communities. This is unsurprising, given the low density of population, limited intercommunity commerce, and importance of direct city‐market relations. Social relations among rural communities are also important in shaping rural livelihoods and lifeways. We report on the findings of a large‐scale census of communities in the Napo River basin in northeastern Peru (n=174). Data were gathered on intervillage crop seed acquisition and cooperative labor sharing as two key inputs in agriculture, and on intervillage soccer matches, which are integral to rural social life. We analyze the socio‐spatial networks of each practice, paying attention to settlement patterns, community ethnicity, and differential access to the uplands. We find that seeds and labor flow along soccer network lines. Rural social networks appear to be structured strongly by ethnicity (homophily) and reflect important complementarities between upland and lowland communities (weak ties).