Calculated sea surface temperatures of two sediment cores

The reconstruction of ocean history employs a large variety of methods with origins in the biological, chemical, and physical sciences, and uses modern statistical techniques for the interpretation of extensive and complex data sets. Various sediment properties deliver useful information for reconstructing environmental parameters. Those properties that have a close relationship to environmental parameters are called ''proxy variables'' (''proxies'' for short). Proxies are measurable descriptors for desired (but unobservable) variables. Surface water temperature is probably the most important parameter for describing the conditions of past oceans and is crucial for climate modelling. Proxies for temperature are: abundance of microfossils dwelling in surface waters, oxygen isotope composition of planktic foraminifers, the ratio of magnesium or strontium to calcium in calcareous shells or the ratio of certain organic molecules (e.g. alkenones produced by coccolithophorids). Surface water salinity, which is important in modelling of ocean circulation, is much more difficult to reconstruct. At present there is no established method for a direct determination of this parameter. Measurements associated with the paleochemistry of bottom waters to reconstruct bottom water age and flow are made on benthic foraminifers, ostracodes, and deep-sea corals. Important geochemical tracers are d13C and Cd/Ca ratios. When using benthic foraminifers, knowledge of the sediment depth habitat of species is crucial. Reconstructions of productivity patterns are of great interest because of important links to current patterns, mixing of water masses, wind, the global carbon cycle, and biogeography. Productivity is reflected in the flux of carbon into the sediment. There are a number of fluxes other than those of organic carbon that can be useful in assessing productivity fluctuations. Among others, carbonate and opal flux have been used, as well as particulate barite. Furthermore, microfossil assemblages contain clues to the intensity of production as some species occur preferentially in high-productivity regions while others avoid these. One marker for the fertility of sub-surface waters (that is, nutrient availability) is the carbon isotope ratio within that water (13C/12C, expressed as d13C). Carbon isotope ratios in today's ocean are negatively correlated with nitrate and phosphate contents. Another tracer of phosphate content in ocean waters is the Cd/Ca ratio. The correlation between this ratio and phosphate concentrations is quite well documented. A rather new development to obtain clues on ocean fertility (nitrate utilization) is the analysis of the 15N/14N ratio in organic matter. The fractionation dynamics are analogous to those of carbon isotopes. These various ratios are captured within the organisms growing within the tagged water. A number of reconstructions of the partial pressure of CO2 have been attempted using d13C differences between planktic and benthic foraminifers and d13C values of bulk organic material or individual organic components. To define the carbon system in sea water, two elements of the system have to be known in addition to temperature. These can be any combination of total CO2 , alkalinity, or pH. To reconstruct pH, the boron isotope composition of carbonates has been used. Ba patterns have been used to infer the distribution of alkalinity in past oceans. Information relating to atmospheric circulationand climate is transported to the ocean by wind or rivers, in the form of minerals or as plant andanimal remains. The most useful tracers in this respect are silt-sized particles and pollen.

Numerical Investigation of the Flow Features around Heave Plates Oscillating Close to a Free Surface or Seabed

Performance of heave plates used in offshore structures is strongly influenced by their added mass and damping, which are affected by proximity to a boundary. A previous paper by the authors presented numerical simulations of the flow around a circular solid disk oscillating at varying elevations from seabed [1]. The force calculated was used to evaluate the added mass and damping coefficients for the disk. The simulations suggest that as the structure moves closer to the seabed the added mass and damping coefficients (Ca and Cb) increases continuously. In order to understand the physics behind the added mass and damping trends, when a heave plate is moving near a seabed or closer to the free surface, the flow characteristics around the heave plate are examined numerically in this paper. Flow around oscillating disks is dominated by generation and development of phase-dependent vortical structures, characterized by the KC number and the distance from the seabed or free surface to the heave plate. Numerical calculations presented in this paper have comprised the qualitative analysis of the vortex shedding and the investigation of the links between such vortex shedding and, on one hand the damping coefficient, and on the other hand, pairing mechanisms such as the shedding angle.

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated. (c) 2004 Elsevier Ltd. All rights reserved.

