Contribution of Advective and Non-advective Heat Fluxes to the Heat Budget of a Shallow Lagoon

The heat budget in a shallow lagoon has been established from field measurements at a bihourly scale. Information on the main advective and non-advective heat fluxes were collected during year 2003 at Nueva lagoon (Almería, Southern Spain). Heat storage data was obtained from a thermistor chain located in the deepest part of the lagoon and meteorological information was acquired using an automatic meteorological station placed near the lagoon’s shore. In addition, estimation of evaporation was inferred from climatic approaches. Inputs of heat energy were dominated by radiative fluxes, with received net radiation accounting on average for around 95% of the non-advective total gains and radiation losses accounting for around 70% of the non-advective total losses. Sensible heat transfer from/to the atmosphere constituted the second energy input (4%) and output (20%), although heat losses by evaporation were also significant. Conduction of heat into the sediments was a relatively constant form of energy loss but constitutes a minor contribution on the overall heat budget. Considerable variability was evident in nonadvective heat fluxes at different time scales, from diel to seasonal. In relation to advective heat fluxes, groundwater and irrigation surpluses added to the heat storage of Nueva lagoon, whereas heat advected via precipitation was negligible.


INTRODUCTION
One of the major energy sinks of wetland ecosystems is latent heat flux, which is equivalent, in energy terms, to evaporation from a water surface. This term it is also one of the main components of the water balance in aquatic systems and, in fact, constitutes the main water sink in closed lakes. Consequently, an accurate estimation of evaporation rates will lead both to a validation of surface energy balance models and water budget studies. Methods used to measure or model evaporation vary considerably, with significant implications for the accuracy of the results [1] . The majority of these techniques requires high-quality instrumentation and are most applicable to short-term studies rather than long-term monitoring. In addition, surface energy fluxes vary considerably in a daily scale, so short-term data are required to achieve accurate estimations. Consequently, these techniques have not been widely applied in wetland environments during a long period of time, as several months or even a complete year.
As part of an integrated EU-LIFE project on the conservation and restoration of Adra lagoons (Andalusia, Spain), a set of high temporal resolution measurements (air temperature, temperature of the water column, relative air humidity, wind speed and solar radiation) were taken to establish the energy balance in Nueva lagoon during 2003. Latent heat flux data were calculated using the Penman approach in a daily scale from a nearby meteorological station in order to validate the model. Groundwater temperature was measured from a set of piezometers located around the whole perimeter of the lagoons. Energy budgets for coastal wetlands, may take into account the relationship between groundwater (coastal aquifer) and surface water (lagoon), as groundwater recharge could modify the energy inputs to the surface system [2] . Despite the importance of groundwater to the thermal regime of shallow lakes, there are very few studies on the fundamental heat fluxes that influence wetlands dominated by groundwater recharges. The present study reports a field investigation in which heat budgets have been quantified using detailed measurements (bihourly averages). Data were collected during a complete year in order to detect any seasonal contrast. The objectives of this study are (i) to evaluate two different climatic methods for evaporation estimates from a water surface, one of them integrated in an energy balance model (ii) to estimate energy interchange fluxes from a model using micrometeorological data and (iii) to calculate the water budget in the lagoon, in order to estimate fluxes due to advection in the system.

