Analysis of long‐term water balance in the source area of the Yellow River basin

To analyse the long‐term water balance of the Yellow River basin, a new hydrological model was developed and applied to the source area of the basin. The analysis involved 41 years (1960–2000) of daily observation data from 16 meteorological stations. The model is composed of the following three sub‐models: a heat balance model, a runoff formation model and a river‐routing network model. To understand the heat and water balances more precisely, the original model was modified as follows. First, the land surface was classified into five types (bare, grassland, forest, irrigation area and water surface) using a high‐resolution land‐use map. Potential evaporation was then calculated using land‐surface temperatures estimated by the heat balance model. The maximum evapotranspiration of each land surface was calculated from potential evaporation using functions of the leaf area index (LAI). Finally, actual evapotranspiration was estimated by regulating the maximum evapotranspiration using functions of soil moisture content. The river discharge estimated by the model agreed well with the observed data in most years. However, relatively large errors, which may have been caused by the overestimation of surface flow, appeared in some summer periods. The rapid decrease of river discharge in recent years in the source area of the Yellow River basin depended primarily on the decrease in precipitation. Furthermore, the results suggested that the long‐term water balance in the source area of the Yellow River basin is influenced by land‐use changes. Copyright © 2007 John Wiley & Sons, Ltd.     

Observed and simulated sensitivities of summertime urban surface air temperatures to anthropogenic heat in downtown areas of two Japanese Major Cities,Tokyo and Osaka

In this study, the sensitivities of surface air temperatures to anthropogenic heat (AH) were investigated in downtowns of the two Japanese major cities, Tokyo and Osaka. First, meteorological measurements were made with the simultaneous monitoring of electricity demand in a contrastive couple of a downtown commercial area (C-area) and a residential area (R-area) within each city in summer 2007. From the measurements, the areal-mean surface air temperatures were obtained as \( {\overline{T}}_{\mathrm{C}} \) and \( {\overline{T}}_{\mathrm{R}} \) for each of the C-area and R-area, respectively. Using the actual electricity demand and the estimated motor fuels consumption, their areal total was evaluated as the energy-consumption-basis AH. The estimated C-areas' AH indicated greater values up to 220 W/m² on weekdays and remarkable decrease about by half on weekends, whereas that in the R-areas showed less values of 10–20 W/m² stably. Then, \( {\overline{T}}_{\mathrm{C}}-{\overline{T}}_{\mathrm{R}} \) on calm and fine days were found to be systematically decreased from weekdays to weekends in both cities roughly indicating a proportional relationship with the reductions in the C-areas' AH on weekends. The result suggested a common afternoon sensitivity for both C-areas of around 1.0°C/100 W/m², which indicated an intensity of the AH impact on surface air temperature there. Next, to simulate the observed AH impact, the authors' CM-BEM (a multilayer urban canopy model coupled with a building energy model) was newly implemented in the mesoscale Weather Research and Forecasting (WMF) model. This new system, WRF-CM-BEM, was applied to Tokyo and almost reasonably validated from the aspects of the reproducibility of urban surface air temperature and electricity demand in the observation areas. The simulations also suggested that WRF-CM-BEM underestimated the observed air temperature sensitivity to AH in the Tokyo C-area roughly by half but still in the same order of magnitude.     

Numerical computation of incomplete elliptic integrals of a general form

We present an algorithm to compute the incomplete elliptic integral of a general form. The algorithm efficiently evaluates some linear combinations of incomplete elliptic integrals of all kinds to a high precision. Some numerical examples are given as illustrations. This enables us to numerically calculate the values and the partial derivatives of incomplete elliptic integrals of all kinds, which are essential when dealing with many problems in celestial mechanics, including the analytic solution of the torque-free rotational motion of a rigid body around its barycenter.     

Seasonal variation in the energy and water exchanges above and below a larch forest in eastern Siberia

The water and energy exchanges in forests form one of the most important hydro‐meteorological systems. There have been far fewer investigations of the water and heat exchange in high latitude forests than of those in warm, humid regions. There have been few observations of this system in Siberia for an entire growing season, including the snowmelt and leaf‐fall seasons. In this study, the characteristics of the energy and water budgets in an eastern Siberian larch forest were investigated from the snowmelt season to the leaf‐fall season. The latent heat flux was strongly affected by the transpiration activity of the larch trees and increased quickly as the larch stand began to foliate. The sensible heat dropped at that time, although the net all‐wave radiation increased. Consequently, the seasonal variation in the Bowen ratio was clearly ‘U’‐shaped, and the minimum value (1·0) occurred in June and July. The Bowen ratio was very high (10–25) in early spring, just before leaf opening. The canopy resistance for a big leaf model far exceeded the aerodynamic resistance and fluctuated over a much wider range. The canopy resistance was strongly restricted by the saturation deficit, and its minimum value was 100 s m^?1 (10 mm s^?1 in conductance). This minimum canopy resistance is higher than values obtained for forests in warm, humid regions, but is similar to those measured in other boreal conifer forests. It has been suggested that the senescence of leaves also affects the canopy resistance, which was higher in the leaf‐fall season than in the foliated season. The mean evapotranspiration rate from 21 April 1998 to 7 September 1998 was 1·16 mm day^?1, and the maximum rate, 2·9 mm day^?1, occurred at the beginning of July. For the growing season from 1 June to 31 August, this rate was 1·5 mm day^?1. The total evapotranspiration from the forest (151 mm) exceeded the amount of precipitation (106 mm) and was equal to 73% of the total water input (211 mm), including the snow water equivalent. The understory evapotranspiration reached 35% of the total evapotranspiration, and the interception evaporation was 15% of the gross precipitation. The understory evapotranspiration was high and the interception evaporation was low because the canopy was sparse and the leaf area index was low. Copyright © 2001 John Wiley & Sons, Ltd.     

