Carbon dioxide and methane exchange between a boreal pristine lake and the atmosphere

Lakes serve as sites for terrestrially fixed carbon to be remineralized and transferred back to the atmosphere. Their role in regional carbon cycling is especially important in the Boreal Zone, where lakes can cover up to 20% of the land area. Boreal lakes are often characterized by the presence of a brown water colour, which implies high levels of dissolved organic carbon from the surrounding terrestrial ecosystem, but the load of inorganic carbon from the catchment is largely unknown. Organic carbon is transformed to methane (CH4) and carbon dioxide (CO2) in biological processes that result in lake water gas concentrations that increase above atmospheric equilibrium, thus making boreal lakes as sources of these important greenhouse gases. However, flux estimates are often based on sporadic sampling and modelling and actual flux measurements are scarce. Thus, the detailed temporal flux dynamics of greenhouse gases are still largely unknown. ----- One aim here was to reveal the natural dynamics of CH4 and CO2 concentrations and fluxes in a small boreal lake. The other aim was to test the applicability of a measuring technique for CO2 flux, i.e. the eddy covariance (EC) technique, and a computational method for estimation of primary production and community respiration, both commonly used in terrestrial research, in this lake. Continuous surface water CO2 concentration measurements, also needed in free-water applications to estimate primary production and community respiration, were used over two open water periods in a study of CO2 concentration dynamics. Traditional methods were also used to measure gas concentration and fluxes. The study lake, Valkea-Kotinen, is a small, humic, headwater lake within an old-growth forest catchment with no local anthropogenic disturbance and thus possible changes in gas dynamics reflect the natural variability in lake ecosystems. CH4 accumulated under the ice and in the hypolimnion during summer stratification. The surface water CH4 concentration was always above atmospheric equilibrium and thus the lake was a continuous source of CH4 to the atmosphere. However, the annual CH4 fluxes were small, i.e. 0.11 mol m-2 yr-1, and the timing of fluxes differed from that of other published estimates. The highest fluxes are usually measured in spring after ice melt but in Lake Valkea-Kotinen CH4 was effectively oxidised in spring and highest effluxes occurred in autumn after summer stratification period. CO2 also accumulated under the ice and the hypolimnetic CO2 concentration increased steadily during stratification period. The surface water CO2 concentration was highest in spring and in autumn, whereas during the stable stratification it was sometimes under atmospheric equilibrium. It showed diel, daily and seasonal variation; the diel cycle was clearly driven by light and thus reflected the metabolism of the lacustrine ecosystem. However, the diel cycle was sometimes blurred by injection of hypolimnetic water rich in CO2 and the surface water CO2 concentration was thus controlled by stratification dynamics. The highest CO2 fluxes were measured in spring, autumn and during those hypolimnetic injections causing bursts of CO2 comparable with the spring and autumn fluxes. The annual fluxes averaged 77 (±11 SD) g C m-2 yr-1. In estimating the importance of the lake in recycling terrestrial carbon, the flux was normalized to the catchment area and this normalized flux was compared with net ecosystem production estimates of -50 to 200 g C m-2 yr-1 from unmanaged forests in corresponding temperature and precipitation regimes in the literature. Within this range the flux of Lake Valkea-Kotinen yielded from the increase in source of the surrounding forest by 20% to decrease in sink by 5%. The free water approach gave primary production and community respiration estimates of 5- and 16-fold, respectively, compared with traditional bottle incubations during a 5-day testing period in autumn. The results are in parallel with findings in the literature. Both methods adopted from the terrestrial community also proved useful in lake studies. A large percentage of the EC data was rejected, due to the unfulfilled prerequisites of the method. However, the amount of data accepted remained large compared with what would be feasible with traditional methods. Use of the EC method revealed underestimation of the widely used gas exchange model and suggests simultaneous measurements of actual turbulence at the water surface with comparison of the different gas flux methods to revise the parameterization of the gas transfer velocity used in the models.

Convergence feedback and unstable low frequency oscillations in a simple coupled ocean‐atmosphere model

The role of convergence feedback on the stability of a coupled ocean‐atmosphere system is studied using model III of Hirst (1986). It is shown that the unstable coupled mode found by Hirst is greatly modified by the convergence feedback. If the convergence feedback strength exceeds a critical value, several new unstable intraseasonal modes are also introduced. These modes have very weak dependence on the wave number. These results may explain the behaviour of some coupled models and to some extent provide a mechanism for the observed aperiodicity of the El‐Nino and Southern Oscillation (ENSO) events.

Propagation of quasi-simple waves in a compressible rotating atmosphere