MATERIALS AND METHODS
Study site: The study lake is situated in the Southwest of the province of Almería (Spain), on the Adra River basin (746 km 2 of which the delta of the Adra River takes 22 km 2 ) (Fig. 1). The Nueva lagoon (3.5 m maximum depth; 29 ha surface) forms part of the Adra Lagoons Natural Reserve, a wetland located in the easternmost part of the delta. The whole delta is nowadays occupied with intensive agriculture (greenhouses), so the lagoon bears a great environmental impact due to lixiviation of pesticides and fertilizers through permeable materials. Several irrigation canals distribute water from the Adra River and pumped from the Adra River Delta aquifer. Surplus of irrigation runs into the Albufera de Adra lagoons, making up a considerable fraction of the water budget [3] . The lagoons are located over the Internal Zones of the Betic mountain range in the domain area of the Alpujárride Complex (Sierra de Gádor unit). These materials are basically Triassic phillites and quartzites over which lies a limestone-dolomite block [4] .
The area has an annual average temperature of 18°C, being the minimum in January (11°C) and the maximum in august (26°C). Many hours of sunshine and the protection of the surrounding sierras are the main factors involved in the warm average temperatures, which are the highest of the province. The annual precipitation is very low, around 400 mm on average and presents a typical Mediterranean regime with most of the rainfall occurring in autumn (102 mm) and winter (120 mm).
Heat budget: Daily evaporation was calculated using the Penman approach, which combined energy-budget and mass-balance methods to estimate evaporation from freewater surfaces, as follows: Where E is the daily evaporation rate (mm day 1 ), the slope of the saturation vapour pressure curve at the air temperature, the psychrometric constant and E a is the drying power of the air, given as a daily rate (mm day 1 ) by: Where (e w -e) is the saturation deficit, difference between the saturated (e w ) and the actual (e) vapour pressure (mbars) and f(u), the Penman's function of wind [5,1] . These data where acquired from Adra meteorological station (30S UG 500710.4066978, 42 m asl) in order to compare and validate evaporation results obtained by the heat budget methodology. Data from Adra station are public and accessible through a web-link [6] .
Heat budget was established by detailed monitoring of hydro-meteorological variables using in-lake equipment (a fixed buoy situated in the centre of the lagoon, to obtain water column temperatures from a chain of 10 submersed thermistors) and an adjacent to Nueva lagoon Squirrel MiniMet automatic station (Fig. 1).
The energy budget calculated for the lake, involve the vertical surface exchange that result from the portioning of incident solar (short wave) radiation (Qsa), long-wave radiation absorbed (Qla) and emitted (Qle), sensible heat flux, which is a convective energy flux in response to a temperature gradient (Qsh) and latent heat flux (Qec), which is a convective energy flux in response to a change of state among solid, liquid and gas phases. In the present heat budget, Qec is the energy gained by the lagoon through condensation (Qec > 0) or loosed by evaporation (Qec < 0). Conductive heat exchange with the stream bed (Qhb) was also calculated. Advective heat fluxes considered in the energy budget were heat fluxes due to precipitation (Qp), heat fluxes due to irrigation ditchs (Qi) and heat fluxes due to groundwater (Qg). Heat fluxes resulting from biological and chemical reactions were neglected due to morphology of the lake and time scale under consideration. Total net heat exchange (Qn), can be expressed as: Incident solar short-wave radiation (Qsa) absorbed by the water surface is the region of the spectrum between 0.3 and 3 m. The main factors affecting this variable are surface albedo and cloud cover, both reducing the theoretical total incident solar radiation. Bihourly data of incident short-wave radiation was directly measured using a radiometer.
Long-wave radiation absorbed by the lagoon (Qla) comes principally from water vapour, cloud droplets, carbon dioxide and ozone. The radiation is a function of the air temperature and was computed according to empirical relations.
The above equation is the Stefan-Boltzmann expression modified by several environmental factors. Here A = albedo, T a = absolute air temperature (ºK), the Stefan-Boltzmann constant and 0 the atmospheric emissivity, a function of the water vapour pressure in n a´' the atmosphere (%) and air temperature, computed following Idso [7] . Bihourly values of long wave radiation emitted from the lake (Qle) were computed taking into account that water radiates as a grey body rather than as a black one. Thus, Where emissivity ( ) = 0.96 and T w is the absolute surface water temperature determined at a depth of 0.5 m. Therefore, this term can be considered as an approximation of the lake surface temperature and this was also used to calculate the surface saturation water vapour pressure.
Finally, net radiation (Qr) can be obtained as follows: Energy used for evaporation/condensation, or latent heat flux, (Qec) is the result of heat budget itself as it measures the energy resulting from evaporation/ condensation processes. The sign indicates the energy transfer, being negative if evaporation predominates and vice versa. Empirical approximations [8,9] were used to compute this term of the heat budget: Qec = -f e (e w -e a ) (7) Here f e is an empirical transfer function [10] based on the air and water temperatures and the wind speed (u 12 ) 12 m above the ground (ms 1 ); f e = 4.8 + 1.92 u 12 + (T w -T a ) (8) e w and e a are the partial vapour pressure in the water and in the atmosphere, respectively and were calculated using Richards' equation [3] . The sensible heat flux (Qsh) is the energy flux resulting from molecular and turbulent conduction when there is a movement of the fluid (i.e. wind stress). Again, empirical relations were used to compute this variable: Where, the constant 0.63 is Bowen's ratio (mbar K 1 ).
Conductive heat exchange (Qhb) with the lagoon bed was calculated as: where K = thermal conductivity of bed material (J m 1 s 1°C 1 ), T g = temperature gradient in bed material (°Cm 1 ). Sensors were inserted at a depth of c.5 and 30 cm below the lagoon surface in order to determine the temperature gradient in the bed material. Results were then multiplied by the area of the lagoon's bed and units were changed to Wm 2 , as 1 cal cm 2 h 1 = 4.19 10 4 J m 2 h 1 = 11.63 Wm 2 .
More detailed information [11] about the aforementioned methodology is available in the literature.
Advective fluxes were obtained at a monthly scale by means of the estimation of the lagoon's water budget, given as: Where G i is groundwater inflow, G o is groundwater outflow, V is change in lagoon s water storage, P is precipitation, E is evaporation, S i is surface water inflow and S o is surface water outflow. As the lagoon is considered as terminal, S o cancels out. Surface water inflow has been periodically measured in the irrigation ditch flowing into the lagoon and daily evaporation and precipitation data were obtained from Adra weather station.
In order to obtain relative incidence of advective fluxes to the heat budget, change in storage of thermal energy in the water column was obtained. The total heat content (change in storage) of the water column [12] was quantified following the equation: Where Q t = change in heat content of the lake in cal cm 2 h 1 , A 0 = lagoon area (cm 2 ), z 0 = surface of the lagoon, z m = maximum depth, t z = average temperature in °C of a unit layer of water of thickness h z in cm, with the midpoint at depth z and A z = the area at depth z in cm 2 . Volume in cm 3 was then multiplied directly by temperature in °C to obtain caloric content. The hypsographic curve of Nueva lagoon was used to obtain A z . Information on the temperature at different depths (t z ) in the lake was obtained from a thermistors chain which encompassed 10 depths from 0.5 m to 2.5 m, with a constant interval of 20 cm between the probes. The thermistor chain was attached to a fixed buoy in the deepest point of the lagoon and the bihourly temperature data (maximum error in the readings was 0.05°C) were stored in a Squirrel Grant data logger.
There are many equations available for the estimation of these heat fluxes, especially for the Qsa, Qec and Qsh [13] ; choosing another equation may easily result in a difference of 20-50 Wm 2 in the heat flux estimate [14] . Nevertheless, the results provided using these techniques are in accordance with previous studies in other Spanish aquatic ecosystems [15,11] .