Thermal conductivity of ODP Leg 110 holes

Thirty-four sediment and mudline temperatures were collected from six drill holes on ODP Leg 110 near the toe of the Barbados accretionary complex. When combined with thermal conductivity measurements these data delineate the complicated thermal structure on the edge of this convergent margin. Surface heat-flow values from Leg 110 (calculated from geothermal gradients forced through the bottom-water temperature at mudline) of 92 to 192 mW/m**2 are 80% to 300% higher than values predicted by standard heat flow vs. age models for oceanic crust, but are compatible with earlier surface measurements made at the same latitude. Measured heat flow tends to decrease downhole at four sites, suggesting the presence of heat sources within the sediments. These results are consistent with the flow of warm fluid through the complex along sub-horizontal, high-permeability conduits, including thrust faults, the major decollement zone, and sandy intervals. Simple calculations suggest that this flow is transient, occurring on time scales of tens to tens of thousands of years. High heat flow in the vicinity of 15°30'N and not elsewhere along the deformation front suggests that the Leg 110 drill sites may be situated over a fluid discharge zone, with dewatering more active here than elsewhere along the accretionary complex.

Air Entrainment in the Turbulent Ship Hull Boundary Layer

Turbulent fluctuations in the vicinity of the water free surface along a flat, vertically oriented surface-piercing plate are studied experimentally using a laboratory-scale experiment. In this experiment, a meter-wide stainless steel belt travels horizontally in a loop around two rollers with vertically oriented axes, which are separated by 7.5 meters. This belt device is mounted inside a large water tank with the water level set just below the top edge of the belt. The belt, rollers, and supporting frame are contained within a sheet metal box to keep the device dry except for one 6-meter-long straight test section between rollers. The belt is launched from rest with an acceleration of up to 3-g in order to quickly reach steady state velocity. This creates a temporally evolving boundary layer analogous to the spatially evolving boundary layer created along a flat-sided ship moving at the same velocity, with a length equivalent to the length of belt that has passed the measurement region since the belt motion began. Surface profile measurements in planes normal to the belt surface are conducted using cinematic Laser Induced Fluorescence and quantitative surface profiles are extracted at each instant in time. Using these measurements, free surface fluctuations are examined and the propagation behavior of these free surface ripples is studied. It is found that free surface fluctuations are generated in a region close to the belt surface, where sub-surface velocity fluctuations influence the behavior of these free surface features. These rapidly-changing surface features close to the belt appear to lead to the generation of freely-propagating waves far from the belt, outside the influence of the boundary layer. Sub-surface PIV measurements are performed in order to study the modification of the boundary layer flow field due to the effects of the water free surface. Cinematic planar PIV measurements are performed in horizontal planes parallel to the free surface by imaging the flow from underneath the tank, providing streamwise and wall-normal velocity fields. Additional planar PIV experiments are performed in vertical planes parallel to the belt surface in order to study the bahvior of streamwise and vertical velocity fields. It is found that the boundary layer grows rapidly near the free surface, leading to an overall thicker boundary layer close to the surface. This rapid boundary layer growth appears to be linked to a process of free surface bursting, the sudden onset of free surface fluctuations. Cinematic white light movies are recorded from beneath the water surface in order to determine the onset location of air entrainment. In addition, qualitative observations of these processes are made in order to determine the mechanisms leading to air entrainment present in this flow.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Geomodelling of a fluvial system with semi-supervised support vector regression

Fluvial deposits are a challenge for modelling flow in sub-surface reservoirs. Connectivity and continuity of permeable bodies have a major impact on fluid flow in porous media. Contemporary object-based and multipoint statistics methods face a problem of robust representation of connected structures. An alternative approach to model petrophysical properties is based on machine learning algorithm ? Support Vector Regression (SVR). Semi-supervised SVR is able to establish spatial connectivity taking into account the prior knowledge on natural similarities. SVR as a learning algorithm is robust to noise and captures dependencies from all available data. Semi-supervised SVR applied to a synthetic fluvial reservoir demonstrated robust results, which are well matched to the flow performance

