Paleoclimatic record of a long deep sea core from the eastern Mediterranean

A deep-sea core over 16 m long from the crestal area of the Mediterranean Ridge has been investigated with different techniques, including quantitative micropaleontology, stable isotopes (measured on the epipelagic species Globigerinoides ruber and on the mesopelagic species Globorotalia inflata), and clay mineralogy. The resulting record of climatic fluctuations can be cross correlated to other Mediterranean cores by means of isochronous lithologies (tephra layers and sapropels). The climatic record of the Mediterranean is similar in character, phase, and chronology to the records investigated in the equatorial Pacific and in the Caribbean. Isotope stages 1 to 17 have been recognized. Cyclically repeated stagnant cycles resulting in sapropel deposition complicate both the isotopic and the faunal signal. The isotopic investigations reveal that the temperature change in the surface layers of the eastern Mediterranean was no greater than 8°C in the late “glacial” Pleistocene. The chronostratigraphic and biostratigraphic interpretation of Core KS09 indicate that the mean sedimentation rate was 2.4 cm/1000 years, a value very close to the 2.5 cm/1000 years calculated for the entire Quaternary section at DSDP Site 125, also located in the crestal area of the Mediterranean Ridge in the Ionian Basin. The base of KS09 is likely to be very close to the Brunhes/Matuyama boundary dated at 0.7 my.     

Texture and anisotropy analysis of Qusaiba shales

Scanning and transmission electron microscopy, synchrotron X‐ray diffraction, microtomography and ultrasonic velocity measurements were used to characterize microstructures and anisotropy of three deeply buried Qusaiba shales from the Rub’al‐Khali basin, Saudi Arabia. Kaolinite, illite‐smectite, illite‐mica and chlorite show strong preferred orientation with (001) pole figure maxima perpendicular to the bedding plane ranging from 2.4–6.8 multiples of a random distribution (m.r.d.). Quartz, feldspars and pyrite crystals have a random orientation distribution. Elastic properties of the polyphase aggregate are calculated by averaging the single crystal elastic properties over the orientation distribution, assuming a nonporous material. The average calculated bulk P‐wave velocities are 6.2 km/s (maximum) and 5.5 km/s (minimum), resulting in a P‐wave anisotropy of 12%. The calculated velocities are compared with those determined from ultrasonic velocity measurements on a similar sample. In the ultrasonic experiment, which measures the effects of the shale matrix as well as the effects of porosity, velocities are smaller (P‐wave maximum 5.3 km/s and minimum 4.1 km/s). The difference between calculated and measured velocities is attributed to the effects of anisotropic pore structure and to microfractures present in the sample, which have not been taken into account in the matrix averaging.     

Ecotoxicological effects of a semi-submerged municipal dump (Castle harbour, Bermuda) on the Calico scallop Argopecten gibbus

A biomarker study was undertaken using the Calico scallop Argopecten gibbus to assess the ecotoxicological effects of a semi-submerged municipal dump on the adjacent patch reef lagoon ecosystem (Castle harbour, Bermuda). Caged scallops were deployed in situ for 2 months at various distances from the dump (50 m, 900 m and 2.7 km) and at a reference site (14 km). A suite of biomarkers comprising metallothionein (MT), lipid peroxidation (LPO), vitellin-like proteins (Vn), glutathione S-transferase (GST), DNA strand breaks and condition factor (CF) were investigated in various tissues of the scallop (gill, digestive gland and gonad). Levels of heavy metals were also measured in the whole scallop soft tissue. While there was some variation in response between tissues, in general the results indicated that the dump was negatively impacting scallops deployed in the adjacent marine environment: elevated levels of MT, DNA strand breaks, Vn and GST and reduced condition factor were found for scallops deployed nearest to the dump and at the site 1.5 km from this point source of contamination (Tuckers town) in Castle harbour, with respect to the reference site, North Rock (although this exhibited some degree of metal contamination). The gills from scallops deployed at the dump site were the most responsive tissue, with the highest expression of MT, LPO and DNA damage. This study indicates the potential of the Calico scallop as a convenient bioindicator species in the marine tropical benthos.     

Ultrasonic velocities of North Sea chalk samples: influence of porosity, fluid content and texture

We have studied 56 unfractured chalk samples of the Upper Cretaceous Tor Formation of the Dan, South Arne and Gorm Fields, Danish North Sea. The samples have porosities of between 14% and 45% and calcite content of over 95%. The ultrasonic compressional‐ and shear‐wave velocities (V_P and V_S) for dry and water‐saturated samples were measured at up to 75 bar confining hydrostatic pressure corresponding to effective stress in the reservoir. The porosity is the main control of the ultrasonic velocities and therefore of the elastic moduli. The elastic moduli are slightly higher for samples from the South Arne Field than from the Dan Field for identical porosities. This difference may be due to textural differences between the chalk at the two locations because we observe that large grains (i.e. filled microfossils and fossil fragments) that occur more frequently in samples from the Dan Field have a porosity‐reducing effect and that samples rich in large grains have a relatively low porosity for a given P‐wave modulus. The clay content in the samples is low and is mainly represented by either kaolinite or smectite; samples with smectite have a lower P‐wave modulus than samples with kaolinite at equal porosity. We find that ultrasonic V_P and V_S of dry chalk samples can be satisfactorily estimated with Gassmann's relationships from data for water‐saturated samples. A pronounced difference between the V_P/V_S ratios for dry and water‐saturated chalk samples indicates promising results for seismic amplitude‐versus‐offset analyses.     

