Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.     

A method for determination of epicenters of near earthquakes

Summary A method is developed, which permits a rapid and reliable determination of the location of near epicenters (within about 10^o) using the differences of arrival times of a common phase (usually Sg) to a net of three stations. The resulting equation for computation of epicentral distances is solved numerically once for all by an electronic computer for every station triangle and for all possible time differences for every wave used. The calculated epicentral distances and coordinates are arranged in tables with the two travel-time differences for each triangle as entries. Such tables permit an immediate solution of the problem in every special case. The accuracy of the method is investigated, and necessary precautions in its use are emphasized. Application is made to a few cases of Sg recorded at the Swedish stations.

---

Zusammenfassung Eine Methode wird entwickelt, die eine schnelle und zuverlässige Bestimmung der Lage von nahegelegenen Epizentren (innerhalb rund 10^o) mittels der Laufzeitdifferenzen einer gegebenen Welle (gewöhnlichSg) in einem Netz von drei Stationen ermöglicht. Die erhaltene Gleichung für die Berechnung der Herdentfernungen wird mittels einer elektronischen Rechenmaschine für jedes Stations-dreieck und für alle möglichen Zeitdifferenzen für jede verwendete Welle ein für allemal numerisch gelöst. Die berechneten Herdentfernungen und-Koordinaten werden in Tabellenform dargestellt mit den beiden Laufzeitdifferenzen für jedes Dreieck als Eingangswerte. Solche Tabellen ermöglichen eine unmittelbare Lösung des Problems in jedem einzelnen Fall. Die Genauigkeit der Methode wird untersucht, und die notwendigen Vorsichtsmassregeln bei ihrer Verwendung werden hervorgehoben. Die Methode wird auf einige Registrierungen vonSg an den schwedischen Stationen angewendet.

     

Correlation between regional and global seismic activity

The regional seismic activity in Sweden is significantly correlated to the global seismic activity. This result suggests that the regional Swedish activity is linked to the global activity by intraplate tectonics, rather than being caused exclusively by regional effects, such as the land uplift after the last glaciation.     

The connection between Labrador Sea buoyancy loss, deep western boundary current strength, and Gulf Stream path in an ocean circulation model

The sensitivity of the North Atlantic gyre circulation to high latitude buoyancy forcing is explored in a global, non-eddy resolving ocean general circulation model. Increased buoyancy forcing strengthens the deep western boundary current, the northern recirculation gyre, and the North Atlantic Current, which leads to a more realistic Gulf Stream path. High latitude density fluxes and surface water mass transformation are strongly dependent on the choice of sea ice and salinity restoring boundary conditions. Coupling the ocean model to a prognostic sea ice model results in much greater buoyancy loss in the Labrador Sea compared to simulations in which the ocean is forced by prescribed sea ice boundary conditions. A comparison of bulk flux forced hindcast simulations which differ only in their sea ice and salinity restoring forcings reveals the effects of a mixed thermohaline boundary condition transport feedback whereby small, positive temperature and salinity anomalies in subpolar regions are amplified when the gyre spins up as a result of increased buoyancy loss and convection. The primary buoyancy flux effects of the sea ice which cause the simulations to diverge are ice melt, which is less physical in the diagnostic sea ice model, and insulation of the ocean, which is less physical with the prognostic sea ice model. Increased salinity restoring ensures a more realistic net winter buoyancy loss in the Labrador Sea, but it is found that improvements in the Gulf Stream simulation can only be achieved with the excessive buoyancy loss associated with weak salinity restoring.     

Eddy formation in the bays of Kamchatka and fluxes to the open ocean

Ocean Dynamics - The Eastern Kamchatka Current (EKC) is the western boundary current of the North Pacific subpolar gyre. Southeast of the Kamchatka Peninsula lies a large anticyclonic eddy, the...     

Using horizontal and vertical building structure to constrain indirect sensor orientation

This paper presents a method to integrate linear horizontal, vertical and right-angled scene structures into the bundle adjustment of image sequences. An increasing number of airborne image acquisition systems is available and equipped with non-metric small- or medium-frame cameras and no or insufficiently accurate INS devices. In cases where the data is to be used for the production of geo-spatial data, where a certain accuracy and precision is required, an indirect sensor orientation, possibly including self-calibration, needs to be performed. The idea which led to the presented approach is to reduce the number of GCPs necessary for this task by applying the mentioned scene structures. The method directly uses the linear structures, visible at man-made objects as fictive observations within the adjustment, while self-calibration of intrinsic camera parameters and lens distortion is included as well.Experiments with two datasets demonstrate that, through this method, only limited GCP information is required to obtain satisfactory results. In fact, in one experiment using oblique images, several scene constraints were provided and only the datum was defined by ground control. The residuals at check points from this setup were similar to the traditional case where several well-distributed GCPs were available in the scene. In the second experiment the ability of this approach to support the bundle adjustment was shown for a UAV dataset. Although no GCP and camera calibration information was available, the visual inspection of adjusted object points and the residuals at horizontal structures confirmed the ability of the method to align an image block with the structure, as embodied in the defined scene constraints. Despite the convincing outcome of the experiments, it needs to be mentioned that some manual work is still involved in defining the constraints. In future work the issue of automation will be addressed.     

