The spawning and early life-history of a New Zealand endemic amphidromous eleotrid,bluegill bully (Gobiomorphus hubbsi)

The larval stage of fishes is critical in determining their dispersal, survival and recruitment, but little is known of the larval behaviours and tolerances of amphidromous fishes, particularly in New Zealand. We report the results of a series of observational and experimental studies on bluegill bully (Gobiomorphus hubbsi), including spawning sites and behaviours, larval characteristics at hatch, phototactic responses of larvae, and larval survival at different salinity levels. Spawning primarily occurred in the lower reaches of the river, and larvae from different nests exhibited marked differences in, and trade-offs between, larval characteristics at hatch, potentially affecting larval success. Larvae were positively phototactic to intense light, an unexpected result based on diel drift patterns and international research. Finally, larvae exhibited markedly higher survival rates when reared at intermediate salinities compared to freshwater or seawater, suggesting estuaries may play an important role as nursery grounds for bluegill bully and other amphidromous fish.     

A new equation of state for the Aral Sea water

A new equation of state for the Aral Sea is proposed. The equation expresses the dependency of the density on the temperature and salinity. It corresponds to the contemporary salinity range and the ionic salt composition as observed in 2007. The equation was obtained based on direct density measurements using a high precision vibration beam densimeter; therefore, it does not rely on any a priory assumptions about the ionic composition. The accuracy of the density calculations from the new equation of state significantly exceeds that from the extrapolation of the oceanic equation.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

Lower Cambrian patch reefs and associated sediments: southern Labrador, Canada

The complex pattern of biological accretion, internal sedimentation, early lithification, and biological destruction, that characterizes modern reefs and many fossil reefs has been recognized in archaeocyathid-rich patch reefs of Lower Cambrian age in the Forteau Formation, southern Labrador. Patch reefs occur as isolated masses or complex associations of many discrete masses of archaeocyathid-rich limestone and skeletal lime sands, surrounded by well-bedded skeletal limestones and shales. Each reef is composed of many loafshaped mounds stacked on top of one another. The limestone of each mound comprises archaeocyathids and Renalcis or Renalcis-like structures in a matrix of argillaceous lime mud rich in sponge spicules, trilobite and salterellid skeletons. Numerous growth cavities roofed by pendant Renalcis-like organisms and Renalcis are partially to completely filled with geopetal sediment indicating that much of the matrix was deposited as internal sediment. Two stages of diagenetic alteration are recognized: (1) syn-depositional, which affected only the reefs, and (2) post-depositional, which affected both reefs and inter-reef sediments. On the sea floor reef sediments were pervasively cemented and fibrous carbonate was precipitated in intraskeletal and growth cavities. These limestones and cements as well as archaeocyathid skeletons, were subsequently bored by endolithic organisms. Later post-depositional subaerial diagenesis resulted first in dissolution of certain skeletons and precipitation of calcite cement above the water table, followed by extensive precipitation of pore-filling calcite below the water table. These carbonate reefs are similar in structure to the basal pioneer accumulations of much younger lower and middle Palaeozoic reefs. They did not develop into massive ‘ecologic’ reefs because archaeocyathids never developed the necessary large, massive, hemispherical skeletons. This occurrence indicates that reefs developed more or less coincident with, and not long after, the appearance of skeletal metazoans in the Lower Cambrian.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Application of time- and magnitude-predictable model for long-term earthquake prediction in Iran

The regional time- and magnitude-predictable model has been applied successfully in diverse regions of the world to describe the occurrence of main shocks. In the current study, the model has been calibrated against the historical and instrumental catalog of Iranian earthquakes. The Iranian plateau is divided into 15 seismogenic provinces; then, the interevent times for strong main shocks have been determined for each one. The empirical relations reported by Papazachos et al. (Tectonophysics 271:295–323, 1997a) for the Alpine–Himalayan belt (including Iran) were adopted except for the constant terms that were calculated separately for every seismotectonic area. By using the calibrated equations developed for the study area and taking into account the occurrence time and magnitude of the last main shocks in each seismogenic source, the time-dependent conditional probabilities of occurrence P(?t) of the next main shocks during next 10, 20, 30, 40 and 50 years as well as the magnitude of the expected main shocks (M _f) have been estimated. The immediate probability (within next 10 years) of a large main shock is estimated to be high and moderate (>35 %) in all regions except zones 9 (M _f = 5.8) and 15 (M _f = 6.1). However, it should be noted that the probabilities have been estimated for different M _f values in 15 regions. Comparing the model predictions with the actual earthquake occurrence rates shows the good performance of the model for Iranian plateau.     

Electrical and electromagnetic surveys to locate possible causes of water seepage to ground surface at a quarry open pit near Helwan city,Egypt

Water seepage to ground surface at a limestone quarry located at Wadi Garawy about 20 km south-east of Helwan city in Egypt posed a real threat to the mining activity at the quarry. The quarry area is known to be very dry for decades and away from water utilities and infrastructures that may cause water leaks to the quarry. Geophysical investigation including 1D Vertical Electrical Sounding (VES), 2D Electrical Resistivity Tomography (ERT) and 1D Transient Electromagnetic (TEM) surveys were conducted to characterize the rock sequence and locate what could be a possible source of water seepage to the quarry. The resistivity profiles generated from the VES and TEM surveys mapped the rock units in the area down to depths exceeded 100 m. The ERT profiles acquired from the quarrying zone close to the water seepage spot have imaged the top of groundwater level at few meters below the ground surface at the quarry open pit. The spot of groundwater seepage seemed to occur at an area of limestone dissolution that were filled by finer sediments. The finer sediments acted as a hydrological conduit that allowed an upward seepage of groundwater to ground surface under the capillary action effect.     

The marine Lockne impact structure, Jämtland, Sweden: a review

The Lockne impact structure in Jämtland (63°00'20"N, 14°49'30"E) formed in the Middle Ordovician at approximately 455 Ma. The structure is a concentric crater with a total diameter of 13.5 km. The impact took place in a marine environment. Seawater played an important role in the cratering process and in crater morphology and the amount of melt remaining in the structure. Seawater rushed back into the crater in a resurge, eroding and redepositing the ejecta among the resurge deposit. Seawater furthermore facilitated the hydrothermal system, which was driven by the residual heat in the structure. The Lockne structure hosts shocked quartz and an iridium anomaly. The rim wall round the crater collapsed in the modification stage of the crater and was annihilated by the resurge. The fractured basement and the impact breccia were initially rich in open cavities. These became partly filled with dominantly calcite. The filling contributed to a low-density contrast, generating a negative gravity anomaly of 22 gu. The gravity model indicates a central uplift and a NW-directed tilt of the structure. This tilt is also seen in the magnetic models. The apparent absence of any impact melt is probably real and related to the environment of impact.     

Search for variable sources using data of “Cold” surveys

We search for variable sources, using the data of the surveys conducted on the RATAN-600 radio telescope in 1980–1994 at 3.94 GHz. To test the radio sources of the RCR (RATAN Cold Refined) catalog for variability, we estimated the long-term variability indices V of the studied objects, their relative variability amplitudes V _χ , and the χ ² probabilities p. Out of about two hundred considered sources, 41 proved to have positive long-term variability indices, suggesting that these sources may be variable. Fifteen objects can be considered to be reliably variable according to the χ ² criterion p > 0.98, three of these sources have χ ² probabilities p ≥ 0.999. The corresponding probabilities for six sources lie in the 0.95 < p < 0.98 interval, and those of the remaining 20 objects in the 0.73 ≤ p < 0.95 interval. Twenty four of 41 objects are variable or possibly variable in the optical range, and five objects are known variable radio sources. We construct the light curves and spectra for the sources with positive long-term variability indices.     

Integer least-squares theory for the GNSS compass

Global navigation satellite system (GNSS) carrier phase integer ambiguity resolution is the key to high-precision positioning and attitude determination. In this contribution, we develop new integer least-squares (ILS) theory for the GNSS compass model, together with efficient integer search strategies. It extends current unconstrained ILS theory to the nonlinearly constrained case, an extension that is particularly suited for precise attitude determination. As opposed to current practice, our method does proper justice to the a priori given information. The nonlinear baseline constraint is fully integrated into the ambiguity objective function, thereby receiving a proper weighting in its minimization and providing guidance for the integer search. Different search strategies are developed to compute exact and approximate solutions of the nonlinear constrained ILS problem. Their applicability depends on the strength of the GNSS model and on the length of the baseline. Two of the presented search strategies, a global and a local one, are based on the use of an ellipsoidal search space. This has the advantage that standard methods can be applied. The global ellipsoidal search strategy is applicable to GNSS models of sufficient strength, while the local ellipsoidal search strategy is applicable to models for which the baseline lengths are not too small. We also develop search strategies for the most challenging case, namely when the curvature of the non-ellipsoidal ambiguity search space needs to be taken into account. Two such strategies are presented, an approximate one and a rigorous, somewhat more complex, one. The approximate one is applicable when the fixed baseline variance matrix is close to diagonal. Both methods make use of a search and shrink strategy. The rigorous solution is efficiently obtained by means of a search and shrink strategy that uses non-quadratic, but easy-to-evaluate, bounding functions of the ambiguity objective function. The theory presented is generally valid and it is not restricted to any particular GNSS or combination of GNSSs. Its general applicability also applies to the measurement scenarios (e.g. single-epoch vs. multi-epoch, or single-frequency vs. multi-frequency). In particular it is applicable to the most challenging case of unaided, single frequency, single epoch GNSS attitude determination. The success rate performance of the different methods is also illustrated.