A class of self-propagating linear and nonlinear travelling wave solutions for compressible rotating fluid is studied using both numerical and analytical techiques. It is shown that, in general, a three dimensional linear wave is not periodic. However, for some range of wave numbers depending on rotation, horizontally propagating waves are periodic. When the rotation ohgr is equal to $$\sqrt {(\gamma - 1)/(4\gamma )}$$ , all horizontal waves are periodic. Here, gamma is the ratio of specific heats. The analytical study is based on phase space analysis. It reveals that the quasi-simple waves are periodic only in some plane, even when the propagation is horizontal, in contrast to the case of non-rotating flows for which there is a single parameter family of periodic solutions provided the waves propagate horizontally. A classification of the singular points of the governing differential equations for quasi-simple waves is also appended.

Long waves in inviscid compressible atmospheres: Isothermal atmosphere

Solitary waves and cnoidal waves have been found in an adiabatic compressible atmosphere which, under ambient conditions, has winds, and is isothermal. The theory is illustrated with an example for which the background wind is linearly increasing. It is found that the number of possible critical speeds of the flow depends crucially on whether the Richardson number is greater or less than one‐fourth.

Long waves in inviscid compressible atmospheres: Isothermal atmosphere

Solitary waves and cnoidal waves have been found in an adiabatic compressible atmosphere which, under ambient conditions, has winds, and is isothermal. The theory is illustrated with an example for which the background wind is linearly increasing. It is found that the number of possible critical speeds of the flow depends crucially on whether the Richardson number is greater or less than one‐fourth.

On the Existence of Hydromagnetic Interface Waves in a Structured Atmosphere

The conditions under which the hydromagnetic interface waves can exist at a magnetic interface is deduced. Using these conditions, it is shown that a slow interface wave with a phase velocity about 5Km/s and a fast interface wave with a phase velocity 6.5 to 8km/s at the photospheric level can exist.

Mechanisms of variability and predictability of the tropical coupled ocean-atmosphere system

A conceptual model is proposed to explain the observed aperiodicity in the short term climate fluctuations of the tropical coupled ocean-atmosphere system. This is based on the evidence presented here that the tropical coupled ocean-atmosphere system sustains a low frequency inter-annual mode and a host of higher frequency intra-seasonal unstable modes. At long wavelengths, the low frequency mode is dominant while at short wavelengths, the high frequency modes are dominant resulting in the co-existence of a long wave low frequency mode with some short wave intra-seasonal modes in the tropical coupled system. It is argued that due to its long wavelength, the low frequency mode would behave like a linear oscillator while the higher frequency short wave modes would be nonlinear. The conceptual model envisages that an interaction between the low frequency linear oscillator and the high frequency nonlinear oscillations results in the observed aperiodicity of the tropical coupled system. This is illustrated by representing the higher frequency intra-seasonal oscillations by a nonlinear low order model which is then coupled to a linear oscillator with a periodicity of four years. The physical mechanism resulting in the aperiodicity in the low frequency oscillations and implications of these results on the predictability of the coupled system are discussed.

Prediction of the Indian summer monsoon rainfall using a state-of-the-art coupled ocean-atmosphere model

A state-of-the-art model of the coupled ocean-atmosphere system, the climate forecast system (CFS), from the National Centres for Environmental Prediction (NCEP), USA, has been ported onto the PARAM Padma parallel computing system at the Centre for Development of Advanced Computing (CDAC), Bangalore and retrospective predictions for the summer monsoon (June-September) season of 2009 have been generated, using five initial conditions for the atmosphere and one initial condition for the ocean for May 2009. Whereas a large deficit in the Indian summer monsoon rainfall (ISMR; June-September) was experienced over the Indian region (with the all-India rainfall deficit by 22% of the average), the ensemble average prediction was for above-average rainfall during the summer monsoon. The retrospective predictions of ISMR with CFS from NCEP for 1981-2008 have been analysed. The retrospective predictions from NCEP for the summer monsoon of 1994 and that from CDAC for 2009 have been compared with the simulations for each of the seasons with the stand-alone atmospheric component of the model, the global forecast system (GFS), and observations. It has been shown that the simulation with GFS for 2009 showed deficit rainfall as observed. The large error in the prediction for the monsoon of 2009 can be attributed to a positive Indian Ocean Dipole event seen in the prediction from July onwards, which was not present in the observations. This suggests that the error could be reduced with improvement of the ocean model over the equatorial Indian Ocean.

Diurnal variation of atmospheric CO2 concentration and delta C-13 in an urban atmosphere during winter-role of the Nocturnal Boundary Layer