Impact of Local Groundwater Flow Model Errors on Transport and a Practical Solution for the Issue

A groundwater flow model is typically used to provide the flow field for conducting groundwater solute transport simulations. The advection term of the mass conserved formulation for groundwater transport assumes that the flow field is perfectly balanced and that all water flowing into a numerical grid cell is exactly balanced by outflows after accounting for sources/sinks or internal storage. However, in many complicated regional or site‐scale models, there may be localized flow balance errors that may be difficult to eliminate through tighter flow convergence tolerances due to simulation time constraints or numerical limits on convergence tolerances. Thus, if water is erroneously gained or lost within a grid cell during the flow computation, the solutes within it will also be numerically affected in the associated transport simulation. Transport solutions neglect this error in groundwater flow as the transport equations that are solved assume no error in flow. This flow imbalance error can however have consequences on the transport solution ranging from unnoticeable errors in the resulting concentrations to spurious oscillations that can grow in time and hinder further solution. An approach has been suggested here, to explicitly handle these flow imbalances during mass conserved advective transport computations and report them in the corresponding transport mass balance output, as corrections that are needed to handle errors originating in the flow solution. Example problems are provided to explain the concepts and demonstrate the impacts.     

Definition of Physical Height Systems for Telluric Planets and Moons

In planetary sciences, the geodetic (geometric) heights defined with respect to the reference surface (the sphere or the ellipsoid) or with respect to the center of the planet/moon are typically used for mapping topographic surface, compilation of global topographic models, detailed mapping of potential landing sites, and other space science and engineering purposes. Nevertheless, certain applications, such as studies of gravity-driven mass movements, require the physical heights to be defined with respect to the equipotential surface. Taking the analogy with terrestrial height systems, the realization of height systems for telluric planets and moons could be done by means of defining the orthometric and geoidal heights. In this case, however, the definition of the orthometric heights in principle differs. Whereas the terrestrial geoid is described as an equipotential surface that best approximates the mean sea level, such a definition for planets/moons is irrelevant in the absence of (liquid) global oceans. A more natural choice for planets and moons is to adopt the geoidal equipotential surface that closely approximates the geometric reference surface (the sphere or the ellipsoid). In this study, we address these aspects by proposing a more accurate approach for defining the orthometric heights for telluric planets and moons from available topographic and gravity models, while adopting the average crustal density in the absence of reliable crustal density models. In particular, we discuss a proper treatment of topographic masses in the context of gravimetric geoid determination. In numerical studies, we investigate differences between the geodetic and orthometric heights, represented by the geoidal heights, on Mercury, Venus, Mars, and Moon. Our results reveal that these differences are significant. The geoidal heights on Mercury vary from ? 132 to 166 m. On Venus, the geoidal heights are between ? 51 and 137 m with maxima on this planet at Atla Regio and Beta Regio. The largest geoid undulations between ? 747 and 1685 m were found on Mars, with the extreme positive geoidal heights under Olympus Mons in Tharsis region. Large variations in the geoidal geometry are also confirmed on the Moon, with the geoidal heights ranging from ? 298 to 461 m. For comparison, the terrestrial geoid undulations are mostly within ± 100 m. We also demonstrate that a commonly used method for computing the geoidal heights that disregards the differences between the gravity field outside and inside topographic masses yields relatively large errors. According to our estimates, these errors are ? 0.3/+ 3.4 m for Mercury, 0.0/+ 13.3 m for Venus, ? 1.4/+ 125.6 m for Mars, and ? 5.6/+ 45.2 m for the Moon.     

Experiences with the use of mass-density maps in residual gravity forward modelling

In many modern local and regional gravity field modelling concepts, the short-wavelength gravitational signal modeled by the residual terrain modelling (RTM) technique is used to augment global geopotential models, or to smooth observed gravity prior to data gridding. In practice, the evaluation of RTM effects mostly relies on a constant density assumption, because of the difficulty and complexity of obtaining information on the actual distribution of density of topographic masses. Where the actual density of topographic masses deviates from the adopted value, errors are present in the RTM mass-model, and hence, in the forward-modelled residual gravity field. In this paper we attempt to overcome this problem by combining the RTM technique with a high-resolution mass-density model. We compute RTM gravity quantities over New Zealand, with different combinations of elevation models and mass-density assumptions using gravity and GPS/levelling measurements, precise terrain and bathymetry models, a high-resolution mass-density model and constant density assumptions as main input databases. Based on gravity observations and the RTM technique, optimum densities are detected for North Island of ~2500 kg m^?3, South Island of ~2600 kg m^?3, and the whole New Zealand of ~2590 kg m^?3. Comparison among the three sets of residual gravity disturbances computed from different mass-density assumptions show that, together with a global potential model, the high-resolution New Zealand density model explains ~89.5% of gravitational signals, a constant density assumption of 2670 kg m^?3 explains ~90.2%, while a regionally optimum mass-density explains ~90.3%. Detailed comparison shows that the New Zealand density model works best over areas with small residual heights. Over areas with larger residual heights, subsurface density variations appear to affect the residual gravity disturbance. This effect is found to reach about 30 mGal over Southern Alpine Fault. In order to improve the RTM modelling with mass-density maps, a higher-quality mass-density model that provides radially varying mass-density data would be desirable.     