RESULTS
Air and water temperatures are the main variables that affect energy exchanges from water surface on an aquatic system. Plotting of average diel cycles for surface water in Nueva lagoon and air temperature in winter and summer (Fig. 2) reveals that water temperature has a modest diel variation compared with that of air temperature. Average daily range for surface water temperature was 0.59°C in the winter period and 1.82°C in the summer period. Daily ranges for air temperature in winter and summer were 10.16°C and 11.79°C, respectively. From the obtained results, a relative insensitivity of water temperature to changes in air temperature is observed in the lake, probably due to the groundwater recharge from the detritic aquifer.  (Table 1) shows that the variation of surface water temperature in the lagoon was much less  than that of air temperature. Mean water temperature was slightly higher (c.0.73°C) than mean air temperature in winter, spring and autumn, but nearly 1.5°C lower in summer.
Apart from water and air temperature, meteorological data involved in the heat budget can be observed in Table  2. Relative humidity was high during the whole year (c.70% on average) due to the proximity of the lagoon to the sea shore, no more than 250 m and high temperatures. Wind velocity was higher in January, April and December (> 2 ms 1 ) but didn't descend below 1 ms 1 during the rest of the year. Average typical deviation computed for wind speed was 0.7.
High solar radiation during the whole year is a key factor with a great deal of influence in the heat budget in any ecosystem; Nueva lagoon received a c.190 Wm 2 average radiation during 2003, being higher in June and July and lower in November and December.
Monthly averages for non-advective heat fluxes obtained in Nueva lagoon during the study period are shown in Table 3. In addition, Penman daily evaporation results obtained from Adra weather station are also shown. It can be observed that radiative fluxes, Qsa (+161 Wm 2 ), Qla (+270 Wm 2 ) and Qle (-403 Wm 2 ) are the most important components of the heat budget; the overall effect of non-advective radiative fluxes is to add heat energy to the lagoon, although this phenomenon is more pronounced in the summer. For the duration of four months of the year, radiation budget is negative.
Convection (Qec and Qsh) and conduction (Qhb) are, on an annual scale, major sinks of thermal energy of the system. The data shown in Table 3 depict that the net global balance of non-advective heat fluxes is negative, but change in thermal energy stored in the lagoon is not observed. Compensation of the energy balance is due to advective heat fluxes, mostly groundwater recharge and inflow from irrigation ditches, as will be discussed later on this study.
The Penman combination approach provides an annual evaporation rate of 1239.6 mm (Table 2), whereas the heat budget would lead to an annual evaporation rate of 1226.1. The energy balance and the Penman method yield mean annual evaporation rates that constitute a difference of 13.5 mm, that is 1.08 and 1.10% of the annual evaporation values. Rigorous uncertainty calculations would be unrealistic because all methods include assumptions on empirical equations whose associated errors are difficult or impossible to estimate. However, the consistency of the results can be discussed through sensitivity analysis of the two methods to their respective input variables: measured data (air temperature, relative humidity, wind speed and short-wave radiation) and the Bowen ratio.
Analysis of the influences of individual input variables on evaporation estimated for both methods (by varying their value by 10%) reflects that they are mainly sensitive to the short-wave radiation, which is the energy source of evaporation. The energy balance method is the most sensitive to this parameter. Air temperature is also a common variable to both methods. In the energy balance method it appears only through the calculation of longwave radiation and it does not significantly influence the evaporation rate. In the Penman method the sensitivity to temperature is larger. Uncertainties on wind speed data affect the Penman evaporation rate, but this variable has low influence. A 10% variation of the Bowen ratio, used in the energy balance method, only slightly changes the evaporation estimate. Assuming a 10% error in all the input variables, the cumulative errors in the evaporation rates would range between 10 and 20%, so both methods provide consistent inter-annual evaporation rates.
From the results obtained in the heat budget, evaporation/condensation heat fluxes have a net effect of removing energy from the lake during both the winter and the summer periods, although more strongly in the latter. Similar situation occurs with sensible heat. In any case, a small fraction of the overall heat budget is positive both for sensible heat and condensation.
The average absolute, percentage gain and losses and duration (fraction of the year) of non-advective heat energy in the study site are detailed in Table 4. Radiation was the major input and accounted for > 94% of the total non-advective energy gains. Condensation and bed conduction from the lagoon's floor were the minor inputs, both totalizing less than 0.5% of the total non-advective energy inputs. Nevertheless, some gain in heat was computed as a result of condensation as discussed earlier, c.7% of the duration of the year. Significant amounts of heat energy were added to the lagoon by sensible transfer, especially on the summer period, although total contribution on the heat budget was less than 2%. The average sensible heat gain was greater in both absolute and percentage terms for the summer period, when the cooling influence of groundwater flow was most marked and air temperature was higher than water temperature. Table 3: Monthly averages of non-advective heat fluxes and evaporation rates obtained from the energy balance and the Penman equation Heat losses from the lagoon were dominated by radiative fluxes, as a consequence of the warm climate of this region and the high radiation, although sensible heat transfer from the water body to the overlying air was more than 20% of the total losses during the study period. Evaporative losses were constant throughout the year (93% of the time of the year), enhanced by windy conditions. This component of the heat budget constitutes more than 5% of the global heat budget (Table 4). Energy losses from the lagoon s floor through conduction was not a very significant form of non-advective output in the study site, although a quasi-constant (99.5% of the time of the year) flux of c.4.5 Wm 2 was obtained. This situation is probably due to the covering of the floor of the lagoon by macrophyte vegetation.

