Bayesian power-law regression with a location parameter,with applications for construction of discharge rating curves

This paper presents a Bayesian approach for fitting the standard power-law rating curve model to a set of stage-discharge measurements. Methods for eliciting both regional and at-site prior information, and issues concerning the determination of prior forms, are discussed. An efficient MCMC algorithm for the specific problem is derived. The appropriateness of the proposed method is demonstrated by applying the model to both simulated and real-life data. However, some problems came to light in the applications, and these are discussed.     

Urban cadastral mapping using very high resolution remote sensing data

With growing urban expanses, one of the pre-requisites for effective governance is Urban Information Systems (UIS) with content down to individual properties (and individuals). The basic input i.e., a map, in UIS should show individual property boundaries showing the plan outline of all structures existing within, at a scale of 1:1000 and larger with sub-metre to centimeters planimetric and geometric accuracy. With very high resolution remote sensing data of the order of 1m available in hand, it is possible to prepare maps with high resolution spatial content. The present exercise demonstrates a method of preparing a geometrically and planimetrically accurate urban cadastral map on very large scale for a small area of about 5 sq km. IKONOS merged data with 1m resolution is used for the purpose. Mapping was done in conjunction with on-site measurements and sketches. Guides are used to maintain shape symmetry and accuracy of buildings and other features. Working out cost of mapping per unit area is another objective in the present exercise. For want of fully or semi-automatic methods of information extraction from very high resolution remote sensing data, it is imperative that mapping should be carried out in conjunction with some on-site measurements wherever necessary.     

Tracking embodied carbon flows in the Belt and Road regions

In the past few decades, economic globalization has driven rapid growth of cross-border trade and a new international division of labor, leading to increasing inter- country embodied carbon flows. Multi-region input-output (MRIO) analysis is used to identify embodied carbon flows between major world regions, including seven regions along the Belt and Road (BR), and the spatial distribution of production- and consumption-based carbon intensities. The results show that current embodied carbon flows are virtually all from BR regions to developed countries, with more than 95% of world net embodied carbon exports coming from BR regions. Consumption in the United States and European Union countries induce about 30% of the carbon emissions in most BR regions, indicating that the former bear a high proportion of consumers’ responsibility for the carbon emitted in the latter. For this reason, measuring environmental responsibilities from consumption rather than a production- based perspective is more equitable, while developing countries should be given a louder voice in the construction through dialogue and cooperation, in part in the context of the Belt and Road Initiative, of an inclusive global climate governance system.     

Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.     

Enrichment characteristics and risk assessment of Hg in bird feathers from Caohai wetland in Guizhou Province,China

Based on the analysis of the enrichment characteristics of Hg and MeHg in bird feathers from Caohai National Nature Reserve in Guizhou, the risks of Hg pollution to the birds from Caohai wetland have been evaluated. The total Hg content of bird feathers ranges from 40 to 5058 ng/g with an average of 924 ng/g. The content of MeHg is significantly correlated with total Hg (r = 0.68, p < 0.01), and the content are among 0.75 and 113 ng/g. The total Hg content in the birds feathers is significantly dependent on their feeding habits, which is mainly in accordance with the following rule: carnivorous birds > omnivorous birds that are mainly carnivorous > omnivorous birds that are mainly herbivorous. There are also differences in the Hg enrichment ability in the different parts of bird feathers, and the total Hg and MeHg content in the wing feathers are significantly higher than that in the other parts of feathers. The bioaccumulation coefficients of aqueous Hg and MeHg by bird feathers are 0.9 × 10⁴–112.13 × 10⁴ (mean value is 20.47 × 10⁴) and 0.47 × 10⁴–70.4 × 10⁴ (mean value is 9.52 × 10⁴), respectively. Although the whole Hg level in Caohai bird feathers is not too high, the Hg content in some carnivorous birds exceeds over or approaches the abnormal threshold when birds are breeding (5 μg/g), which indicates that the birds in Caohai wetland are faced with some risks of ecological Hg pollution.     

Evaluation of CBM reservoir fracture characteristics by azimuthal AVA inversion

Exploitation of coalbed methane (CBM) is directly affected by the strike and density of fractures, particularly because fractures can act as conduits for water and gas (after desorption). To realize the azimuthal amplitude versus angle (AVA) inversion for a CBM reservoir, we treated reservoirs whose fractures are aligned vertically as being equivalent to a horizontal transverse isotropic medium and simplified the P-P wave reflection coefficient formula to an approximately linear equation. Then, we adopted a modified damping least squares algorithm to solve the azimuthal AVA inversion problem. With edited macro-bins, corrected residual moveout, and normalized amplitude, the azimuthal gathers exhibited unambiguous azimuthal AVA characteristics. The #3 CBM reservoir in the Qinshui Basin was selected for a case study. The density and strike of the reservoir fractures were obtained and analyzed, and the reservoir was divided into three classes according to fracture density. Finally, in comparison with the gas production data, we confirmed that the inversion results and reservoir classification were both reasonable and credible.     