Forests and water in South America

South America is experiencing rapid change in forest cover, of both native and planted forest. Forest cover loss is primarily attributable to fire, logging, and conversion of native forest to agriculture, pasture, and forest plantations, and types of change vary within and among the many diverse types of forests in South America. Major changes in forest cover and growing policy concerns underscore an urgent need for research on sustainable forest management and water ecosystem services in South America. Differences in land ownership and management objectives create trade‐offs between wood production and water ecosystem services from forests. Work is needed to quantify how forest change and management affect ecosystem services, such as wood production versus water provision. Current scientific understanding of forest management effects on water ecosystem services in South America has important limitations, including a scarcity of long‐term records and few long‐term integrated watershed studies. Industry, government, universities, and local communities should collaborate on integrated applied studies of forests and water. Data archiving and publically available data are required. The creation of national networks and a multi‐country South America network to identify and implement common water research protocols, share results, and explore their implications would promote common and well‐supported policies. Hydrologists working in South America are well placed to tackle the challenges and opportunities for collaborative research that will maintain the intrinsic values and water ecosystem services provided by South America's forests.     

Modeling surface water–groundwater interaction in New Zealand: Model development and application

Most rivers worldwide have a strong interaction with groundwater when they leave the mountains and flow over alluvial plains before flowing into the seas or disappearing in the deserts, and in New Zealand, typically, rivers lose water to the groundwater in the upper plains and generally gain water from the groundwater in the lower plains. Aiming at simulating surface water–groundwater interaction nationally in New Zealand, we developed a conceptual groundwater module for the national hydrologic model TopNet to simulate surface water–groundwater interaction, groundwater flow, and intercatchment groundwater flow. The developed model was applied to the Pareora catchment in South Island of New Zealand, where there are concurrent spot gauged flows. Results show that the model simulations not only fit quite well to flow measurement but also to concurrent spot gauged flows, and compared to the original TopNet, it has a significant improvement in the low flows. Sensitivity analysis shows river flow is sensitive to the river losing/gaining rate instead of groundwater characteristic, while groundwater storage is sensitive to both river losing/gaining rate and groundwater characteristic. This indicates our conceptual approach is promising for nationwide modeling without the large amount of geology and aquifer data typically required by physically‐based modeling approaches.     

Assessment of spatial transferability of process‐based hydrological model parameters in two neighbouring catchments in the Himalayan Region

Estimating the hydrological regime of ungauged catchments in the Himalayan region is challenging due to a lack of sufficient monitoring stations. In this paper, the spatial transferability of the model parameters of the process‐oriented J2000 hydrological model was investigated in 2 glaciated subcatchments of the Koshi river basin in eastern Nepal. The catchments have a high degree of similarity with respect to their static landscape features. The model was first calibrated (1986–1991) and validated (1992–1997) in the Dudh Koshi subcatchment. The calibrated and validated model parameters were then transferred to the nearby Tamor catchment (2001–2009). Sensitivity and uncertainty analyses were carried out for both subcatchments to discover the sensitivity range of the parameters in the two catchments. The model represented the overall hydrograph well in both subcatchments, including baseflow, rising and falling limbs; however, the peak flows were underestimated. The efficiency results according to both Nash–Sutcliffe (ENS) and the coefficient of determination (r²) were above 0.84 in both catchments (1986–1997 in Dudh Koshi and 2001–2009 in Tamor). The ranking of the parameters in respect to their sensitivity matched well for both catchments while taking ENS and log Nash–Sutcliffe (LNS) efficiencies into account. However, there were some differences in sensitivity to ENS and LNS for moderately and less‐sensitive parameters, although the majority (13 out of 16 for ENS and 16 out of 16 for LNS) had a sensitivity response in a similar range. The generalized uncertainty likelihood estimation results suggest that the parameter uncertainty are most of the time within the range and the ensemble mean matches very good (ENS: 0.84) with observed discharge. The results indicate that transfer of the J2000 parameters to a neighbouring catchment in the Himalayan region with similar physiographic landscape characteristics is viable. This indicates the possibility of applying a calibrated process‐based J2000 model to other ungauged catchments in the Himalayan region, which could provide important insights into the hydrological system dynamics and provide much needed information to support water resources planning and management.