Evaporation rate (mm) ---------------------------------------
Transfer of energy by mass transfer (advective fluxes) to the lagoon is produced mainly by direct rainfall and water inflows to the system, via surface (irrigation surplus) or groundwater recharge. The percentage contribution of rainfall, groundwater and ditch contribution to the heat storage of Nueva lagoon is listed in Table 5. Contribution of each component is presented for winter, spring, summer and autumn 2003 in order to analyze seasonal patterns. The stored energy evolution in the lake was obtained from bihourly data of water temperature from a thermistors chain and then multiplied by the volume of each horizontal layer to obtain total calorie content in the water column. A difference in the heat storage of c.4500 cal cm 1 was observed between summer (c.6500 cal cm 1 ) and winter (c.3000 cal cm 1 ) in Nueva lagoon, being the minimum (2311 cal cm 1 ) registered at 10.00 h on 1 February and the maximum (6854 cal cm 1 ) registered at 18.00 h on 15 August 2003. Average heat content in the lagoon was 4876 cal cm 1 . In addition, no change in heat storage was observed from 1st January 2003 to 31 December 2003.
The amount of heat added by rainfall falling onto the water surface was small, being zero in the summer due to null precipitation. On an annual scale, average heat contribution by rainfall was c. 2%. Groundwater inflow contributes more significantly to the heat storage of the lagoon, this influence was more pronounced in the summer period (c. 8% of total energy content in lake), when hydraulic gradients were higher. In this situation, groundwater discharge from the Delta of Adra River aquifer to the Mediterranean sea must be extremely small in that sector of the aquifer and the lagoon may behave as a groundwater sink, this process induced by high evaporation rates (159 mm month 1 in August 2003). Heat added by irrigation surplus was small in winter, but relatively constant in spring, summer and autumn (c. 6%) being a considerable source of heat energy to the study site.
In relation to the thermal response of the lake to the diel variations in meteorological conditions and the resulting energy budget, we have selected three defined thermal phases: pre-stratification (from January to June), stratification (from July to August) and poststratification/overturn (from September to December). Here, thermal stratification is defined as a temperature difference greater than 1.5 °C between surface and deep waters [16] . Figure 3 shows daily variations in nonadvective energy fluxes, energy stored, water-column temperature, wind speed and precipitation. Five-day periods have been selected in order to examine daily changes in each of the three thermal phases established.
During the pre-stratification phase, rainfall episodes could lead to a decrease on the air temperature and, consequently, to a reduction on the temperature of the water column of more than 3°C during a 24-h period, as can be observed in Fig. 3a. In addition, wind speed exhibited considerable diel variations, the highest wind speed value occurred in the interval from 12.00 to 14.00 h (14 ms 1 ). Heat flux due to net radiation shows a strong periodic variation, with maximum midday values ranging from nearly 0 Wm 2 (15 March) to more than 550 Wm 2 (20 March), affecting the changes in storage of thermal energy. Minimum night-time radiative fluxes remain fairly constant (c. -130 Wm 2 ) during the pre-stratification phase. In cloudy and windy periods the lake could lose thermal energy during the whole day (Fig. 3a, see 17  march). Diel cooling and heating effects of the upper layers of Nueva lagoon were observed during the stratification period (Fig. 3b). As a consequence, a secondary shallow thermocline appears from 10.00 to 14.00 h almost every day, contributing to daily increments of the water column stability in the system. The average value of wind speed during stratification was 1.2 ms 1 , but small-scale windy episodes (from 14.00 to 22.00 h, 15 July) with an average wind speed of c.10 ms 1 lead to the gradual extension of the warmer surface water to greater depths and decreased the water column stability. From 16 July to 20 July, both gains and losses of thermal energy from net radiation were higher than that obtained in the pre-stratification period, +600 and -180 Wm 2 , respectively. The consequence of the net energy exchange within the lake's surface resulted in changes in the heat energy storage which fits a sinusoidal curve. On a daily basis, the system was loosing energy to the atmosphere during night and early morning and gaining thermal energy during the sunlight period of the day. Wind speed conditions during the post-stratification period (September to December 2003) varied from those measured during the stratification period. Average wind-speed was c. 2.5 ms 1 , with maximum values rising  episodes is the source of colder water to the system that cools the water column. Energy storage reflects diel variations, but also an annual pattern, similar to that of net radiation absorbed or emitted by the lagoon, although with a delay of two months above 10 ms 1 (14.00 h on 18 and 19 October). The mixing effect of the wind on the lagoon is evident from the data obtained, so wind speed is the main factor involved in mixing processes in water surfaces. Heat gains from net radiation were also appreciably lower, maximum midday values ranging from 480 to 200 Wm 2 . Diel evolution of water column temperature during this phase was very homogeneous and the differences between the temperature of surface and bottom water was not higher than 1.5°C, so the water column can be considered as mixed. Average diel variation of non-advective heat fluxes in winter, spring, summer and autumn are depicted in Fig. 4. The net radiation changed from negative to positive values at c.6.30 h in summer, whereas this change occurs at c.9.30 h in winter. Change of net radiation from positive to negative values occurs at c.18.00 h in summer and at c.17.00 h in winter. Evaporative heat fluxes also show a diel cycle, at least in spring and summer, with maximum losses in the middle of the day. Conditions that promoted heat gain by condensation during short periods in the early morning and evening are not represented in Fig. 4, although they were occasionally detected within the study period. Such processes were not continuous in any season, so averaged values did not reproduce the situation in a seasonal scale. The diel changes of sensible heat flux have a similar pattern of values rising to a peak in the middle of the day, but there are differences within the different seasons. Heat gain through the air was greater in winter and took place for all the diel cycle, except for the period between c.20.00 and 8.00 h. The seasonal diel variation of bed conduction ranged from -1 to -3.5 Wm 2 during the study period in Nueva lagoon. The variation in the course of 24h periods was not very strong, with constant losses throughout the day. The overall energy balance in Wm -2 and energy storage evolution in Nueva lagoon (cal cm 1 ) for the study period is shown in Fig. 5. Drastic reductions in energy stored in the lagoon are related to rainfall episodes. Groundwater recharge associated with rainfall

