Modelling of Ocean Shallow Convection under Fast Ice in the Arctic

Ocean convection in the Antarctic has been studied many times and has been revealed to be responsible for ice-cover reduction. In the Arctic, proof of that phenomenon has not been documented. It is believed that this phenomenon happens on a smaller scale in the Arctic when local circulation of deep warmer water melts and slows ice production. An example of this is the North Water (NOW) polynya in northern Baffin Bay. A polynya is an area of open water in an otherwise ice-covered area. As ice forms under the fast ice near the boundary of the polynya, ocean salts (brine) are ejected from the newly formed ice. This water, which has an increased concentration of salt, sinks and is replaced by warmer water from below, and this slows ice formation. In our study a coupled one-dimensional thermodynamic snow–fast ice model incorporating ocean heat flux input via a shallow convection model was used. Ice thickness was calculated using a thermodynamic model that included a current-induced entrainment model and a convection model to account for brine rejection during ice growth. Atmospheric observations from Grise Fiord and Thule and ocean profiles around the NOW polynya near these sites were used as input to the model. This purely thermodynamic study enables us to obtain ice thickness values that can be compared with qualitative observations. This modelling study compares two sites related to the NOW polynya. The results indicate that the shallow convection model simulates the reduction of fast ice near Thule but not near Grise Fiord.     

Discharge and suspended sediment flux estimated along the mainstream of the Amazon and the Madeira Rivers (from in situ and MODIS Satellite Data)

Water and suspended sediment fluxes are considered during the period 2000–2008 in a region including the full Amazon River from the confluence of the Negro River to Santarém, the end part of the Solimões River, and the lower part of the Madeira River. Three types of data are used: water discharge estimated from field measurements, and suspended sediment obtained from field measurements and derived from MODIS satellite data. A generalized least square method including a propagating term is developed in order to propagate the signal upward and downward the river. The approach is introduced and tested. Several experiments are considered in order, first, to estimate the ability to propagate the signal from stations located before the confluences of Negro and Madeira Rivers to stations located on the Amazon River; second to investigate the possibility to propagate the signal along the Amazon River which dynamics is coupled with floodplains dynamics; and third to produce optimal solutions of water and sediment fluxes. For each experiment, the influence of field and satellite data is compared. The approach is efficient in the upper part of the region of study where the Solimões, the Negro and the Madeira Rivers meet and fails in the lower part of the region where interactions between Amazon River and floodplains play an important role on the fluxes’ dynamics. The optimal experiment includes in situ and satellite data from all the stations available and is used to analyse the recent evolution of suspended sediment flux along the Amazon River and its interaction with the large coupled floodplains. A high accumulation rate is observed during the 2000–2002 period, followed by decreasing rates until 2005 and by increasing values in 2006 and 2007. Our results suggest that floodplains extending along a river reach of 390 km-long between Itacoatiara and Óbidos trap about 15% of the suspended sediment flux passing at Óbidos. The simulated deposition rate is of about 0.3 Mt km⁻¹ yr⁻¹ corresponding to an accretion rate of about 27 mm yr⁻¹.     

Eclogite thermobarometry: The consistency between conventional thermobarometry and forward phase-equilibrium modelling

Eclogite thermobarometry is crucial for constraining the depths and temperatures to which oceanic and continental crust subduct. However, obtaining the pressure and temperature (P–T) conditions of eclogites is complex as they commonly display high-variance mineral assemblages, and the mineral compositions only vary slightly with P–T. In this contribution, we present a comparison between two independent and commonly used thermobarometric approaches for eclogites: conventional thermobarometry and forward phase-equilibrium modelling. We assess how consistent the thermobarometric calculations are using the garnet–clinopyroxene–phengite barometer and garnet–clinopyroxene thermometer with predictions from forward modelling (i.e. comparing the relative differences between approaches). Our results show that the overall mismatch in methods is typically ±0.2–0.3 GPa and ±29–42°C although differences as large as 80°C and 0.7 GPa are possible for a few narrow ranges of P–T conditions in the forward models. Such mismatch is interpreted as the relative differences among methods, and not as absolute uncertainties or accuracies for either method. For most of the investigated P–T conditions, the relatively minor differences between methods means that the choice in thermobarometric method itself is less important for geological interpretation than careful sample characterization and petrographic interpretation for deriving P–T from eclogites. Although thermobarometry is known to be sensitive to the assumed X_Fe³⁺ of a rock (or mineral), the relative differences between methods are not particularly sensitive to the choice of bulk-rock X_Fe³⁺, except at high temperatures (>650°C, amphibole absent) and for very large differences in assumed X_Fe³⁺ (0–0.5). We find that the most important difference between approaches is the activity–composition (a–x) relations, as opposed to the end-member thermodynamic data or other aspects of experimental calibration. When equivalent a–x relations are used in the conventional barometer, P calculations are nearly identical to phase-equilibrium models (ΔP < 0.1). To further assess the implications of these results for real rocks, we also evaluate common mathematical optimizations of reaction constants used for obtaining the maximum P–T with conventional thermobarometric approaches (e.g. using the highest aGrs² × aPrp in garnet and Si content in phengite, and the lowest aDi in clinopyroxene). These approaches should be used with caution, because they may not represent the compositions of equilibrium mineral assemblages at eclogite facies conditions and therefore systematically bias P–T calculations. Assuming method accuracy, geological meaningful P_max at a typical eclogite facies temperature of ~660°C will be obtained by using the greatest aDi, aCel, and aPrp and lowest aGrs and aMs; garnet and clinopyroxene with the lowest Fe²⁺/Mg ratios may yield geological meaningful T_max at a typical eclogite facies pressure of 2.5 GPa.     

