3-D BAROCLINIC NUMERICAL SIMULATION OF THE SOUTH CHINA SEA I. UPPER CIRCULATION

A three-dimensional baroclinic shelf sea model was employed to simulate the seasonal characteristics of the South China Sea (SCS) upper circulation. The results showed that: in summer, an anticyclonic eddy, after its formation between the Bashi Channel and Dongsha Islands in the northeastern SCS, moves southwestward until it disperses slowly. There exists a northward western boundary current along the east shore of the Indo-China Peninsula in the western SCS and an anticyclonic gyre in the southern SCS. But at the end of summer and beginning of autumn, a weak local cyclonic eddy forms in the Nansha Trough, then grows slowly and moves westward till it becomes a cyclonic gyre in the southern SCS in autumn. At the beginning of winter, there exists a cyclonic gyre in the northern and southern SCS, and there is a southward western boundary current along the east shore of the Indo-China Peninsula. But at the end of winter, an anticyclonic eddy grows and moves toward the western boundary after forming in the Nansha Trough. The eddy‘s movement induces a new opposite sign eddy on its eastern side, while the strength of the southward western boundary current gets weakened. This phenomenon continues till spring and causes eddies in the southern SCS.     

Ensemble data assimilation in a simple coupled climate model: The role of ocean-atmosphere interaction

A conceptual coupled ocean-atmosphere model was used to study coupled ensemble data assimilation schemes with a focus on the role of ocean-atmosphere interaction in the assimilation. The optimal scheme was the fully coupled data assimilation scheme that employs the coupled covariance matrix and assimilates observations in both the atmosphere and ocean. The assimilation of synoptic atmospheric variability that captures the temporal fluctuation of the weather noise was found to be critical for the estimation of not only the atmospheric, but also oceanic states. The synoptic atmosphere observation was especially important in the mid-latitude system, where oceanic variability is driven by weather noise. The assimilation of synoptic atmospheric variability in the coupled model improved the atmospheric variability in the analysis and the subsequent forecasts, reducing error in the surface forcing and, in turn, in the ocean state. Atmospheric observation was able to further improve the oceanic state estimation directly through the coupled covariance between the atmosphere and ocean states. Relative to the mid-latitude system, the tropical system was influenced more by ocean-atmosphere interaction and, thus, the assimilation of oceanic observation becomes more important for the estimation of the ocean and atmosphere.     

Hydrometeorological behaviour of pine and larch forests in eastern Siberia

Seasonal changes in the water and energy exchanges over a pine forest in eastern Siberia were investigated and compared with published data from a nearby larch forest. Continuous observations (April to August 2000) were made of the eddy‐correlation sensible heat flux and latent heat flux above the canopy. The energy balance was almost closed, although the sum of the turbulent fluxes sometimes exceeded the available energy flux (R_n ? G) when the latent heat flux was large; this was related to the wind direction. We examined the seasonal variation in energy balance components at this site. The seasonal variation and magnitude of the sensible heat flux (H) was similar to that of the latent heat flux (λE), with maximum values occurring in mid‐June. Consequently, the Bowen ratio was around 1·0 on many days during the study period. On some clear days just after rainfall, λE was very large and the sum of H and λE exceeded R_n ? G. The evapotranspiration rate above the dry canopy from May to August was 2·2 mm day^?1. The contributions of understory evapotranspiration (E_u) and overstory transpiration (E_o) to the evapotranspiration of the entire ecosystem (E_t) were both from 25 to 50% throughout the period analysed. These results suggest that E_u plays a very important role in the water cycle at this site. From snowmelt through the tree growth season (23 April to 19 August 2000), the total incoming water, comprised of the sum of precipitation and the water equivalent of the snow at the beginning of the melt season, was 228 mm. Total evapotranspiration from the forest, including interception loss and evaporation from the soil when the canopy was wet, was 208–254 mm. The difference between the incoming and outgoing amounts in the water balance was from +20 to ?26 mm. The water and energy exchanges of the pine and larch forest differed in that λE and H increased slowly in the pine forest, whereas λE increased rapidly in the larch forest and H decreased sharply after the melting season. Consequently, the shape of the Bowen ratio curves at the two sites differed over the period analysed, as a result of the differences in the species in each forest and in soil thawing. Copyright © 2003 John Wiley & Sons, Ltd.     

Polynomial expansions in likelihoods for spatial data: A case study

This paper illustrates the computational benefits of polynomial representations for quantities in the likelihood function for the spatial linear model based on the power covariance scheme. These benefits include a comprehensive study of likelihoods and maximum likelihood estimators for data. For simplicity, we focus on a relatively simple covariance scheme and data observed at equal intervals along a transect; we briefly indicate how generalizations to more complicated covariance functions and higher dimensions will operate.     

Estimating evapotranspiration for a temperate salt marsh,Newcastle, Australia

Evapotranspiration was studied at a salt marsh site in the Hunter River estuary, NSW, Australia, during 1996–8. Estimates of actual evapotranspiration (E_a) were obtained for three sites using the eddy correlation method. These values were compared with results obtained with the Penman and Penman–Monteith equations, and with pan evaporation. The Penman–Monteith method was found to be most reliable in estimating daily and hourly evapotranspiration. Surface resistance values averaging 12 s m^?1 were derived from the eddy correlation estimates. Recent tidal flooding and rainfall were found to decrease surface resistance and increase E_a/E_p ratios. Estimates of evapotranspiration obtained using the Penman–Monteith method were shown to be sensitive to changes in surface resistance, canopy height and the method used to estimate net radiation from incoming solar radiation. These results underline the importance of accurately estimating such parameters based on site‐specific data rather than relying on empirical equations, which are derived primarily for crops and forests. Copyright © 2001 John Wiley & Sons, Ltd.     

Methodological Development of the Conditional Sampling Method. Part II: Quality Control Criteria of Relaxed Eddy Accumulation Flux Measurements

Determination of biosphere–atmosphere exchanges requires accurate quantification of the turbulent fluxes of energy and of a wide variety of trace gases. Relaxed Eddy Accumulation (REA) is a method that has received increasing attention in recent years, because it does not require any rapid sensor for the scalar measurements as the Eddy Correlation method (EC) does. As in all micrometeorological studies, REA measurements in the atmospheric surface layer are valid under some restrictive conditions so as to be representative of a specific ecosystem. These conditions are the homogeneity of the underlying surface, stationary and horizontally homogeneous turbulence. For most experiments these conditions are not fully satisfied. Data uncertainties can also be related to not fulfilling the method principles or to the technical characteristics of the REA system itself. In order to assess REA measurement quality, a methodological approach of data analysis is developed in this study. This methodological analysis is based on the establishment of criteria for data quality control. A set of them, deduced from the mean and turbulent flow, are called ‘Dynamic criteria’ and are designated to control the stationarity and homogeneity of the w function and the validation of Taylor’s hypothesis. A second set (‘REA operational criteria’) is designed to check the sampling process and, more precisely, the homogeneity of the negative and positive selection process throughout the sampling period. A third set of criteria (‘Chemical scalar criteria’) concerns the scalar measurements. Results of the criteria application to data measured at two different experimental sites are also presented. Cut-off limits of criteria are defined based on their statistical distribution and shown to be specific for each site. Strictness of each criterion, defined by the percentage of flagged samples, is analysed in conjunction with the meteorological conditions and atmospheric stability. It is found that flagged samples mainly correspond to neutral and stable nocturnal conditions. During daytime, nearly free convection conditions can also yield poor quality data.     

On ensemble representation of the observation-error covariance in the Ensemble Kalman Filter

Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003.