Elements of spin motion

For use in numerical studies of rotational motion, a set of elements is introduced for the torque-free rotational motion of a rigid body around its barycenter. The elements are defined as the initial values of a modification of the Andoyer canonical variables. A computational procedure is obtained for determining these elements from the combination of the spin angular momentum vector and a triad defining the orientation of the rigid body. A numerical experiment shows that the errors of transformation between the elements and variables are sufficiently small. The errors increase linearly with time for some elements and quadratically for some others.     

Semi-empirical estimation of ground motion using observed records at a site in Shikoku, Japan

A semi-empirical approach using fore- or after-shockrecords as Green's functions is applicable to thesimulation of strong ground motion, however suchrecords are obviously not available for predictionpurposes. Thus we have predicted ground motion fora hypothetical large earthquake from other minorevents by adopting a distance correction based ongeometrical spreading. Another difficulty inprediction is fault modeling. Surface traces weresimplified as fault models 27, 46, 55, and 77 km inlength. Further, the actual fault rupture may beinhomogeneous, so an asperity distribution isassumed. This asperity model assumes thatdislocation and stress drop are double than theaverage values. Although, the near field term isneglected in our simulation, no significantdifference was seen in the motions estimated byindividual models for periods up to 2.0 seconds. This indicates that the dependence of source size issmall for strong motion, perhaps as a result of therandom summation of high-frequency phases.     

Rectangular rotation of spherical harmonic expansion of arbitrary high degree and order

In order to move the polar singularity of arbitrary spherical harmonic expansion to a point on the equator, we rotate the expansion around the y-axis by \(90^{\circ }\) such that the x-axis becomes a new pole. The expansion coefficients are transformed by multiplying a special value of Wigner D-matrix and a normalization factor. The transformation matrix is unchanged whether the coefficients are \(4 \pi \) fully normalized or Schmidt quasi-normalized. The matrix is recursively computed by the so-called X-number formulation (Fukushima in J Geodesy 86: 271–285, 2012a). As an example, we obtained \(2190\times 2190\) coefficients of the rectangular rotated spherical harmonic expansion of EGM2008. A proper combination of the original and the rotated expansions will be useful in (i) integrating the polar orbits of artificial satellites precisely and (ii) synthesizing/analyzing the gravitational/geomagnetic potentials and their derivatives accurately in the high latitude regions including the arctic and antarctic area.     

Rapid and Simple Measurement of H2 Emission from Active Faults Using Compact Sampling Equipments

Smaller equipment has been developed for sampling and measuring H₂ in fault zones. A considerable volume of hydrogen emission, suggesting the presence of fluid paths along fault zones, was detected within a 1‐h simple field procedure at an active fault in Central Japan, the Atotsugawa Fault. The equipment enabled measurements in a thin fault gouge with a thickness of 10 mm; this led to pervasive sampling from small fault zones hosted in harder rocks. A rapid evaluation of the spatio‐temporal heterogeneity of hydrogen emissions along the faults, using the present method, would increase knowledge of fluid circulation around faults.     

Dissolved organic carbon in a eutrophic lake; dynamics,biodegradability and origin

The seasonal and spatial changes in dissolved organic carbon (DOC) in Lake Kasumigaura, a shallow, eutrophic lake, were analyzed and the lability of DOC was tested by long-term incubations. There was a nearly 1 mgCl^–1 downstream increase in refractory DOC in the lake; at the center it fluctuated little seasonally. The characteristic UV-absorbance: DOC ratios were determined for samples from the influent rivers (pedogenic: used interchangeably with allochthonous) and outdoor experimental ponds (autochthonous) during incubations. These ratios were then used to calculate the proportion of total measured lake water DOC in each of four components: pedogenic-refractory (PR), pedogenic-labile (PL), autochthonous-refractory (AR) and autochthonous-labile (AL). PR was uniform (around 1.5 mgCl^–1) or diminished very slightly over time. AR increased from nearly zero at the station closest to an influent river to 1 mgCl^–1 at the lake center. PL declined downstream from 0.3 mgCl^–1 to zero. AL was virtually constant at 0.8 mgCl^–1 except at the station closest to the influent river. The constancy of the UV-absorbance: DOC ratio during the biodegradation process was confirmed for Lake Kasumigaura; hence a two-component model (pedogenic-autochthonous) could be applied here without consideration of DOC lability. However, this assumption is not always met for other water bodies, and therefore it should be checked before applying a two-component model elsewhere.     

In-water algorithms for ADEOS/OCTS

In-water algorithms for OCTS standard products were developed using in situ data and installed for operationally processing at NASDA/EOC. This paper describes the in-water algorithms Version 1.0 for chlorophylla concentration, pigment concentration, and attenuation coefficient at a wavelength of 490 nm. The selected OCTS standard algorithms (Ver. 1.0) are as follows: