Impacts of land use change and climate variations on annual inflow into Miyun Reservoir,Beijing, China

The Miyun Reservoir, the only surface water source for Beijing city, has experienced water supply decline in recent decades. Previous studies suggest that both land use change and climate contribute to the changes of water supply in this critical watershed. However, the specific causes of the decline in the Miyun Reservoir are debatable under a non-stationary climate in the past 4 decades. The central objective of this study was to quantify the separate and collective contributions of land use change and climate variability to the decreasing inflow into the Miyun Reservoir during 1961–2008. Different from previous studies on this watershed, we used a comprehensive approach to quantify the timing of changes in hydrology and associated environmental variables using the long-term historical hydrometeorology and remote-sensing-based land use records. To effectively quantify the different impacts of the climate variation and land use change on streamflow during different sub-periods, an annual water balance model (AWB), the climate elasticity model (CEM), and a rainfall–runoff model (RRM) were employed to conduct attribution analysis synthetically. We found a significant (p  <  0.01) decrease in annual streamflow, a significant positive trend in annual potential evapotranspiration (p  <  0.01), and an insignificant (p  >  0.1) negative trend in annual precipitation during 1961–2008. We identified two streamflow breakpoints, 1983 and 1999, by the sequential Mann–Kendall test and double-mass curve. Climate variability alone did not explain the decrease in inflow to the Miyun Reservoir. Reduction of water yield was closely related to increase in actual evapotranspiration due to the expansion of forestland and reduction in cropland and grassland, and was likely exacerbated by increased water consumption for domestic and industrial uses in the basin. The contribution to the observed streamflow decline from land use change fell from 64–92 % during 1984–1999 to 36–58 % during 2000–2008, whereas the contribution from climate variation climbed from 8–36 % during the 1984–1999 to 42–64 % during 2000–2008. Model uncertainty analysis further demonstrated that climate warming played a dominant role in streamflow reduction in the most recent decade (i.e., 2000s). We conclude that future climate change and variability will further challenge the water supply capacity of the Miyun Reservoir to meet water demand. A comprehensive watershed management strategy needs to consider the climate variations besides vegetation management in the study basin.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.