Stable carbon isotopes (δ<Superscript>13</Superscript>C) of total organic carbon and long-chain <Emphasis Type="Italic">n</Emphasis>-alkanes as proxies for climate and environmental change in a sediment core from Lake Petén-Itzá, Guatemala

Sediment core PI-6 from Lake Petén Itzá, Guatemala, possesses an ~85-ka record of climate and environmental change from lowland Central America. Variations in sediment lithology suggest large and abrupt changes in precipitation during the last glacial and deglacial periods, and into the early Holocene. We measured stable carbon isotope ratios of total organic carbon and long-chain n-alkanes from the core, the latter representing a largely allochthonous (terrestrial) source of organic matter, to reveal past shifts in the relative proportion of C₃–C₄ terrestrial biomass. We sought to test whether stable carbon isotope results were consistent with other paleoclimate proxies measured in the PI-6 core, and if extraction and isotope analysis of n-alkanes is warranted. The largest δ¹³C variations are associated with Heinrich Events. Carbon isotope values in sediments deposited during the last glacial maximum indicate moderate precipitation with little fluctuation. The deglacial was a period of pronounced climate variability, e.g. a relatively warm and moist Bølling–Allerød, but a cool and dry Younger Dryas. Arid periods of the deglacial were inferred from samples with high δ¹³C values in total organic carbon, which reflect times of greater proportions of C₄ plants. These inferences are supported by stable isotope measurements on ostracod shells and relative abundance of grass pollen from the same depths in core PI-6. Similar trends in carbon stable isotopes measured on bulk organic carbon and n-alkanes suggest that carbon isotope measures on bulk organic carbon in sediments from this lake are sufficient to infer past climate-driven shifts in local vegetation.     

Data-Driven Predictive Modeling of Mineral Prospectivity Using Random Forests: A Case Study in Catanduanes Island (Philippines)

The Random Forests (RF) algorithm is a machine learning method that has recently been demonstrated as a viable technique for data-driven predictive modeling of mineral prospectivity, and thus, it is instructive to further examine its usefulness in this particular field. A case study was carried out using data from Catanduanes Island (Philippines) to investigate further (a) if RF modeling can be used for data-driven modeling of mineral prospectivity in areas with few (i.e., <20) mineral occurrences and (b) if RF modeling can handle predictor variables with missing values. We found that RF modeling outperforms evidential belief (EB) modeling of prospectivity for hydrothermal Au–Cu deposits in Catanduanes Island, where 17 hydrothermal Au–Cu prospects are known to exist. Moreover, just like EB modeling, RF modeling allows analysis of the spatial relationships between known prospects and individual layers of predictor data. Furthermore, RF modeling can handle missing values in predictor data through an RF-based imputation technique whereas in EB modeling, missing values are simply represented by maximum uncertainty. Therefore, the RF algorithm is a potentially useful method for data-driven predictive modeling of mineral prospectivity in regions with few (i.e., <20) occurrences of mineral deposits of the type sought. However, further testing of the method in other regions with few mineral occurrences is warranted to fully determine its usefulness in data-driven predictive modeling of mineral prospectivity.     

Characteristics of pegmatoidal granite exposed near Bayalan,Ajmer district,Rajasthan

The study involves the characterization of pegmatoidal granite, southeast of Beawar, Ajmer district, Rajasthan. Earlier researchers had described this granite as part of the BGC, basement to the Bhim Group of the Delhi Super Group rocks. However, the present study indicates that it is younger than the rocks of Bhim Group of South Delhi Fold Belt, into which it is intrusive. The intrusion is structurally controlled and the outcrop pattern is phacolithic. The granite had intruded post-D₂ deformation of the Delhi orogeny along the axial planes of D₂ folds. The intrusion has also resulted in the formation of a contact aureole about the calc gneisses.     

Stochastic evaluation of simple pairing approaches to reconstruct incomplete rainfall time series

Two-station pairing approaches are routinely used to infill missing information in incomplete rainfall databases. We evaluated the performance of three simple methodologies to reconstruct incomplete time series in presence of variable nonlinear correlation between data pairs. Nonlinearity stems from the statistics describing the marginal peak-over-threshold (POT) values of rainfall events. A Monte Carlo analysis was developed to quantitatively assess expected errors from the use of chronological pairing (CP) with linear and nonlinear regression and frequency pairing (FP). CP is based on a priori selection of regression functions, while FP is based on matching the probability of non-exceedance of an event from one time series with the probability of non-exceedance of a similar event from another time series. We adopted a generalized Pareto (GP) model to describe POT events, and a t-copula algorithm to generate reference nonlinearly correlated pairs of random temporal distributions distributed according with the GP model. The results suggest that the optimal methodology strongly depends on GP statistics. In general, CP seems to provide the lowest errors when GP statistics were similar and correlation became linear; we found that a power-2 function performs well for the selected statistics when the number of missing points is limited. FP outperforms the other methods when POT statistics are different and variables are markedly nonlinearly correlated. Ensemble-based results seem to be supported by the analysis of observed precipitation at two real-world gauge stations.     

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.     

