A computational approach to unbiased districting

In the context of discrete districting problems with geographical constraints, we demonstrate that determining an (ex post) unbiased districting, which requires that the number of representatives of a party should be proportional to its share of votes, turns out to be a computationally intractable (NP-complete) problem. This raises doubts as to whether an independent jury will be able to come up with a “fair” redistricting plan in case of a large population, that is, there is no guarantee for finding an unbiased districting (even if such exists). We also show that, in the absence of geographical constraints, an unbiased districting can be implemented by a simple alternating-move game among the two parties.

An LP-based inconsistency monitoring of pairwise comparison matrices

A distance-based inconsistency indicator, defined by the third author for the consistency-driven pairwise comparisons method, is extended to the incomplete case. The corresponding optimization problem is transformed into an equivalent linear programming problem. The results can be applied in the process of filling in the matrix as the decision maker gets automatic feedback. As soon as a serious error occurs among the matrix elements, even due to a misprint, a significant increase in the inconsistency index is reported. The high inconsistency may be alarmed not only at the end of the process of filling in the matrix but also during the completion process. Numerical examples are also provided.

On optimal completions of incomplete pairwise comparison matrices

An important variant of a key problem for multi-attribute decision making is considered. We study the extension of the pairwise comparison matrix to the case when only partial information is available: for some pairs no comparison is given. It is natural to define the inconsistency of a partially filled matrix as the inconsistency of its best, completely filled completion. We study here the uniqueness problem of the best completion for two weighting methods, the Eigen-vector Method and the Logarithmic Least Squares Method. In both settings we obtain the same simple graph theoretic characterization of the uniqueness. The optimal completion will be unique if and only if the graph associated with the partially defined matrix is connected. Some numerical experiences are discussed at the end of the paper.

Numerical Investigation of the Effect of Magnetic Field on Natural Convection in a Curved-Shape Enclosure

This investigation reports the magnetic field effect on natural convection heat transfer in a curved-shape enclosure. The numerical investigation is carried out using the control volume-based-finite element method (CVFEM). The numerical investigations are performed for various values of Hartmann number and Rayleigh number. The obtained results are depicted in terms of streamlines and isotherms which show the significant effects of Hartmann number on the fluid flow and temperature distribution inside the enclosure. Also, it was found that the Nusselt number decreases with an increase in the Hartmann number.

Downhole logging of the AND-2A borehole, Victoria Land Basin, McMurdo Sound, Antarctica

Under the framework of the ANDRILL Southern McMurdo Sound (SMS) Project successful downhole experiments were conducted in the 1138.54 metre (m)-deep AND-2A borehole. Wireline logs successfully recorded were: magnetic susceptibility, spectral gamma ray, sonic velocity, borehole televiewer, neutron porosity, density, calliper, geochemistry, temperature and dipmeter. A resistivity tool and its backup both failed to operate, thus resistivity data were not collected. Due to hole conditions, logs were collected in several passes from the total depth at ~1138 metres below sea floor (mbsf) to ~230 mbsf, except for some intervals that were either inaccessible due to bridging or were shielded by the drill string. Furthermore, a Vertical Seismic Profile (VSP) was created from ~1000 mbsf up to the sea floor. The first hydraulic fracturing stress measurements in Antarctica were conducted in the interval 1000-1138 mbsf. This extensive data set will allow the SMS Science Team to reach some of the ambitious objectives of the SMS Project. Valuable contributions can be expected for the following topics: cyclicity and climate change, heat flux and fluid flow, seismic stratigraphy in the Victoria Land Basin, and structure and state of the modern crustal stress field.

Fault-controlled and stratabound dolostones in the Late Aptian-earliest Albian Benassal Formation (Maestrat Basin, E Spain) : petrology and geochemistry constrains