DISCUSSION AND CONCLUSION
Annual water temperature in The Nueva lagoon exhibit a typical pattern for many coastal shallow lakes in the Mediterranean area [17] . In addition, groundwater recharge to the lagoon has also been reported [18,3] and constitutes a relatively considerable part both to the water and heat budgets. Published studies on heat budgets using detailed measurements with such a high temporal resolution are scarce and only confined for 24-h periods [16] . However, there are some literature available based on seasonal observations [14] or about heat budgets made at a high temporal resolution in groundwater dominated rivers with continuous measurements during several days in winter and summer seasons [2] . The scarceness of published studies about heat budgets in wetlands is not observed in oceanic research. A recent study about the depth of formation of tropical thermoclines and its relation with the ocean heat budget stresses out the relevance of the timescale involved in the process [19] . In Spain, an investigation by the authors on a high mountain lake (La Caldera, Sierra Nevada), where groundwater is a minor component of the water budget, revealed that short-wave radiation, long-wave radiation absorbed, long-wave radiation emitted, evaporation/condensation and sensible heat fluxes contributed on average 40, 34, 33.5, 36 and 28% respectively to non-advective heat budget. The measuring period (74 days) takes a complete ice-free episode [11] . Consideration of the data collected at two hour intervals reveals that the non-advective fluxes affecting the study site had a wide range of variation within the study period. Many of the main heat fluxes have a systematic diel fluctuation, whereas only heat gain due to radiation has also an annual fluctuation with a summer maximum. Besides, net radiation is a source of energy during the daylight hours, but is responsible for heat loss during the night-time when shortwave inputs from the sun ceases and back radiation of long-wave energy from the water surface dominates. Heat flux due to evaporation also shows a diel cycle with maximum looses in the midday and early afternoon. Heat loss due to evaporation is reduced at nigh time. The diel changes of sensible heat flux are different within different seasons, but also exhibited a similar pattern, with values rising to a peak in the middle of the day. Changes in heat fluxes through bed conduction show little variations in a daily scale. This is a minor component of the of the heat budget as the lagoon floor is totally covered with submersed vegetation and this can be a factor that reduce bed conduction.
In terms of advective fluxes, the minor role of precipitation confirmed results gathered in water courses [2] . Advection due to groundwater and irrigation surplus was a more constant and more significant source of heat in Nueva lagoon. No loss of groundwater to the Mediterranean Sea was recorded in the lagoon; this situation confirmed the result obtained in previous studies about the water budget in this system [18] . Although groundwater contributes heat to the lagoon in both winter and summer, temperature differences between inflows and the lagoon water column mean that this input has a warming effect in winter but a cooling one in summer. Diel energy fluxes from and to the lake determine the thermal/density stratification, which is the major controlling factor in the plankton dynamics. The proper knowledge of the physical processes occurring in the aquatic ecosystem and its biological coupling is relevant to properly understand the ecosystem structure and dynamics. In this sense, the study of temporal scales governing different physical processes is of vital importance. As the heat exchange from the water surface to the atmosphere is both positive and negative during a 24-h period, this kind of environmental modeling has to be applied with at least bihourly data.
The present study has revealed considerable fluctuations in heat energy between and within days and the influence of groundwater recharge on the heat budget in Nueva lagoon during year 2003, although there is a need for further and longer-term research involving biological variables on Nueva lagoon and other aquatic ecosystems. In this sense, future research directions take account on the study of migrating patterns in phytoplankton and zooplankton communities in Nueva lagoon with a bihourly scale in order to detect the influence of physical limnology in biological dynamics.