Horizontal Grid Size Selection and its Influence on Mesoscale Model Simulations

The use of two-dimensional spectral analysis applied to terrain heights in order to determine characteristic terrain spatial scales and its subsequent use for the objective definition of an adequate grid size required to resolve terrain forcing are presented in this paper. In order to illustrate the influence of grid size, atmospheric flow in a complex terrain area of the Spanish east coast is simulated by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical model using different horizontal grid resolutions. In this area, a grid size of 2 km is required to account for 95% of terrain variance. Comparison among results of the different simulations shows that, although the main wind behavior does not change dramatically, some small-scale features appear when using a resolution of 2 km or finer. Horizontal flow pattern differences are significant both in the nighttime, when terrain forcing is more relevant, and in the daytime, when thermal forcing is dominant. Vertical structures also are investigated, and results show that vertical advection is influenced highly by the horizontal grid size during the daytime period. The turbulent kinetic energy and potential temperature vertical cross sections show substantial differences in the structure of the planetary boundary layer for each model configuration

Investigation of the occurrence of pesticide residues in rural wells and surface water following application to tobacco

In this work the exposure of wells and surface water to pesticides, commonly used for tobacco cropping, was assessed. Water consumption wells and surface water flows were sampled at different times. After a preconcentration step with solid phase extraction (SPE), the selected pesticides were determined by gas chromatography with electron capture detection (GC-ECD) or high performance liquid chromatography with diode array detection (HPLC-DAD). No pesticides were detected in the well water samples and surface water flow in the winter season. However, in the spring and summer higher concentrations of chlorpyrifos and imidacloprid were found in the water source samples. Atrazine, simazine and clomazone were also found. The occurrence of pesticides in collected water samples was related with the application to tobacco.

Two-Phase (Solid-Liquid) Flow in Inclined Pipes

This paper presents an experimental research about the behavior of two-phase flows in inclined pipes. The inclination angle varied from 5° to 45° and the slurry solid concentration varied up to 15%. It was concluded that the head losses of the downward sloping pipe flow are always lower than the head losses of the horizontal flow and these are always lower than the head losses of the upward sloping pipe flow, regardless the concentration and inclination angle. It was possible to develop empirical equations to calculate the head losses of the horizontal flow and the upward and downward sloping pipe flows.

A Coupled 2 × 2D Babcock–Leighton Solar Dynamo Model. I. Surface Magnetic Flux Evolution

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.

An integrated study on the hydrogeology of Bharathapuzha river basin, South West coast of India