Benthic photosynthesis in submerged Wadden Sea intertidal flats

In this study we compare benthic photosynthesis during inundation in coarse sand, fine sand, and mixed sediment (sand/mud) intertidal flats in the German Wadden Sea. In situ determinations of oxygen-, DIC- and nutrient fluxes in stirred benthic chamber incubations were combined with measurements of sedimentary chlorophyll, incident light intensity at the sediment surface and scalar irradiance within the sediment. During submergence, microphytobenthos was light limited at all study sites as indicated by rapid response of gross photosynthesis to increasing incident light at the sea floor. However, depth integrated scalar irradiance was 2 to 3 times higher in the sands than in the mud. Consequently, gross photosynthesis in the net autotrophic fine sand and coarse sand flats during inundation was on average 4 and 11 times higher than in the net heterotrophic mud flat, despite higher total chlorophyll concentration in mud. Benthic photosynthesis may be enhanced in intertidal sands during inundation due to: (1) higher light availability to the microphytobenthos in the sands compared to muds, (2) more efficient transport of photosynthesis-limiting solutes to the microalgae with pore water flows in the permeable sands, and (3) more active metabolic state and different life strategies of microphytobenthos inhabiting sands.     

Flexible management of fishing rights and a sustainable fisheries industry in Europe

Although overexploitation of commercial fish stocks in European waters has been in the public debate now for more than 20 years, the European Union has so far failed to implement sustainable fisheries management. Millions in subsidies paid to the fishing industry have led to significant excess capacity in the fishing fleet. Various feeble attempts to stop overexploitation of marine resources have failed. The cause is that fishing policy is highly dominated by short-term socioeconomic interests. There is an urgent need for a new fisheries management system in Europe that supports reductions in the fishing fleet, increases responsibility among fishers and guarantees long-term conservation of natural marine resources.Transferable rights to fish have proved a reliable and effective means of creating incentives to conserve marine resources. By strengthening individual fishing rights under flexible quota management systems, the EU Member States could, within the Common Fisheries Policy, make a significant contribution to conserving fish stocks, to reducing excess capacity and to raising the profitability of the fisheries industry. A closer look at existing reservations against a flexible management system shows most of the objections to be overstated or capable of resolution.     

Frequency distributions and correlations of solar X-ray flare parameters

We have determined frequency distributions of flare parameters from over 12000 solar flares recorded with the Hard X-Ray Burst Spectrometer (HXRBS) on the Solar Maximum Mission (SMM) satellite. These parameters include the flare duration, the peak counting rate, the peak hard X-ray flux, the total energy in electrons, and the peak energy flux in electrons (the latter two computed assuming a thick-target flare model). The energies were computed above a threshold energy between 25 and 50 keV. All of the distributions can be represented by power laws above the HXRBS sensitivity threshold. Correlations among these parameters are determined from linear regression fits as well as from the slopes of the frequency distributions. Variations of the frequency distributions were investigated with respect to the solar activity cycle.Theoretical models for the frequency distribution of flare parameters depend on the probability of flaring and the temporal evolution of the flare energy build-up. Our results are consistent with stochastic flaring and exponential energy build-up, with an average build-up time constant that is 0.5 times the mean time between flares. The measured distributions of flares are also consistent with predicted distributions of flares from computer simulations of avalanche models that are governed by the principle of self-organized criticality.     

Sedimentary evolution of a late Pleistocene wetland indicating extreme coastal uplift in southern Tanzania

Facies analyses of Pleistocene deposits from southern coastal Tanzania (Lindi District) document that sediments formed in a wetland evolving on a coastal terrace in the Lindi Fracture Zone foreland. The exposed succession shows a marked sedimentary change from tidal to terrestrial facies. ¹⁴C analyses on gastropod shells indicate the emergence of the Lindi coast at ∼ 44 ¹⁴C ka BP. Emergence and subsequent elevation of terraces to 21 m above present-day sea level was linked to the falling eustatic sea level prior to the last glacial maximum, and to a periodic elevation due to extensional tectonic episodes in the eastern branch of the East African Rift System (EARS). Since ∼ 44 ¹⁴C ka BP tectonic uplift at the coast was 80-110 m, comparable to that in the extreme uplift areas of the EARS.