This study was developed under the ExxonMobil FC2 Alliance (Fundamental Controls on Flow in Carbonates). The authors wish to thank ExxonMobil Production Company and ExxonMobil Upstream Research Company for providing funding. The views in this article by Sherry L. Stafford are her own and not necessarily those of ExxonMobil. This research was supported by the Sedimentary Geology Research Group of the Generalitat de Catalunya (2014SGR251). We would like to thank Andrea Ceriani and Paola Ronchi for their critical and valuable reviews, and Associated Editor Piero Gianolla for the editorial work.

Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Biogeochemical investigation of sediments at the Hakon Mosby Mud Volcano (Barents Sea)

We have investigated if in a cold seep methane or sulfide is used for chemosynthetic primary production and if significant amounts of the sulfide produced by anaerobic oxidation of methane are oxidized geochemically and hence are not available for chemosynthetic production. Geochemically controlled redox reactions and biological turnover were compared in different habitats of the Håkon Mosby Mud Volcano. The center of the mud volcano is characterized by the highest fluid flow, and most primary production by the microbial community depends on oxidation of methane. The small amount of sulfide produced is oxidized geochemically with oxygen or is precipitated with dissolved iron. In the medium flow peripheral Beggiatoa habitat sulfide is largely oxidized biologically. The oxygen and nitrate supply is high enough that Beggiatoa can oxidize the sulfide completely, and chemical sulfide oxidation or precipitation is not important. An internally stored nitrate reservoir with average concentrations of 110 mmol L-1 enables the Beggiatoa to oxidize sulfide anaerobically. The pH profile indicates sequential sulfide oxidation with elemental sulfur as intermediate. Gray thiotrophic mats associated with perturbed sediments showed a high heterogeneity in sulfate turnover and high sulfide fluxes, balanced by the opposing oxygen and nitrate fluxes so that biological oxidation dominates over geochemical sulfide removal processes. The three habitats indicate substantial small-scale variability in carbon fixation pathways either through direct biological use of methane or through indirect carbon fixation of methane-derived carbon dioxide by chemolithotrophic sulfide oxidation.

Strontium and neodymium concentrations and isotopic compositions of Arctic marine sediments

Sr and Nd isotopic compositions of Arctic marine sediments characterize changes of sediment source regions and trace shelf-ocean particle pathways during glacial-interglacial transitions in the eastern Arctic Ocean. In the 140-ka sedimentary record of a marine core from Yermak Plateau, north of Svalbard, 87Sr/86Sr ratios and epsion-Nd values vary between 0.717 and 0.740 and 39.3 and 314.9, respectively. Sr and Nd isotopic composition both change characteristically during glacial-interglacial cycles and are correlated with the extension of the Svalbard/Barents Sea ice sheet (SBIS). The downcore variation in Sr and Nd isotopic composition indicates climatically induced changes in sediment provenance from two isotopically distinct end-members: (1) Eurasian shelf sediments as a distal source; and (2) Svalbard bedrock as a proximal source that coincide with a change in transport mechanism from sea ice to glacial ice. During glacier advance from Svalbard and intensified glacial bedrock erosion, epsion-Nd values decrease gradually to a minimum value of 314.9 due to increased input of crystalline Svalbard bedrock material. During glacial maxima, the SBIS covered the entire Barents Sea shelf and supplied increasing amounts of Eurasian shelf material to the Arctic Ocean as ice rafted detritus (IRD). Epsion-Nd values in glacial sediments reach maximum values that are comparable to the average value of modern Eurasian shelf and sea ice sediments (epsion-Nd = 310.3). This confirms ice rafting as a major sediment transport mechanism for Eurasian shelf sediments into the Arctic Ocean and trace a sediment origin from the Kara Sea/Laptev Sea shelf area. After the decay of the shelf-based SBIS, the glacial shelf sediment spikes during glacial terminations I (epsion-Nd = 310.6) and II (epsion-Nd = 310.1) epsion-Nd values rapidly decrease to values of 312.5 typical for interglacial averages. The downcore Sr isotopic composition is anticorrelated to the Nd isotopic composition, but may be also influenced by grain-size effects. In contrast, the Nd isotopic composition in clay- to silt-size fractions of one bulk sediment sample is similar to within 0.3-0.8 epsion-Nd units and seems to be a grain-size independent provenance tracer.