The present work deals with the An integrated study on the hydrogeology of Bharathapuzha river basin ,south west coast of india. To study the spatial and temporal behaviour of the groundwater system of the Bharathapuzha river basin.To discover the sub-surface parameter by ground resistivity surveys.T o determine the groundwater quality of the Bharathapuzha river basin for the different seasons {pre monsoon and post monsoon with reference to the domestic and irrigational water quality standards.Present study will provide a good database on the hydrogeological aspects within the river basin.The study area covers l7 block Panchayats. Of these, Chitoor block is ‘over exploited’, Kollengode, Trithala, and Palakkad are ‘critical’ in category and Kuttippuram and Sreekrishnapuram blocks are ‘semi critical’ in terms of groundwater development.Comparison of Geomorphology map with drainage map shows that the geomorphology has a clear control on the drainage net work of the basin. The structural hill area shows a highest drainage network, where as pediment shows lowest drainage network.There are many discontinuous lineament in the Bharathapuzha river basin which can be connected by a straight line.Ground water flow directions are generally towards the western portions of the study area. From the northern region Water flows towards the central and also water from the eastern and southern side confluences at the centre and move towards western side of the basin.The positive correlation of transmissivity and storativity values show good aquifer conditions exists in the present study area .

Sensitivity of radar backscatter to desert surface roughness

Synthetic aperture radar (SAR) data have proved useful in remote sensing studies of deserts, enabling different surfaces to be discriminated by differences in roughness properties. Roughness is characterized in SAR backscatter models using the standard deviation of surface heights (sigma), correlation length (L) and autocorrelation function (rho(xi)). Previous research has suggested that these parameters are of limited use for characterizing surface roughness, and are often unreliable due to the collection of too few roughness profiles, or under-sampling in terms of resolution or profile length (L-p). This paper reports on work aimed at establishing the effects of L-p and sampling resolution on SAR backscatter estimations and site discrimination. Results indicate significant relationships between the average roughness parameters and L-p, but large variability in roughness parameters prevents any clear understanding of these relationships. Integral equation model simulations demonstrate limited change with L-p and under-estimate backscatter relative to SAR observations. However, modelled and observed backscatter conform in pattern and magnitude for C-band systems but not for L-band data. Variation in surface roughness alone does not explain variability in site discrimination. Other factors (possibly sub-surface scattering) appear to play a significant role in controlling backscatter characteristics at lower frequencies.

Evaluating the potential for statistical decadal predictions of sea surface temperatures with a perfect model approach

We explore the potential for making statistical decadal predictions of sea surface temperatures (SSTs) in a perfect model analysis, with a focus on the Atlantic basin. Various statistical methods (Lagged correlations, Linear Inverse Modelling and Constructed Analogue) are found to have significant skill in predicting the internal variability of Atlantic SSTs for up to a decade ahead in control integrations of two different global climate models (GCMs), namely HadCM3 and HadGEM1. Statistical methods which consider non-local information tend to perform best, but which is the most successful statistical method depends on the region considered, GCM data used and prediction lead time. However, the Constructed Analogue method tends to have the highest skill at longer lead times. Importantly, the regions of greatest prediction skill can be very different to regions identified as potentially predictable from variance explained arguments. This finding suggests that significant local decadal variability is not necessarily a prerequisite for skillful decadal predictions, and that the statistical methods are capturing some of the dynamics of low-frequency SST evolution. In particular, using data from HadGEM1, significant skill at lead times of 6–10 years is found in the tropical North Atlantic, a region with relatively little decadal variability compared to interannual variability. This skill appears to come from reconstructing the SSTs in the far north Atlantic, suggesting that the more northern latitudes are optimal for SST observations to improve predictions. We additionally explore whether adding sub-surface temperature data improves these decadal statistical predictions, and find that, again, it depends on the region, prediction lead time and GCM data used. Overall, we argue that the estimated prediction skill motivates the further development of statistical decadal predictions of SSTs as a benchmark for current and future GCM-based decadal climate predictions.

The basic ingredients of the North Atlantic Storm Track. Part II: Sea surface temperatures

The impact of North Atlantic SST patterns on the storm track is investigated using a hierarchy of GCM simulations using idealized (aquaplanet) and “semirealistic” boundary conditions in the atmospheric component (HadAM3) of the third climate configuration of the Met Office Unified Model (HadCM3). This framework enables the mechanisms determining the tropospheric response to North Atlantic SST patterns to be examined, both in isolation and in combination with continental-scale landmasses and orography. In isolation, a “Gulf Stream” SST pattern acts to strengthen the downstream storm track while a “North Atlantic Drift” SST pattern weakens it. These changes are consistent with changes in the extratropical SST gradient and near-surface baroclinicity, and each storm-track response is associated with a consistent change in the tropospheric jet structure. Locally enhanced near-surface horizontal wind convergence is found over the warm side of strengthened SST gradients associated with ascending air and increased precipitation, consistent with previous studies. When the combined SST pattern is introduced into the semirealistic framework (including the “North American” continent and the “Rocky Mountains”), the results suggest that the topographically generated southwest–northeast tilt in the North Atlantic storm track is enhanced. In particular, the Gulf Stream shifts the storm track south in the western Atlantic whereas the strong high-latitude SST gradient in the northeastern Atlantic enhances the storm track there.