Solar-, monsoon- and Kuroshio-influenced thermocline depth and sea surface salinity in the southern Okinawa Trough during the past 17,300 years

Factors influencing millennial-scale variability in the thermocline depth (vertical mixing) and sea surface salinity (SSS) of the southern Okinawa Trough (OT) during the past 17,300 years were investigated based on foraminifer oxygen isotope records of the surface dweller Globigerinoides ruber sensu stricto and the thermocline dweller Pulleniatina obliquiloculata in the AMS ¹⁴C dated OKT-3 core. The thermocline depth is influenced by surface thermal buoyancy (heat) flux, in turn controlled by the annual mean insolation at 30°N and the strength of the East Asian winter monsoon (EAWM). Strong insolation and weak EAWM tend to increase buoyancy gain (decrease buoyancy loss), corresponding to shallow thermocline depths, and vice versa. Regional SSS is influenced by the global ice volume, the Kuroshio Current (KC), and vertical mixing. A deep thermocline coincides with a high SSS because strong vertical mixing brings more, saltier subsurface KC water to the surface, and vice versa. Local SSS (excluding the global ice volume effect) became lower in the northern OT than in the southern OT after ~9.2 ka, implying that Changjiang diluted water had stronger influence in the northern sector. SSS show no major changes during the Bølling/Allerød and Younger Dryas events, probably because the KC disturbed the North Atlantic signals. This argues against earlier interpretations of sea surface temperature records of this core. Wavelet and spectral analyses of the Δδ¹⁸O_P-G (δ¹⁸O of P. obliquiloculata minus G. ruber s.s.) and δ¹⁸O_local records display 1,540-, 1,480-, 1,050-, 860-, 640-, and 630-year periods. These are consistent with published evidence of a pervasive periodicity of 1,500 years in global climate as well as EAWM and KC signatures, and a fundamental solar periodicity of 1,000 years and intermediary derived periodicity of 700 years.     

BVOC emissions from mechanical wounding of leaves and branches of <Emphasis Type="Italic">Eucalyptus sideroxylon</Emphasis> (red ironbark)

Emissions of biogenic volatile organic compounds (BVOCs) from mechanical wounding of leaves and branches of plants can contribute to the atmospheric burden of volatile organic compounds (VOCs) in both (a) urban airsheds (from urban garden maintenance) and (b) the global atmosphere (from large scale forest harvesting). These emissions of BVOCs are poorly understood and quantified, and their role in urban and global emissions inventories neglected. This paper presents measurements of the magnitude, duration and composition of emissions of BVOCs, carbon dioxide (CO₂) and methane (CH₄) from freshly cut leaf mulch and wood chips derived from a common eucalypt tree, Eucalyptus sideroxylon (red ironbark), found in southeastern Australian forests and gardens. The emissions of BVOCs from freshly cut and shredded leaves and wood of E. sideroxylon were found to be 2.3 ± 0.6 and 0.05 ± 0.04 mg g^-1 DM (Dry Mass) from leaf mulch and wood chips respectively and to last typically for 1 day following cutting. Three sampling techniques were used for VOC speciation and the 12 most abundant BVOCs released from the mulch materials were identified. The specific BVOCs emitted in order of decreasing abundance from leaf mulch are: (a) stored plant oils, 1,8-cineole, α–pinene and o-cymene which make up the major part of the emissions, (b) a minor contribution from chemicals associated with environmental stress and wound defence, (Z)–3–hexenyl acetate, (E)-2-hexenal and (Z)-3-hexen-1-ol, and (c) a second minor contribution from metabolic products, acetaldehyde and acetone. The observed integrated emissions of BVOCs from leaves following mulching are equivalent to more than half and perhaps all of the likely stored plant oils in the leaves. For the two comparable studies available, one of a plant with stored oils (this study) and one of a plant without stored plant oils, the emissions of leaf wound defence BVOCs are in the same range for both plants. In the plant with stored plant oils, the plant oil emissions are about a factor of 11 larger in emission rate than the plant wound defence BVOCs. A compilation of available leaf wounding BVOC emission studies indicates that for plants with stored plant oils, plant oil emissions dominate, whereas with other plants, leaf wound defence BVOCs dominate the emissions.     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.