Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Simulations of multimode bolometric interferometers

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.     

Merging of low-mass systems and the origin of the fundamental plane

We present the preliminary results of a study of how small stellar systems merge to form larger ones. As we display the families of galaxies in the μ_e - R_e plane (effective surface brightness versus effective radius) we realize that different morphological types occupy different loci, evidencing the different physical mechanisms operating in each family. As proposed by Capaccioli et al. (1992) this diagram is the logical equivalent of the HR diagram for stars. Here we take some initial steps in understanding of how we can establish the evolutionary tracks, solely due to dynamical processes, in the μ_e - R_e plane, ultimately making a dwarf elliptical to turn into a normal elliptical galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Seismic Inference using Genetic Algorithms

A flood of reliable seismic data will soon arrive. The migration to largertelescopes on the ground may free up 4-m class instruments for multi-sitecampaigns, and several forthcoming satellite missions promise to yieldnearly uninterrupted long-term coverage of many pulsating stars. We willthen face the challenge of determining the fundamental properties of thesestars from the data, by trying to match them with the output of ourcomputer models. The traditional approach to this task is to make informedguesses for each of the model parameters, and then adjust them iterativelyuntil an adequate match is found. The trouble is: how do we know that oursolution is unique, or that some other combination of parameters will notdo even better? Computers are now sufficiently powerful and inexpensivethat we can produce large grids of models and simply compare all ofthem to the observations. The question then becomes: what range ofparameters do we want to consider, and how many models do we want tocalculate? This can minimize the subjective nature of the process, but itmay not be the most efficient approach and it may give us a false sense ofsecurity that the final result is correct, when it is really justoptimal. I discuss these issues in the context of recent advances inthe asteroseismological analysis of white dwarf stars.