Isotopic composition of water in fog, rain, throughfall and stemflow in the Monte Verde, Costa Rica

Fog deposition, precipitation, throughfall and stemflow were measured in a windward tropical montane cloud forest near Monteverde, Costa Rica, for a 65-day period during the dry season of 2003. Net fog deposition was measured directly using the eddy covariance (EC) method and it amounted to 1.2 ± 0.1 mm/day (mean ± standard error). Fog water deposition was 5-9% of incident rainfall for the entire period, which is at the low end of previously reported values. Stable isotope concentrations (d18O and d2H) were determined in a large number of samples of each water component. Mass balance-based estimates of fog deposition were 1.0 ± 0.3 and 5.0 ± 2.7 mm/day (mean ± SE) when d18O and d2H were used as tracer, respectively. Comparisons between direct fog deposition measurements and the results of the mass balance model using d18O as a tracer indicated that the latter might be a good tool to estimate fog deposition in the absence of direct measurement under many (but not all) conditions. At 506 mm, measured water inputs over the 65 days (fog plus rain) fell short by 46 mm compared to the canopy output of 552 mm (throughfall, stemflow and interception evaporation). This discrepancy is attributed to the underestimation of rainfall during conditions of high wind.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.