Beryllium 7 measured during POLARSTERN cruise ARK-XXII/2

The natural cosmogenic radionuclide 7Be (T1/2 = 53.4 d) is supplied to the surface ocean from the atmosphere and, in the Arctic Ocean, can be used as a tracer of the efficiency with which sea ice intercepts the atmospheric fluxes of chemical species and of the importance of ice as a transport mechanism for particulate matter and chemical species. Analyses of 7Be in samples of surface water, surface sea ice, water beneath the ice, sea ice sediments, and precipitation from the Eurasian Basin of the Arctic Ocean show that the fraction of sea ice coverage determines the amount of 7Be in the surface water. When sea ice coverage is <40%, the 7Be inventory in the upper ocean (130 ± 19 Bq m**-2) is in good agreement with that expected from the inventory from 7Be atmospheric flux (128 ± 21 Bq m**-2). In contrast, when ice coverage is >80%, the water column inventory drops to 58 ± 20 Bq m**-2. The 7Be inventory in sea ice is 39 ± 23 Bq m**-2, and mass balance calculations show that sea ice can intercept 30 ± 18% of the atmospheric flux of 7Be during the studied period. We suggest that other atmospherically transported contaminants should be similarly intercepted. 7Be in the ice also can be used to estimate that the annual transport and release of sediment to the ablation area of the Fram Strait is -500 g m**-2, a value comparable to previously measured fluxes in sediment traps deployed in the area.

Seawater carbonate chemistry and impacts of calcium and carbonate ion concentration on Mg and Sr incorporation in cultured benthic foraminifera (Ammonia tepida and Heterostegina depressa), 2011

Laboratory culture experiments were conducted to determine effects of seawater carbonate ion concentration ([CO32-]), and thereby calcite saturation state, on Mg and Sr incorporation into calcite of two species of shallow-water benthic foraminifera: Ammonia tepida and Heterostegina depressa. Impact on Mg and Sr incorporation by increased seawater [CO32-] and thereby higher calcite saturation state, is absent in either species. Comparison to results from a similar culturing experiment, in which calcite saturation state was varied as a function of [Ca2+], reveals that saturation state affects incorporation of Mg and Sr through calcium- rather than carbonate availability. The similarity in response by both species is surprising since the average Mg/Ca ratio is ~ 70 times higher in H. depressa than in A. tepida. Furthermore, these results suggest that the ions involved in biomineralization (i.e. Ca2+ and DIC) are processed by separate cellular transport mechanisms. The similar response of Mg and Sr incorporation in this study suggests that only differences in the Ca2+ transport mechanism affect divalent cation partitioning.

Terrigenous input and marine productivity for the Pacific sector of the Southern Ocean based on lipid biomarkers

In this study, we present a new multiproxy data set of terrigenous input, marine productivity and sea surface temperature (SST) from 52 surface sediment samples collected along E-W transects in the Pacific sector of the Southern Ocean. Allochtonous terrigenous input was characterized by the distribution of plant wax n-alkanes and soil-derived branched glycerol dialkyl glycerol tetraethers (brGDGTs). 230Th-normalized burial rates of both compound groups were highest close to the potential sources in Australia and New Zealand and are strongly related to lithogenic contents, indicating common sources and transport. Detection of both long-chain n-alkanes and brGDGTs at the most remote sites in the open ocean strongly suggests a primarily eolian transport mechanism to at least 110°W, i.e. by prevailing westerly winds. Two independent organic SST proxies were used, the UK'37 based on long-chain alkenones, and the TEX86 based on isoprenoid GDGTs. Both, UK'37 and TEX86 indices show robust relationships with temperature over a temperature range between 0.5 and 20°C, likely implying different seasonal and regional imprints on the temperature signal. While alkenone-based temperature estimates reliably reflect modern SST even at the low temperature end, large temperature residuals are observed for the polar ocean using the TEX86 index. 230Th-normalized burial rates of alkenones are highest close to the Subtropical Front and are positively related to lithogenic fluxes throughout the study area. In contrast, highest isoGDGT burial south of the Antarctic Polar Front is not related with dust flux but may be largely controlled by diatom blooms, and thus high opal fluxes during austral summer.

