We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling. 相似文献
Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm3). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS 187Os/189Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM 187Os/189Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles. 相似文献
This paper presents a non-destructive, low-cost, photo-based, 3D reconstruction technique for characterizing geo-materials with irregular shapes of a relatively large size. After being validated against two traditional volume measurement methods, namely the vernier caliper method and the fluid displacement method for regular and irregular shapes, respectively, 3D photogrammetry was used to analyse the grout bulbs formed in laboratory pressure grouting tests. The reconstructed 3D mesh model of the sample provides accurate and detailed 3D vertex data, which allowed the volume, densification efficiency and bleeding behaviour of the grout bulbs to be analysed. Comparing the bulb section views at different grouting pressures also offers an intuitive observation of the grout development and propagation process. Moreover, the 3D vertex data and surface area included in the model are of great importance in validating numerical predictions of the pressure grouting process and analysing the interface shear resistance of grouted soil nails or anchors. Compared to existing approaches, the new 3D photogrammetry method possesses several key advantages: (a) it does not require expensive, specialized equipment; (b) samples are not destroyed or modified during testing; (c) it allows to reconstruct objects of various scales and (d) the software is public domain. Therefore, the adoption of this 3D photogrammetry method will facilitate research in the pressure grouting process and can be extended to other problems in geotechnical engineering. 相似文献
Extrusive and intrusive igneous rocks represent different parts of a magmatic system and ultimately provide complementary information about the processes operating beneath volcanoes. To shed light on such processes, we have examined and quantified the textures and mineral compositions of plutonic and cumulate xenoliths and lavas from Bequia, Lesser Antilles arc. Both suites contain assemblages of iddingsitized olivine, plagioclase, clinopyroxene and spinel with rare orthopyroxene and ilmenite. Mineral zoning is widespread, but more protracted in lavas than xenoliths. Plagioclase cores and olivine have high anorthite (An?≤?98) and low forsterite (Fo?≤?84) compositions respectively, implying crystallisation from a hydrous mafic melt that was already fractionated. Xenolith textures range from adcumulate to orthocumulate with variable mineral crystallisation sequences. Textural criteria are used to organize the xenoliths into six groups. Amphibole, notably absent from lavas, is a common feature of xenoliths, together with minor biotite and apatite. Bulk compositions of xenoliths deviate from the liquid line of descent of lavas supporting a cumulate origin with varying degrees of reactive infiltration by evolved hydrous melts, preserved as melt inclusions in xenolith crystals. Volatile saturation pressures in melt inclusions indicate cumulate crystallization over a 162–571 MPa pressure range under conditions of high dissolved water contents (up to 7.8 wt% H2O), consistent with a variety of other thermobarometric estimates. Phase assemblages of xenoliths are consistent with published experimental data on volatile-saturated low-magnesium and high-alumina basalts and basaltic andesite from the Lesser Antilles at pressures of 200–1000 MPa, temperatures of 950–1050 °C and dissolved H2O contents of 4–7 wt%. Once extracted from mid-crustal mushes, residual melts ascend to higher levels and undergo H2O-saturated crystallization in shallow, pre-eruptive reservoirs to form phenocrysts and glomerocrysts. The absence of amphibole from lavas reflects instability at low pressures, whereas its abundance in xenoliths testifies to its importance in mid-crustal differentiation processes. A complex, vertically extensive (6 to at least 21 km depth) magmatic system is inferred beneath Bequia. Xenoliths represent fragments of the mush incorporated into ascending magmas. The widespread occurrence of evolved melts in the mush, but the absence of erupted evolved magmas, in contrast to islands in the northern Lesser Antilles, may reflect the relative immaturity of the Bequia magmatic system. 相似文献
Carbonates in fresh hypabyssal kimberlites worldwide have been studied to understand their origin [i.e. primary magmatic (high T) versus deuteric (‘low T’) versus hydrothermal/alteration (‘low T’)] and identify optimal strategies for petrogenetic studies of kimberlitic carbonates. The approach presented here integrates detailed textural characterisation, cathodoluminescence (CL) imaging, in situ major- and trace-element analysis, as well as in situ Sr-isotope analysis. The results reveal a wide textural diversity. Calcite occurs as fine-grained groundmass, larger laths, segregations, veins or as a late crystallising phase, replacing olivine or early carbonates. Different generations of carbonates commonly coexist in the same kimberlite, each one defined by a characteristic texture, CL response and composition (e.g., variable Sr and Ba concentrations). In situ Sr isotope analysis revealed a magmatic signature for most of the carbonates, based on comparable 87Sr/86Sr values between these carbonates and the coexisting perovskite, a robust magmatic phase. However, this study also shows that in situ Sr isotope analysis not always allow distinction between primary (i.e., magmatic) and texturally secondary carbonates within the same sample. Carbonates with a clear secondary origin (e.g., late-stage veins) occasionally show the same moderately depleted 87Sr/86Sr ratios of primary carbonates and coexisting perovskite (e.g., calcite laths-shaped crystals with 87Sr/86Sr values identical within uncertainty to those of vein calcite in the De Beers kimberlite). This complexity emphasises the necessity of integrating detailed petrography, geochemical and in situ Sr isotopic analyses for an accurate interpretation of carbonate petrogenesis in kimberlites. Therefore, the complex petrogenesis of carbonates demonstrated here not only highlights the compositional variability of kimberlites, but also raises concerns about the use of bulk-carbonate C-O isotope studies to characterise the parental melt compositions. Conversely, our integrated textural and in situ study successfully identifies the most appropriate (i.e. primary) carbonates for providing constraints on the isotopic parameters of parental kimberlite magmas.
On 22 March 2014, a massive, catastrophic landslide occurred near Oso, Washington, USA, sweeping more than 1 km across the adjacent valley flats and killing 43 people. For the following 5 weeks, hundreds of workers engaged in an exhaustive search, rescue, and recovery effort directly in the landslide runout path. These workers could not avoid the risks posed by additional large-scale slope collapses. In an effort to ensure worker safety, multiple agencies cooperated to swiftly deploy a monitoring and alerting system consisting of sensors, automated data processing and web-based display, along with defined communication protocols and clear calls to action for emergency management and search personnel. Guided by the principle that an accelerating landslide poses a greater threat than a steadily moving or stationary mass, the system was designed to detect ground motion and vibration using complementary monitoring techniques. Near real-time information was provided by continuous GPS, seismometers/geophones, and extensometers. This information was augmented by repeat-assessment techniques such as terrestrial and aerial laser scanning and time-lapse photography. Fortunately, no major additional landsliding occurred. However, we did detect small headscarp failures as well as slow movement of the remaining landslide mass with the monitoring system. This was an exceptional response situation and the lessons learned are applicable to other landslide disaster crises. They underscore the need for cogent landslide expertise and ready-to-deploy monitoring equipment, the value of using redundant monitoring techniques with distinct goals, the benefit of clearly defined communication protocols, and the importance of continued research into forecasting landslide behavior to allow timely warning.
Abstract Adiabatic, two-dimensional, steady-state finite-amplitude, hydrostatic gravity waves produced by flow over a ridge are considered. Nonlinear self advection steepens the wave until the streamlines attain a vertical slope at a critical height zc. The height zc, where this occurs, depends on the ridge crest height and adiabatic expansion of the atmosphere. Dissipation is introduced in order to balance nonlinear self advection, and to maintain a marginal state above zc. The approach is to assume that the wave is inviscid except in a thin layer, small compared to a vertical wavelength, where dissipation cannot be neglected. The solutions in each region are matched to obtain a continuous solution for the streamline displacement δ. Solutions are presented for different values of the nondimensional dissipation parameter β. Eddy viscosity coefficients and the thickness of the dissipative layer are expressed as functions of β, and their magnitudes are compared to other theoretical evaluations and to values inferred from radar measurements of the stratosphere. The Fourier spectrum of the solution for z ≫ zc is shown to decay exponentially at large vertical wave numbers n. In comparison, a spectral decay law n?-8/3 characterizes the marginal state of the wave at z = zc. 相似文献
Previous climate risk assessments provide important methodological insights into how to derive tractable research questions and the appropriate use of data under uncertainty, as well as identifying steps that benefit from stakeholder involvement. Here we propose the use of a framework for the systematic and objective exploration of climate risk assessments. The matrix facilitates a breakdown of information about aim and context, main results, methodological choices, stakeholder involvement, sources and characteristics of uncertainties and overall weaknesses. We then apply the matrix to three risk assessments in the water sector to explore some methodological strengths and weaknesses of approaches strongly linked to climate model outputs (top-down) versus those that originate from local knowledge of climate exposures (bottom-up), and demonstrate that closer integration with social and physical sciences is more likely to yield robust climate risk assessments. 相似文献
Abstract The term ‘‘solitary wave'’ is usually used to denote a steadily propagating permanent form solution of a nonlinear wave equation, with the permanency arising from a balance between steepening and dispersive tendencies. It is known that large-scale thermal anomalies in the ocean are subject to a steepening mechanism driven by the beta effect, while at the smaller deformation scale, such phenomena are highly dispersive. It is shown here that the evolution of a physical system subject to both effects is governed by the ‘‘frontal semi-geostrophic equation'’ (FSGE), which is valid for large amplitude thermocline disturbances. Solitary wave solutions of the FSGE (here named planetons) are calculated and their properties are described with a view towards examining the behavior of finite amplitude solitary waves. In contrast, most known solitary wave solutions belong to weakly nonlinear wave equations (e.g., the Korteweg—deVries (KdV) equation). The FSGE is shown to reduce to the KdV equation at small amplitudes. Classical sech2 solitons thus represent a limiting class of solutions to the FSGE. The primary new effect on planetons at finite amplitudes is nonlinear dispersion. It is argued that due to this effect the propagation rates of finite amplitude planetons differ significantly from the ‘‘weak planeton'', or KdV, dispersion relation. Planeton structure is found to be simple and reminiscent of KdV solitons. Numerical evidence is presented which suggests that collisions between finite amplitude solitary waves are weakly inelastic, indicating the loss of true soliton behavior of the FSGE at moderate amplitudes. Lastly, the sensitivity of solitary waves to the existence of a nontrivial far field is demonstrated and the role of this analysis in the interpretation of lab experiments and the evolution of the thermocline is discussed. 相似文献
Here we apply quantitative technique to describe the seafloor seepages based on the multi-beam backscatter and bathymetric investigations to characterize the pockmark morphology. The variable seafloor backscatter strength for coarser seafloor sediments are related to the diagenesis derived from biodegraded seepages. In this regard, box counting method is used to estimate ‘fractal dimension’ for backscatter imagery data of 398 blocks. These blocks are further sub-grouped into six classes depending on the spread of pockmark related seepages. The study area lies 102 km west off Marmagao along the central west coast of India which contains pre-dominantly (70%) gas-charged sediments. Comparison between the estimated self-similar fractals reveals that there is approximately 97% correlation between the box (Dbox) and information (Dinfo) dimensions. Box dimension–derived fractal dimension values, suggest that the seepages are more along the fault trace in deeper waters, in comparison to sparsely distributed shallow water seepages. Besides, this poor seepage is confined within the smooth to moderately rough seafloor. It is established that the high backscatter strength along the upper slope of the pockmark region having higher fractal dimensions reflects multifractal behavior of seepage distribution. Entire area indicates patchy seepage patterns as supported by estimated fractal values showing intermittent fluctuations, which emphasizes non-linear behavior. Estimated self organizing criticality (SOC) parameters for six representative blocks reveal that the nature of pockmark, fault trace, sediment nature coupled with slumping of pockmark’s wall, sediment movement due to bottom currents are controlling the dynamic balance in the area seepage system. Further, our study emphasizing the multifractal behavior of seepage blocks, clearly depicts the drift in the seepage pattern. 相似文献