A numerical method for estimation of impulse response and its application to laboratory data     

Masayuki Kikuchi 《Physics of the Earth and Planetary Interiors》1979,19(3):255-262

Following Burg's algorithm, we develop a new mathematical method of estimating the impulse response of a linear system. Since no data-invariant time window is used in this method, the estimated impulse response gives us reliable information on both the amplitude and the phase spectra of the system.The method is applied to laboratory data to investigate the dynamical response of rock specimens. The experimental results show that the waveform distortion with travel distance may well be described by a standard viscoelastic solid model. However, they also show that there is an additional factor of attenuation which is independent of the frequency and therefore not accompanied by any velocity dispersion.  相似文献   

Statistical analysis of polar motion data — I. Theory for the Chandler wobble analysis     

Isao Okamoto  Naokichi Kikuchi 《Geophysical Journal International》1982,69(3):669-719

  相似文献   

Retardation capacity of altered granitic rock distributed along fractured and faulted zones in the orogenic belt of Japan     

H. Yoshida  R. Metcalfe  Y. Seida  H. Takahashi  T. Kikuchi 《Engineering Geology》2009,106(3-4):116-122

There are wide areas of granitic rocks in the Japanese orogenic belt. These granitic bodies inevitably contain fracture and fault systems associated with alteration zones. However, relatively little attention has been given to the possible influence of such widely distributed alteration zones on the migration of radionuclides from any radioactive waste repository that might in future be sited within granitic rock. In particular, the influences of alteration products and micro-fractures, due to chemical sorption and/or physical retardation require further consideration. In order to understand the retardation capacity of the altered deep granitic rocks, detailed geometrical characterization of pores, geochemical analysis, and batch sorption and flow-through experiments have been carried out. Those results show that the altered granite has a large volume of accessible pores, particularly in potassium-feldspar grains, which would influence nuclide retardation more than the accessible porosity in other minerals present, such as biotite. The distribution coefficients, Kd estimated from batch sorption tests and flow-through experiments suggest that altered granite has a high capability to retard the migration of nuclides. The retardation would probably be due to sorption on altered minerals such as sericite and iron hydroxides formed along grain boundaries and in pores created by dissolution, in addition to sorption on primary sorptive minerals. These results provide confidence that even altered and fractured parts of any granitic rock that might be encountered in a site for the disposal of high level radioactive waste may still retard radionuclide migration and thereby help the geosphere to function as a barrier.  相似文献   

Toward Increased Use of Data Worth Analyses in Groundwater Studies          下载免费PDF全文

Colin Kikuchi 《Ground water》2017,55(5):670-673

  相似文献   

Analysis of subsurface temperature data to quantify groundwater recharge rates in a closed Altiplano basin,northern Chile   总被引：1，自引：0，他引：1  

C. P. Kikuchi  T. P. A. Ferré 《Hydrogeology Journal》2017,25(1):103-121

Quantifying groundwater recharge is a fundamental part of groundwater resource assessment and management, and is requisite to determining the safe yield of an aquifer. Natural groundwater recharge in arid and semi-arid regions comprises several mechanisms: in-place, mountain-front, and mountain-block recharge. A field study was undertaken in a high-plain basin in the Altiplano region of northern Chile to quantify the magnitude of in-place and mountain-front recharge. Water fluxes corresponding to both recharge mechanisms were calculated using heat as a natural tracer. To quantify in-place recharge, time-series temperature data in cased boreholes were collected, and the annual fluctuation at multiple depths analyzed to infer the water flux through the unsaturated zone. To quantify mountain-front recharge, time-series temperature data were collected in perennial and ephemeral stream channels. Streambed thermographs were analyzed to determine the onset and duration of flow in ephemeral channels, and the vertical water fluxes into both perennial and ephemeral channels. The point flux estimates in streambeds and the unsaturated zone were upscaled to channel and basin-floor areas to provide comparative estimates of the range of volumetric recharge rates corresponding to each recharge mechanism. The results of this study show that mountain-front recharge is substantially more important than in-place recharge in this basin. The results further demonstrate the worth of time-series subsurface temperature data to characterize both in-place and mountain-front recharge processes.  相似文献   

Lunar photometric properties at wavelengths 0.5-1.6 μm acquired by SELENE Spectral Profiler and their dependency on local albedo and latitudinal zones     

Yasuhiro Yokota  Tsuneo Matsunaga  Junichi Haruyama  Satoru Yamamoto  Tomokatsu Morota  Kazuto Saiki  Kohei Kitazato  Akira Iwasaki  Naru Hirata  Rie Honda  Hitoshi Mizutani 《Icarus》2011,215(2):639-660

The lunar photometric function, which describes the dependency of the observed radiance on the observation geometry, is used for photometric correction of lunar visible/near-infrared data. A precise photometric correction parameter set is crucial for many applications including mineral identification and reflectance map mosaics. We present, for the first time, spectrally continuous photometric correction parameters for both sides of the Moon for wavelengths in the range 0.5-1.6 μm and solar phase angles between 5° and 85°, derived from Kaguya (SELENE) Spectral Profiler (SP) data. Since the measured radiance also depends on the surface albedo, we developed a statistical method for selecting areas with relatively uniform albedos from a nearly 7000-orbit SP data set. Using the selected data set, we obtained empirical photometric correction parameter sets for three albedo groups (high, medium, and low). We did this because the photometric function depends on the albedo, especially at phase angles below about 20° for which the shadow hiding opposition effect is appreciable. We determined the parameters in 160 bands and discovered a small variation in the opposition effect due to the albedo variation of mafic mineral absorption. The consistency of the photometric correction was checked by comparing observations made at different times of the same area on the lunar surface. Variations in the spectra obtained were lower than 2%, except for the large phase angle data in mare. Lastly, we developed a correction method for low solar elevation data, which is required for high latitude regions. By investigating low solar elevation data, we introduced an additional correction method. We used the new photometric correction to generate a 1° mesh global lunar reflectance map cube in a wavelength range of 0.5-1.6 μm. Surprisingly, these maps reveal that high latitude (?75°) regions in both the north and south have much lower spectral continuum slopes (color ratio r_1547.7nm/r_752.8nm ? 1.8) than the low and medium latitude regions, which implies lower degrees of space weathering.  相似文献   