Flexural and gravity modelling of the Mérida Andes and Barinas–Apure Basin, Western Venezuela

The kinematic evolution of the Barinas–Apure Basin and the southern Mérida Andes from Lower Miocene to the Present is numerically modelled using flexural isostatic theory and geophysical and geological data. Two published regional transects are used to build up a reference section, which is then used to constrain important parameters (e.g. shortenings and sedimentary thicknesses) for the flexural modelling. To control the location of the main fault system in the flexural model earthquake information is also used. The estimated flexural elastic thickness of the South American lithosphere beneath the Barinas–Apure Basin and the Mérida Andes Range is 25 km. The value for the final total shortening is 60 km. The flexural isostatic model shows that the Andean uplift has caused the South American lithosphere subsidence and the development of the Barinas–Apure Basin.In addition, gravity modelling was used to understand deep crustal features that could not be predicted by flexural theory. Consequently, the best-fit flexural model is used to build a gravity model across the Mérida Andes and the Barinas–Apure Basin preserving the best-controlled structural features from the flexural modelling (e.g. basin wavelength and depth) and slightly changing the main bodies density values and deep crustal structures. The final gravity model is intended to be representative of the major features affecting the gravity field in the study area. The predicted morphology in the lower crustal level of the final gravity model favours the hypothesis of a present delamination or megathrust of the Maracaibo crust over the South American Shield. This process would use the Conrad discontinuity as a main detachment surface within an incipient NW dipping continental subduction.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

Numerical simulations of anomalies of precipitation and surface air temperature in china in the summer of 1997

1.IntroductionStatisticalstudiesdemonstratedthatinEINinoyearstheprecipitationinsummerintheChangjiangRiverandHuaiheRiverBasinsisprobablyabovethenormalwhileitispossiblybelowthenormalinthenorthernChinaandtheHetao(theGreatBendoftheHuangheRiver)region.ThetemperatureinsummerisusuallylowerthannormalinEastAsia,especiallyinNortheastChina.Therewere6yearswithseverelowtemperaturesince1951,andtheyare1954,1957,1964,1972,1976and1983,whichareallrelatedtotheEINinoyears(seeHuangetal.,1989,1992;Xiangand…     

三维定量预测方法在内蒙古红石山南金矿中的应用

红石山南金矿位于内蒙古阿木乌苏—鹰嘴红山金、锑、铁、钨Ⅳ级成矿带中段南部。为了对红石山南金矿的资源潜力进行定量预测评价,并为下一步勘探工作提供指导,在该地区开展以三维可视化技术和信息量法为核心的三维成矿预测很有必要。通过建立红石山金矿的找矿模型,圈定出3个成矿有利区,将二维勘探工作与三维成矿预测有效结合,将矿产预测拓展到三维空间,为今后研究定量化预测评价方法提供参考。     

The role of the seismic trigger in the Calitri landslide (Italy): historical reconstruction and dynamic analysis

The study covered by this paper was focused on the historical case of the Calitri landslide, which was repeatedly reactivated by earthquakes, as reported since 1694. The town of Calitri (Southern Italy) is located on a ridge whose southern slope, from its top to the Ofanto river valley floor, has been historically affected by major landsliding. The last record of recurrence of the Calitri landslide leads back to the 1980 Irpinia earthquake, which caused significant damage to the town and had pervasive and visible ground effects. Based on a detailed historical reconstruction of landsliding and seismicity at Calitri, the study analysed the current static and dynamic stability of the landsliding slope by means of a finite-difference numerical analysis, taking into account the various factors of landslide initiation and kinematics. The results of the numerical analysis: (i) were consistent with the roto-translational mechanism observed upon the latest reactivation of the landslide; (ii) demonstrated that excess pore pressure redistribution caused a lag between the seismic trigger and the initiation of landsliding; and (iii) showed the impact of seismic input frequency on propagation and depth of slope instability.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.