ARMA order model detection using minimum of Kurtosis: application on seismic data

The objective of this study is to find the order and coefficients of non-low-phase causal filters for ARMA (auto regressive moving average) filter model, using the Kurtosis minimization criterion. This method is based on the Kurtosis calculation of the treated sample at the input level and its identification at the output of the ARMA model. For this purpose, the order and coefficients of the AR (auto regressive) part are identified using the Yule-Walker algorithm at order two and then extended to order four. To obtain the MA (moving average) part, the AR components are calculated at first from the ARMA filter by deconvolution. Then, spectrally equivalent and minimum phase (SEMP) MA filter is identified using the Durbin algorithm at second and fourth order. Finally, the correct filter is found when the Kurtosis value of the output ARMA filter reconstituted is the closest to the Kurtosis of introduced signal. The proposed method is then tested on simulated processes and applied to real seismic data to perform blind deconvolution and obtain the reflectivity coefficients of subsoil studied.     

Phytoplankton biomass and size fractions in surface waters of the Australian sector of the southern ocean

We collected surface water along the 142nd E meridian from Tasmania to Antarctica in December 1999. We measured temperature, salinity and total chlorophyll a; additionally, we collected suspended particle size fractions and used fluorometric analysis to determine the quantity of chlorophyll a in each of four cell size classes: picoplankton (<3 μm), two nanoplankton fractions (3–10 μm and 10–20 μm) and microplankton (> 20 μm). Changes in temperature and salinity show that we crossed 6 water masses separated by 5 fronts. We found low abundance (<0.2 mg m⁻³) of chlorophyll in all size classes, with the exception of higher values near the continent (0.2 to 0.4 mg m⁻³). Lowest chlorophyll values (<0.1 mg m⁻³) were found in the Polar Frontal Zone (51° to 54°S). Microplankton made up the largest portion of total chlorophyll throughout most of the region. We conclude that biomass of all phytoplankton fractions, especially pico-and nanoplankton, was constrained by limiting factors, most probably iron, throughout the region and that ecosystem dynamics within a zone are not circumpolar but are regionalized within sectors.     

Parametric vibration of submerged floating tunnel tether under random excitation

For the study of the parametric vibration response of submerged floating tunnel tether under random excitation,a nonlinear random parametric vibration equation of coupled tether and tube of submerged floating tunnel is set up.Subsequently,vibration response of tether in the tether-tube system is analyzed by Monte Carlo method.It may be concluded that when the tube is subjected to zero-mean Gaussian white noise random excitation,the displacement and velocity root mean square responses of tether reach the peak if the circular frequency of tube doubles that of tether;the displacement and velocity root mean square responses of tether increase as the random excitation root mean square increases;owing to the damping force of water,the displacement and velocity root mean square responses of tether decrease rapidly compared with tether in air;increasing the damping of the tether or tube reduces the displacement and velocity root mean square responses of tether;the large-amplitude vibration of tether may be avoided by locating dampers on the tether or tube.     

Study on scattering wave force of horizontal and vertical plate type breakwaters

The interaction between wave and horizontal and vertical plates is investigated by the boundary element method,and the relations of wave exciting force with plate thickness,submergence and length are obtained.It is found that:1) The efficient wave exciting force exists while plate submergence is less than 0.5 m,and the plate is very thin with order O(0.005 m).2) The maximum heave wave exciting force exists,and it is the main factor for surface and submerged horizontal plate while the roll force can be ignored.3) The maximum sway wave exciting force exists,it is the main factor for surface or submerged vertical plate,and the roll force is about 20 times of horizontal plate.     

SAR-imaged spiral eddies in Mutsu Bay and their dynamic and kinematic models

A pattern of slick streaks winding into a spiral, known as a spiral eddy, was identified in 5 images taken by the ERS-1/2 synthetic aperture radar (SAR) in Mutsu Bay (Japan); dynamic and kinematic models of these spiral eddies have been proposed. Common characteristics of the five spiral eddies are: 1) an eddy diameter of about 15 km; 2) their location in the western part of the bay; and 3) their cyclonic direction of rotation. Moreover, the wind conditions over the bay were common: prior to acquiring the images, a strong easterly wind continued blowing for more than one day. The wind field on the bay is known to be orographically steered and has strong windstress vorticity, which generates cyclonic circulation. The diameter and location of the circulation simulated with a numerical ocean model corresponded well to those of the identified spiral eddies. Based on these facts, we propose a dynamic model for the movement of a slick streak, and a kinematic model for the formation of a spiral eddy. We have assumed calm air, a microlayer and seawater with a cyclonic circulation in the dynamic model. The balance of forces is established in the microlayer among the frictional force from the seawater, the frictional force from the calm air, the gravitational force, and the Coriolis force. As a result, the velocity vector of the microlayer deflects slightly towards the center of the cyclonic circulation. We have assumed a point source of the microlayer in the kinematic model. The shapes of a slick streak simulated with the models agree well with the identified patterns in the SAR images.