A mass transfer model of bauxite formation

The formation of bauxite due to weathering of a granitic protolith has been simulated by means of a one-dimensional flow and reaction model based on the mass transfer principle. The model couples mineral dissolution and precipitation reactions, speciation in solution, and advective solute transport in a porous medium. A very important aspect of the modeling study is the use of mineral reaction rates determined experimentally in the laboratory. The important effects of solution saturation state and pH have been incorporated into the kinetic rate laws governing the heterogeneous reactions. The values of these parameters have been obtained from the scientific literature to guarantee that realistic reaction rates are used in the simulations. Albite and quartz are the minerals that make up the parent rock in the model. Gibbsite, kaolinite, and a Na-mica (as a surrogate for smectite) are the secondary minerals that have been taken into account. Long-term simulations (>1 Ma) have been run, and the formation of a bauxitic profile, with an upper gibbsite-rich and a lower kaolinite-rich zone, is predicted. In early stages of the process (up to a few hundreds of thousands of years), both gibbsite and kaolinite precipitate directly from solution as a consequence of albite dissolution. In later stages, the bulk of gibbsite precipitation derives from the incongruent dissolution of kaolinite, while kaolinite precipitation is still caused by the dissolution of albite. This is also reflected by the formation of two reaction fronts in the profile. These results are compared with weathering sequences from the Los Pijiguaos bauxite deposit, Venezuela. The overlap between the gibbsite and kaolinite zones and the replacement of kaolinite by gibbsite are consistent with model calculations. Mechanical denudation has to be called upon to explain the limited thicknesses of the weathering profiles in the field. The role of mechanical erosion is supported by the presence of microsedimentary structures in the bauxite and the balance between dissolved and suspended loads in the streams draining the area.     

Carbonate platform depositional sequences, Upper Cretaceous, south-central Pyrenees (Spain)

The Upper Cretaceous (Cenomanian-Maastrichtian) of the south-central Pyrenees shows five carbonate platform sequences where the major parameters are tectonism, relative sea-level fluctuations and inherited depositional profile. Depositional geometries and basin analysis permit an understanding of the depositional history.Five depositional sequences have been recognized: (1) The Santa Fe sequence (Middle-Upper Cenomanian), a ramp to a skeletal rimmed shelf with an escarpment bypass margin. The lower boundary is an angular unconformity and the upper one records a sea-level drop. The platform location of the margin was determined by a listric normal fault. (2) An abrupt sea-level rise drowned the former platform. The Congost sequence (Turonian-Lower Coniacian), a distally steepened ramp with erosional distal deep slope. The depositional model was largely controlled by pre-existing basin morphology. Cessation of platform development was due to a relative sea-level drop. (3) The Sant Corneli sequence (Upper Coniacian-Lower Santonian), a mixed terrigenous-skeletal homoclinal ramp with upright margin, deep slope and dysaerobic basin. The slope results from the backstepping by 24 km of the previous margin and gentle basin tilting. The platform margin remained more or less at the same position, and relief between platform and slope increased indicating continuous relative sea-level rise. The upper boundary is an angular unconformity at the platform margin produced by an abrupt sea-level rise and platform drowning, and by listric normal faulting. (4) The Vallcarga sequence (Upper Santonian-Campanian), a distal-steepened skeletal homoclinal ramp, erosional escarpment and turbidite basin, which corresponds to the Mesozoic maximum marine expansion. A listric normal fault created two depositional areas: a more or less flat footwall block with a north-northwest prograding carbonate ramp.     

Partial Melting and Counterclockwise P T Path of Subducted Oceanic Crust (Sierra del Convento Melange, Cuba)

Partial melting of subducted oceanic crust has been identifiedin the Sierra del Convento mélange (Cuba). This serpentinite-matrixmélange contains blocks of mid-ocean ridge basalt (MORB)-derivedplagioclase-lacking epidote ± garnet amphibolite intimatelyassociated with peraluminous trondhjemitic–tonalitic rocks.Field relations, major element bulk-rock compositions, mineralassemblages, peak metamorphic conditions (c. 750°C, 14–16kbar), experimental evidence, and theoretical phase relationssupport formation of the trondhjemitic–tonalitic rocksby wet melting of subducted amphibolites. Phase relations andmass-balance calculations indicate eutectic- and peritectic-likemelting reactions characterized by large stoichiometric coefficientsof reactant plagioclase and suggest that this phase was completelyconsumed upon melting. The magmatic assemblages of the trondhjemitic–tonaliticmelts, consisting of plagioclase, quartz, epidote, ±paragonite, ± pargasite, and ± kyanite, crystallizedat depth (14–15 kbar). The peraluminous composition ofthe melts is consistent with experimental evidence, explainsthe presence of magmatic paragonite and (relict) kyanite, andplaces important constraints on the interpretation of slab-derivedmagmatic rocks. Calculated P–T conditions indicate counterclockwiseP–T paths during exhumation, when retrograde blueschist-faciesoverprints, composed of combinations of omphacite, glaucophane,actinolite, tremolite, paragonite, lawsonite, albite, (clino)zoisite,chlorite, pumpellyite and phengite, were formed in the amphibolitesand trondhjemites. Partial melting of subducted oceanic crustin eastern Cuba is unique in the Caribbean realm and has importantconsequences for the plate-tectonic interpretation of the region,as it supports a scenario of onset of subduction of a youngoceanic lithosphere during the early Cretaceous (c. 120 Ma).The counterclockwise P–T paths were caused by ensuingexhumation during continued subduction. KEY WORDS: amphibolite; Cuba; exhumation; partial melting; trondhjemite; subduction     

A Monte Carlo Model Of The Nocturnal Surface Temperatures In Urban Canyons

A model for the urban canyon is formulated for meteorologicalconditions of weak winds at night time. Thermal radiation, conductivity and convection are simulated by means of the Monte Carlo method. These are the main physical processesof energy transfer that give rise to the characteristic temperaturedistribution in these systems. The model has been satisfactory tested under ideal conditions for which analytical solutions exist.The predictions of the model under morerealistic conditions accurately reproduce the observationalresults. A strong temperature gradient across streets, with the canyon corners up to 4 °C warmer than the canyon centre, is found for the deepest canyons. This theoretical predictionhas been successfully verified with measurementstaken in a number of streets of the city of Granada in Spain.     

Impact of anthropogenic CO 2 on the next glacial cycle

The model of Paillard and Parrenin (Earth Planet Sci Lett 227(3–4):263–271, 2004) has been recently optimized for the last eight glacial cycles, leading to two different relaxation models with model-data correlations between 0.8 and 0.9 (García-Olivares and Herrero (Clim Dyn 1–25, 2012b)). These two models are here used to predict the effect of an anthropogenic CO ₂ pulse on the evolution of atmospheric CO ₂, global ice volume and Antarctic ice cover during the next 300 kyr. The initial atmospheric CO ₂ condition is obtained after a critical data analysis that sets 1300 Gt as the most realistic carbon Ultimate Recoverable Resources (URR), with the help of a global compartmental model to determine the carbon transfer function to the atmosphere. The next 20 kyr will have an abnormally high greenhouse effect which, according to the CO ₂ values, will lengthen the present interglacial by some 25 to 33 kyr. This is because the perturbation of the current interglacial will lead to a delay in the future advance of the ice sheet on the Antarctic shelf, causing that the relative maximum of boreal insolation found 65 kyr after present (AP) will not affect the developing glaciation. Instead, it will be the following insolation peak, about 110 kyr AP, which will find an appropriate climatic state to trigger the next deglaciation.     

Influence of ENSO and of the Indian Ocean Dipole on the Indian summer monsoon variability

Indian summer monsoon (ISM) variability is forced from external factors (like the El Niño Southern Oscillation, ENSO) but it contains also an internal component that tends to reduce its potential for predictability. Large-scale and local monsoon indices based on precipitation and atmospheric circulation parameters are used as a measure of ISM variability. In a 9-members ensemble of AMIP-type experiments (with same boundary SST forcing and different initial conditions) their potential predictability is comparable using both local and large-scale monsoon indices. In the sample analyzed, about half of more predictable monsoon years coincide with El Niño and/or positive Indian Ocean Dipole (IOD) events. Summer monsoon characteristics during ENSO and IOD years are analyzed through composites computed over a three years period (i.e. one year before and one year after the event peak) to investigate the mutual relationship between the events lagged in time. The connection between ISM and IOD is mostly confined in the summer and autumn, while that with ENSO is stronger and extends more in time. In the coupled model results the IOD influence on the monsoon is large, even because in the model IOD events are intense and easily reproduced due to a strong air-sea feedback in the eastern side of the basin. Monsoon seasons preceding or following an El Niño or a La Niña event are not exactly symmetric, even in terms of their biennial character. In most of the cases, both in reanalysis and model, El Niño and positive IOD events tend to co-occur with larger anomalies either in the Indo-Pacific ocean sector or over India, while La Niña and negative IOD do not. From the observed record, the ENSO-IOD correlation is positive strong and significant since mid-60s and it may correspond with either strong or weak ENSO-monsoon relationship and with strong or weak IOD-monsoon relationship. A main difference between those periods is the relationship between Indian monsoon rainfall and SST in other ocean basins rather than the Indo-Pacific sector alone.     

Simulation of discrete cracks driven by nearly incompressible fluid via 2D combined finite-discrete element method

In this work, we propose a novel hydraulic solver in order to simulate key mechanisms that control fluid-driven cracks in the framework of the combined finite-discrete element method (FDEM). The main innovative aspect of the present work is the independence of the fluid's critical time step size with the fracture opening. This advantage is extremely important because it means that very fine meshes can be used around areas of interest, such as boreholes, without penalizing the computational cost as fractures propagate (ie, open) and the fluid flows through them. This is a great advantage over other recently introduced approaches that exhibit a dependency of the time step in the form of Δt_crit ∝ (l/a)² where l is the element size and a is the fracture aperture. This paper presents a series of benchmark cases for the proposed solver. The rationale adopted by the authors was to benchmark and validate the implementation of the hydraulic solver in an incremental fashion, starting from the simplest cases and building in complexity. The results shown in this work clearly demonstrate that the proposed approach is able to reproduce analytical results for fluid flow through a single crack. The results presented in this paper also demonstrate that the new approach is robust enough to deal with complex fracture patterns and complex geometries; the obtained fluid-driven fracture patterns in the vicinity of a borehole certainly stand to the scrutiny of human visual perception.     

Climatic variations in the coastal zone — Comparison between the Po River delta (Adriatic Sea,Italy) and the Huanghe River delta (Bohai Sea,China)

The aim of this work is to research the influence of natural climatic changes on the evolution of the coastal zone in modern times and the possible implication of human activities on the configuration of the present coastline. Comparison of data of two very far and different areas, the Po River delta, Adriatic Sea and the Huanghe River delta, Bohai Sea, reveals the planetary diffusion of climatic fluctuations and their effects on coastal evolution. Contribution No. 820 from the Marine Geology Institute, C. N. R., Italy. Contribution No.2207 from the Institute of Oceanology, Academia Sinica. Project supported by the National Natural Science Foundation of China.