Global distributions of storm-time ionospheric currents as seen in geomagnetic field variations     

Atsuki Shinbori  Tomoaki Hori  Yoshimasa Tanaka  Yukinobu Koyama  Takashi Kikuchi  Tsutomu Nagatsuma 《极地研究(英文版)》2013,(4):296-314

To investigate temporal and spatial evolution of global geomagnetic field variations from high-latitude to the equator during geomagnetic storms, we analyzed ground geomagnetic field disturbances from high latitudes to the magnetic equator. The daytime ionospheric equivalent current during the storm main phase showed that twin-vortex ionospheric currents driven by the Region 1 field-aligned currents （R1 FACs） are intensified significantly and expand to the low-latitude region of-30~ magnetic latitude. Centers of the currents were located around 70~ and 65~ in the morning and afternoon, respectively. Corresponding to intensification of the R1 FACs, an enhancement of the eastward/westward equatorial electrojet occurred at the daytime/nighttime dip equator. This signature suggests that the enhanced convection electric field penetrates to both the daytime and nighttime equa- tor. During the recovery phase, the daytime equivalent current showed that two new pairs of twin vortices, which are different from two-cell ionospheric currents driven by the R1 FACs, appear in the polar cap and mid latitude. The former led to enhanced north- ward Bz （NBZ） FACs driven by lobe reconnection tailward of the cusps, owing to the northward interplanetary magnetic field （IMF）. The latter was generated by enhanced Region 2 field-aligned currents （R2 FACs）. Associated with these magnetic field variations in the mid-latitudes and polar cap, the equatorial magnetic field variation showed a strongly negative signature, produced by the westward equatorial electrojet current caused by the dusk-to-dawn electric field.  相似文献   

Configuring high frequency radar observations in the Southern Chukchi Sea   总被引：1，自引：1，他引：0  

Gleb Panteleev  Max Yaremchuk  Oceana Francis  Takashi Kikuchi 《Polar Science》2013,7(2):72-81

In recent years, monitoring offshore surface circulation in the Arctic Ocean with high frequency radars has become an issue of increasing practical importance. In this study, radar positions are optimized by minimizing the reconstruction errors of the surface currents in the Southeastern Chukchi Sea. By means of an adjoint sensitivity technique it is shown that in the case of a pair of radars, their optimal (i.e. most favorable) location is at Kivalina, a settlement near the strongest outflow of the Alaskan Coastal Current from the monitored domain. The least favorable location is at Shishmaref, a settlement near relatively weak inflow into the region as observed from the coast. However, if two pairs of radars are available, the best locations are Kivalina and Shishmaref. The results are verified using observational system simulation experiments (OSSEs) performed in the framework of a 4-dimensional variational assimilation of simulated radar observations into a numerical model. It is shown that correct specification of the first guess solution is of primary importance for obtaining realistic results from both adjoint sensitivity analysis and OSSEs. This emphasizes the necessity of obtaining accurate high resolution climatologies for future ice-free offshore regions in the Arctic.  相似文献   

Distribution of local Cs anomalies on the seafloor near the Fukushima Dai-ichi Nuclear Power Plant     

Blair Thornton  Seiki Ohnishi  Tamaki Ura  Naoteru Odano  Shun Sasaki  Tsuneo Fujita  Tomowo Watanabe  Kaoru Nakata  Tsuneo Ono  Daisuke Ambe 《Marine pollution bulletin》2013

An estimated 3.5 ± 0.7 × 10¹⁵ Bq of ¹³⁷Cs is thought to have been discharged into the ocean following the melt down at Fukushima Dai-ichi Nuclear Power Plant (F1NPP). While efforts have been made to monitor seafloor radiation levels, the sampling techniques used cannot capture the continuous distribution of radionuclides. In this work, we apply in situ measurement techniques using a towed gamma ray spectrometer to map the continuous distribution of ¹³⁷Cs on the seafloor within 20 km of the F1NPP. The results reveal the existence of local ¹³⁷Cs anomalies, with levels of ¹³⁷Cs an order of magnitude higher than the surrounding seafloors. The sizes of the anomalies mapped in this work range from a few meters to a few hundreds of meters in length, and it is demonstrated that the distribution of these anomalies is strongly influenced by meter scale features of the terrain.  相似文献   

Lunar gravity field determination using SELENE same-beam differential VLBI tracking data   总被引：2，自引：0，他引：2  

S. Goossens  K. Matsumoto  Q. Liu  F. Kikuchi  K. Sato  H. Hanada  Y. Ishihara  H. Noda  N. Kawano  N. Namiki  T. Iwata  F. G. Lemoine  D. D. Rowlands  Y. Harada  M. Chen 《Journal of Geodesy》2011,85(4):205-228

A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from 1-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080±0.0009 km³/s² (10 sigma). The lunar degree 2 potential Love number k ₂ was also estimated, and has a value of 0.0255 ± 0.0016 (10 sigma as well).  相似文献