Modelling basin-scale distribution of fish occurrence probability for assessment of flow and habitat conditions in rivers

Abstract

Flow regimes play an important role in sustaining biodiversity in river ecosystems. However, the effects of flow regimes on riverine fish have not been clearly described. Therefore, we propose a new methodology to quantitatively link habitat conditions (such as flow indices and physical habitat conditions) to the occurrence probability (OP) of fish species. We developed a basin-scale fish distribution model by integrating the concept of habitat suitability assessment with a distributed hydrological model in order to estimate the OP of fish, with particular attention to flow regime. A generalized linear model was used to evaluate the relationship between the probabilities of fish occurrence and major environmental factors in river sections. A geomorphology-based hydrological model was adopted to simulate river discharge, which was used to calculate 10 flow indices. The occurrence probabilities of 50 fish species in the Sagami River in Japan were modelled. For the prediction accuracy, field survey results that included at least five observations of both the presence and the absence of each species were required to obtain relatively reliable prediction (accuracy > 60%). Using the developed model, important habitat conditions for each species were identified, which showed the importance of low-flow events for more than 10 species, including Hypomesus nipponensis and Rhinogobius fluviatilis. The model also confirmed the positive effects of natural flow and the negative effect of river-crossing structures, such as dams and weirs, on the OP of most species. The suggested approach enables us to evaluate and project the ecological consequences of water resource management policy. The results demonstrate the applicability of the fish distribution model to provide quantitative information on the flow required to maintain fish communities.

Editor Z.W. Kundzewicz; Guest editor M. Acreman

Citation Sui, P., Iwasaki, A., Saavedra, V.O.C., and Yoshimura, C., 2013. Modelling basin-scale distribution of fish occurrence probability for assessment of flow and habitat conditions in rivers. Hydrological Sciences Journal, 59 (3–4), 618–628.     

Sources and geochemical constraints of polycyclic aromatic hydrocarbons (PAHs) in sediments and mussels of two Northern Irish Sea-loughs.

Polycyclic aromatic hydrocarbons (PAHs) were sampled from 30 sediment and 8 mussel (Mytilus edulis) stations in two Northern Irish Sea-toughs. Analysis was performed by gas chromatography coupled to mass spectrometry (GC-MS). Sedimentary organic carbon, % silt/clay and mean grain size were analysed in order to assess the role of geochemistry on PAH distribution. With the exception of two sites in Larne Lough representing localized regions of high contamination, sum(PAHs) in sediments ranged between 83 and 2300 ng g(-1). Regression analysis indicated that particle size and organic C were dominant factors in controlling the distribution of PAHs throughout the sediments. Sources of PAHs in both sea-loughs were dominated by pyrogenic inputs suggesting that diffuse sources such as atmospheric deposition may be a major source to both water bodies. The sum of PAHs in mussels ranged between 95 and 184 ng g(-1). Variations in concentrations may be explained by differences in the condition between mussel populations. Mussels in both sea-loughs exhibit similar metabolic activities towards the PAH compounds which were predominantly bioaccumulated from sediments.     

An Early Pleistocene interglacial deposit at Pingorsuit,North-West Greenland

At the Pingorsuit Glacier in North-West Greenland, an organic-rich deposit that had recently emerged from the retreating ice cap was discovered at an elevation of 480 m above sea level. This paper reports on macrofossil analyses of a coarse detritus gyttja and peaty soil, which occurred beneath a thin cover of till and glacifluvial deposits. The sediments contained remains of vascular plants, mosses, beetles, caddisflies, midges, bryozoans, sponges and other invertebrates. The flora includes black spruce, tree birch, boreal shrubs and wetland and aquatic taxa, which shows that mires, lakes and ponds were present in the area. We describe a new extinct waterwort species Elatine odgaardii. The fossils were deposited in a boreal environment with a mean July air temperature that was at least 9 °C higher than at present. The fossil assemblages show strong similarities with others from Greenland that have been assigned an Early Pleistocene age, and we suggest a similar age for the sediments found at the margin of the Pingorsuit Glacier. At the Pingorsuit Glacier in North-West Greenland, an organic-rich deposit was discovered at an elevation of 480 m above sea level. The sediments contained remains of vascular plants, mosses, beetles, caddisflies, midges, bryozoans, sponges and other invertebrates. The fossils were deposited in a boreal environment with a mean July air temperature that was at least 9 °C higher than at present.     

Efficient reservoir history matching using subspace vectors

In this paper, we describe a method of history matching in which changes to the reservoir model are constructed from a limited set of basis vectors. The purpose of this reparameterization is to reduce the cost of a Newton iteration, without altering the final estimate of model parameters and without substantially slowing the rate of convergence. The utility of a subspace method depends on several factors, including the choice and number of the subspace vectors to be used. Computational gains in efficiency result partly from a reduction in the size of the matrix system that must be solved in a Newton iteration. More important contributions, however, result from a reduction in the number of sensitivity coefficients that must be computed, reduction in the dimensions of the matrices that must be multiplied, and elimination of matrix products involving the inverse of the prior model covariance matrix. These factors affect the efficiency of each Newton iteration. Although computation of the optimal set of subspace vectors may be expensive, we show that the rate of convergence and the final results are somewhat insensitive to the choice of subspace vectors. We also show that it is desirable to start with a small number of subspace vectors and gradually increase the number at each Newton iteration until an acceptable level of data mismatch is obtained.     

Tracking and orbit determination performance of the GRAS instrument on MetOp-A

The global navigation satellite system receiver for atmospheric sounding (GRAS) on MetOp-A is the first European GPS receiver providing dual-frequency navigation and occultation measurements from a spaceborne platform on a routine basis. The receiver is based on ESA’s AGGA-2 correlator chip, which implements a high-quality tracking scheme for semi-codeless P(Y) code tracking on the L1 and L2 frequency. Data collected with the zenith antenna on MetOp-A have been used to perform an in-flight characterization of the GRAS instrument with focus on the tracking and navigation performance. Besides an assessment of the receiver noise and systematic measurement errors, the study addresses the precise orbit determination accuracy achievable with the GRAS receiver. A consistency on the 5 cm level is demonstrated for reduced dynamics orbit solutions computed independently by four different agencies and software packages. With purely kinematic solutions, 10 cm accuracy is obtained. As a part of the analysis, an empirical antenna offset correction and preliminary phase center correction map are derived, which notably reduce the carrier phase residuals and improve the consistency of kinematic orbit determination results.

Antenna phase center calibration for precise positioning of LEO satellites

Phase center variations of the receiver and transmitter antenna constitute a remaining uncertainty in the high precision orbit determination (POD) of low Earth orbit (LEO) satellites using GPS measurements. Triggered by the adoption of absolute phase patterns in the IGS processing standards, a calibration of the Sensor Systems S67-1575-14 antenna with GFZ choke ring has been conducted that serves as POD antenna on various geodetic satellites such as CHAMP, GRACE and TerraSAR-X. Nominal phase patterns have been obtained with a robotic measurement system in a field campaign and the results were used to assess the impact of receiver antenna phase patterns on the achievable positioning accuracy. Along with this, phase center distortions in the actual spacecraft environment were characterized based on POD carrier phase residuals for the GRACE and TerraSAR-X missions. It is shown that the combined ground and in-flight calibration can improve the carrier phase modeling accuracy to a level of 4 mm which is close to the pure receiver noise. A 3.5 cm (3D rms) consistency of kinematic and reduced dynamic orbit determination solutions is achieved for TerraSAR-X, which presumably reflects the limitations of presently available GPS ephemeris products. The reduced dynamic solutions themselves match the observations of high grade satellite laser ranging stations to 1.5 cm but are potentially affected by cross-track biases at the cm-level. With respect to the GPS based relative navigation of TerraSAR-X/TanDEM-X formation, the in-flight calibration of the antenna phase patterns is considered essential for an accurate modeling of differential carrier phase measurements and a mm level baseline reconstruction.

Precision spacecraft navigation using a low-cost GPS receiver

Within the PROBA-2 microsatellite mission, a miniaturized single-frequency GPS receiver based on commercial-off-the-shelf (COTS) technology is employed for onboard navigation and timing. A rapid electronic fuse protects against destructive single-event latch-ups (SEL) and enables a quasi-continuous receiver operation despite the inherent sensitivity to space radiation. While limited to single-frequency C/A-code tracking with a narrow-band frontend, the receiver is able to provide precision navigation services through processing of raw GPS measurements on ground as well as a built-in real-time navigation system. In both cases, ionospheric path delays are eliminated through a combination of L1 pseudorange and carrier phase measurements, which also offers a factor-of-two noise reduction relative to code-only processing. By comparison with satellite laser ranging (SLR) measurements, a 0.3-m (3D rms) accuracy is demonstrated for the PROBA-2 reduced dynamic orbit determinations using post-processed GPS orbit and clock products. Furthermore, the experimental onboard navigation system is shown to provide real-time position information with a 3D rms accuracy of about 1?m, which notably outperforms the specification of the Standard Positioning Service (SPS). In view of their lower hardware complexity, mass budget and power requirements as well as the reduced interference susceptibility, legacy C/A-code receivers can thus provide an attractive alternative to dual-frequency receivers even for demanding navigation applications in low Earth orbit.     

Short-term analysis of GNSS clocks

A characterization of the short-term stability of the atomic frequency standards onboard GNSS satellites is presented. Clock performance is evaluated using two different methods. The first method derives the temporal variation of the satellite’s clock from a polynomial fit through 1-way carrier-phase measurements from a receiver directly connected to a high-precision atomic frequency standard. Alternatively, three-way measurements using inter-station single differences of a second satellite from a neighboring station are used if the receiver’s clock stability at the station tracking the satellite of interest is not sufficient. The second method is a Kalman-filter-based clock estimation based on dual-frequency pseudorange and carrier-phase measurements from a small global or regional tracking network. Both methods are introduced and their respective advantages and disadvantages are discussed. The analysis section presents a characterization of GPS, GLONASS, GIOVE, Galileo IOV, QZSS, and COMPASS clocks based on these two methods. Special focus has been set on the frequency standards of new generation satellites like GPS Block IIF, QZSS, and IOV as well as the Chinese COMPASS/BeiDou-2 system. The analysis shows results for the Allan deviation covering averaging intervals from 1 to 1,000 s, which is of special interest for real-time PPP and other high-rate applications like processing of radio-occultation measurements. The clock interpolation errors for different sampling rates are evaluated for different types of clocks and their effect on PPP is discussed.     

Kalman-filter-based GPS clock estimation for near real-time positioning

In this article, an algorithm for clock offset estimation of the GPS satellites is presented. The algorithm is based on a Kalman-filter and processes undifferenced code and carrier-phase measurements of a global tracking network. The clock offset and drift of the satellite clocks are estimated along with tracking station clock offsets, tropospheric zenith path delay and carrier-phase ambiguities. The article provides a brief overview of already existing near-real-time and real-time clock products. The filter algorithm and data processing scheme is presented. Finally, the accuracy of the orbit and clock product is assessed with a precise orbit determination of the MetOp satellite and compared to results gained with other real-time products.