Temperature measurements, gas chemistry, and gas hydrates of the Dvurechenskii mud volcano

The sediment temperature distribution at mud volcanoes provides insights into their activity and into the occurrence of gas hydrates. If ambient pressure and temperature conditions are close to the limits of the gas hydrate stability field, the sediment temperature distribution not only limits the occurrence of gas hydrates, but is itself influenced by heat production and consumption related to the formation and dissociation of gas hydrates. Located in the Sorokin Trough in the northern Black Sea, the Dvurechenskii mud volcano (DMV) was in the focus of detailed investigations during the M72/2 and M73/3a cruises of the German R/V Meteor and the ROV Quest 4000 m in February and March 2007. A large number of in-situ sediment temperature measurements were conducted from the ROV and with a sensor-equipped gravity corer. Gas hydrates were sampled in pressurized cores using a dynamic autoclave piston corer (DAPC). The thermal structure of the DMV suggests a regime of fluid flow at rates decreasing from the summit towards the edges of the mud volcano, accompanied by intermittent mud expulsion at the summit. Modeled gas hydrate dissociation temperatures reveal that the gas hydrates at the DMV are very close to the stability limits. Changes in heat flow due to variable seepage rates probably do not result in changes in sediment temperature but are compensated by gas hydrate dissociation and formation.

On the feasibility of time-lapse superconducting gravimetry for reservoir monitoring

The feasibility of monitoring fluid flow subsurface processes that result in density changes, using the iGrav superconducting gravimeter, is investigated. Practical targets include steam-assisted gravity drainage (SAGD) bitumen depletion and water pumping from aquifers, for which there is currently a void in low-impact, inexpensive monitoring techniques. This study demonstrates that the iGrav has the potential to be applied to multi-scale and diverse reservoirs. Gravity and gravity gradient signals are forward modeled for a real SAGD reservoir at two time steps, and for surface-fed and groundwater-fed aquifer pumping models, to estimate signal strength and directional dependency of water flow. Time-lapse gravimetry on small-scale reservoirs exhibits two obstacles, namely, a µgal sensitivity requirement and high noise levels in the vicinity of the reservoir. In this study, both limitations are overcome by proposing (i) a portable superconducting gravimeter, and (ii) a pair of instruments under various baseline geometries. This results in improved spatial resolution for locating depletion zones, as well as the cancellation of noise common in both instruments. Results indicate that a pair of iGrav superconducting gravimeters meet the sensitivity requirements and the spatial focusing desired to monitor SAGD bitumen migration at the reservoir scales. For SAGD reservoirs, the well pair separation, reservoir depth, and survey sampling determine the resolvability of individual well pair depletion patterns during the steam chamber rising phase, and general reservoir depletion patterns during the steam chamber spreading phase. Results show that monitoring water table elevation changes due to pumping and tracking whether groundwater or surface water is being extracted are feasible.

Parametric Analysis of Concrete Solar Collectors

Concrete solar collectors offer a type of solar collector with structural, aesthetic and economic advantages over current populartechnologies. This study examines the influential parameters of concrete solar collectors. In addition to the external conditions,the performance of a concrete solar collector is influenced by the thermal properties of the concrete matrix, piping network andfluid. Geometric and fluid flow parameters also influence the performance of the concrete solar collector. A literature review ofconcrete solar collectors is conducted in order to define the benchmark parameters from which individual parameters are thencompared. The numerical model consists of a 1D pipe flow network coupled with the heat transfer in a 3D concrete domain. Thispaper is concerned with the physical parameters that define the concrete solar collector, thus a constant surface temperature isused as the exposed surface boundary condition with all other surfaces being insulated. Results show that, of the parametersinvestigated, the pipe spacing, ps, concrete conductivity, kc, and the pipe embedment depth, demb, are among those parameterswhich have greatest effect on the collector’s performance. The optimum balance between these parameters is presented withrespect to the thermal performance and discussed with reference to practical development issues.