The paper presents the importance of the Nocturnal Boundary Layer in driving the diurnal variability of the atmospheric CO2 mixing ratio and the carbon isotope ratio at ground level from an urban station in India. Our observations are the first of their kind from this region. The atmospheric CO2 mixing ratio and the carbon isotopic ratio were measured for both the morning (05:30-07:30 IST) and afternoon time (16:00-18:00 IST) air samples at 5 m above ground level in Bangalore city, Karnataka State (12 degrees 58' N, 77 degrees 38' E, masl = 920 m) for a 10 day period during the winter of 2008. We observed a change of similar to 7% the in CO2 mixing ratio between the morning and afternoon time air samples. A stable isotope analysis of CO2 from morning samples showed a depletion in the carbon isotope ratio by similar to 2 parts per thousand compared to the afternoon samples. Along with the ground-based measurement of air samples, data of radiosonde measurements were also obtained from the Indian Meteorological Department to identify the vertical atmospheric structure at different time in a day. We proposed the presence or absence of the NBL as a controlling factor for the observed variability in the mixing ratio as well as its isotopic composition. Here we used the Keeling model approach to find out the carbon isotope ratio for the local sources. The local sources have further been characterized as anthropogenic and biological respiration (in %) using a two-component mixing model. We also used a vertical mixing model based on the concept of the mixing of isotopically depleted (carbon isotope) ``polluted air'' (PA) with isotopically enriched ``free atmospheric air'' (FA) above. Using this modeling approach, the contribution of FA at ground level is being estimated for both the morning and afternoon time air samples.

Bispectra of A Tropical Coupled Ocean-Atmosphere System

It has recently been proposed that the broad spectrum of interannual variability in the tropics with a peak around four years results from an interaction between the linear low-frequency oscillatory mode of the coupled system and the nonlinear higher-frequency modes of the system. In this study we determine the bispectrum of the conceptual model consisting of a nonlinear low-order model coupled to a linear oscillator for various values of the coupling constants.

Lyapunov exponents and predictability of the tropical coupled ocean-atmosphere system

It has recently been proposed that the broad spectrum of interannual variability in the tropics with a peak around four years results from an interaction between the linear low-frequency oscillatory mode of the coupled system and the nonlinear higher-frequency modes of the system. In this study we determine the Lyapunov exponents of the conceptual model consisting of a nonlinear low-order model coupled to a linear oscillator for various values of the coupling constants.

Vertical thermal structure of the atmosphere during active and weak phases of convection over the north Bay of Bengal: Observation and model results

Owing to the lack of atmospheric vertical profile data with sufficient accuracy and vertical resolution, the response of the deep atmosphere to passage of monsoon systems over the Bay of Bengal. had not been satisfactorily elucidated. Under the Indian Climate Research Programme, a special observational programme called 'Bay of Bengal Monsoon Experiment' (BOBMEX), was conducted during July-August 1999. The present study is based on the high-resolution radiosondes launched during BOBMEX in the north Bay. Clear changes in the vertical thermal structure of the atmosphere between active and weak phases of convection have been observed. The atmosphere cooled below 6 km height and became warmer between 6 and 13 km height. The warmest layer was located between 8 and 10 km height, and the coldest layer was found just below 5 km height. The largest fluctuations in the humidity field occurred in the mid-troposphere. The observed changes between active and weak phases of convection are compared with the results from an atmospheric general circulation model, which is similar to that used at the National Centre for Medium Range Weather Forecasting, New Delhi. The model is not able to capture realistically some important features of the temperature and humidity profiles in the lower troposphere and in the boundary layer during the active and weak spells.

Some salient features of the atmosphere observed over the north Bay of Bengal during BOBMEX

This paper describes the near surface characteristics and vertical variations based on the observations made at 17.5degreesN and 89degreesE from ORV Sagar Kanya in the north Bay of Bengal during the Bay of Bengal Monsoon Experiment (BOBMEX) carried out in July-August 1999. BOBMEX captured both the active and weak phases of convection. SST remained above the convection threshold throughout the BOBMEX. While the response of the SST to atmospheric forcing was clearly observed, the response of the atmosphere to SST changes was not clear. SST decreased during periods of large scale precipitation, and increased during a weak phase of convection. It is shown that the latent heat flux at comparable wind speeds was about 25-50% lower over the Bay during BOBMEX compared to that over the Indian Ocean during other seasons and tropical west Pacific. On the other hand, the largest variations in the surface daily net heat flux are observed over the Bay during BOBMEX. SST predicted using observed surface fluxes showed that 1-D heat balance model works sometime but not always, and horizontal advection is important. The high resolution Vaisala radiosondes launched during BOBMEX could clearly bring out the changes in the vertical structure of the atmosphere between active and weak phases of convection. Convective Available Potential Energy of the surface air decreased,by 2-3 kJ kg(-1) following convection, and recovered in a time period of one or two days. The mid tropospheric relative humidity and water vapor content, and wind direction show the major changes between the active and weak phases of convection.

Numerical analysis of transpiration cooling in carbon dioxide atmosphere at hypersonic Mach numbers

The numerical solutions are obtained for skin friction, heat transfer to the wall and growth of boundary layer along the flat plate by employing two dimensional Navier-Stokes equations governing the hypersonic flow coupled with species continuity equations. Flow fields have been computed along the flat plate in CO2 atmosphere in the presence of transpiration cooling using air and carbon dioxide.