Monte Carlo simulations of nucleation of colloidal crystals

A crystal nucleus in a finite volume may exhibit phase coexistence with a surrounding fluid. The thermodynamic properties of the coexisting fluid (pressure and chemical potential) are enhanced relative to their coexistence values. This enhancement is uniquely related to the surface excess free energy. rnA model for weakly attractive soft colloidal particles is investigated, the so called Asakura-Oosawa model. In simulations, this model allows for the calculation of the pressure in the liquid using the virial formula directly. The phase coexistence pressure in the thermodynamic limit is obtained from the interface velocity method. We introduce a method by which the chemical potential in dense liquids can be measured. There is neither a need to locate the interface nor to compute the anisotropic interfacial tension to obtain nucleation barriers. Therefore, our analysis is appropriate for nuclei of arbitrary shape. Monte Carlo simulations over a wide range of nucleus volumes yield to nucleation barriers independent from the total system volume. The interfacial tension is determined via the ensemble-switch method, hence a detailed test of classical nucleation theory is possible. The anisotropy of the interfacial tension and the resulting non-spherical shape has only a minor effect on the barrier for the Asakura-Oosawa model.

Perennial Snow and Ice Variations (2000-2008) in the Arctic Circumpolar Land Area from Satellite Observations

Perennial snow and ice (PSI) extent is an important parameter of mountain environments with regard to its involvement in the hydrological cycle and the surface energy budget. We investigated interannual variations of PSI in nine mountain regions of interest (ROI) between 2000 and 2008. For that purpose, a novel MODIS data set processed at the Canada Centre for Remote Sensing at 250 m spatial resolution was utilized. The extent of PSI exhibited significant interannual variations, with coefficients of variation ranging from 5% to 81% depending on the ROI. A strong negative relationship was found between PSI and positive degree-days (threshold 0°C) during the summer months in most ROIs, with linear correlation coefficients (r) being as low as r = −0.90. In the European Alps and Scandinavia, PSI extent was significantly correlated with annual net glacier mass balances, with r = 0.91 and r = 0.85, respectively, suggesting that MODIS-derived PSI extent may be used as an indicator of net glacier mass balances. Validation of PSI extent in two land surface classifications for the years 2000 and 2005, GLC-2000 and Globcover, revealed significant discrepancies of up to 129% for both classifications. With regard to the importance of such classifications for land surface parameterizations in climate and land surface process models, this is a potential source of error to be investigated in future studies. The results presented here provide an interesting insight into variations of PSI in several ROIs and are instrumental for our understanding of sensitive mountain regions in the context of global climate change assessment.

Tree species and diversity effects on soil water seepage in a tropical plantation

Plant diversity has been shown to influence the water cycle of forest ecosystems by differences in water consumption and the associated effects on groundwater recharge. However, the effects of biodiversity on soil water fluxes remain poorly understood for native tree species plantations in the tropics. Therefore, we estimated soil water fluxes and assessed the effects of tree species and diversity on these fluxes in an experimental native tree species plantation in Sardinilla (Panama). The study was conducted during the wet season 2008 on plots of monocultures and mixtures of three or six tree species. Rainfall and soil water content were measured and evapotranspiration was estimated with the Penman-Monteith equation. Soil water fluxes were estimated using a simple soil water budget model considering water input, output, and soil water and groundwater storage changes and in addition, were simulated using the physically based one-dimensional water flow model Hydrus-1D. In general, the Hydrus simulation did not reflect the observed pressure heads, in that modeled pressure heads were higher compared to measured ones. On the other hand, the results of the water balance equation (WBE) reproduced observed water use patterns well. In monocultures, the downward fluxes through the 200 cm-depth plane were highest below Hura crepitans (6.13 mm day−1) and lowest below Luehea seemannii (5.18 mm day−1). The average seepage rate in monocultures (±SE) was 5.66 ± 0.18 mm day−1, and therefore, significantly higher than below six-species mixtures (5.49 ± 0.04 mm day−1) according to overyielding analyses. The three-species mixtures had an average seepage rate of 5.63 ± 0.12 mm day−1 and their values did not differ significantly from the average values of the corresponding species in monocultures. Seepage rates were driven by the transpiration of the varying biomass among the plots (r = 0.61, p = 0.017). Thus, a mixture of trees with different growth rates resulted in moderate seepage rates compared to monocultures of either fast growing or slow growing tree species. Our results demonstrate that tree-species specific biomass production and tree diversity are important controls of seepage rates in the Sardinilla plantation during the wet season.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Trematomus bernacchii experimental data files

As a response to ocean warming, shifts in fish species distribution and changes in production have been reported that have been partly attributed to temperature effects on the physiology of animals. The Southern Ocean hosts some of the most rapidly warming regions on earth and Antarctic organisms are reported to be especially temperature sensitive. While cellular and molecular organismic levels appear, at least partially, to compensate for elevated temperatures, the consequences of acclimation to elevated temperature for the whole organism are often less clear. Growth and reproduction are the driving factors for population structure and abundance. The aim of this study was to assess the effect of long-term acclimation to elevated temperature on energy budget parameters in the high-Antarctic fish Trematomus bernacchii. Our results show a complete temperature compensation for routine metabolic costs after 9 weeks of acclimation to 4°C. However, an up to 84% reduction in mass growth was measured at 2 and 4°C compared with the control group at 0°C, which is best explained by reduced food assimilation rates at warmer temperatures. With regard to a predicted temperature increase of up to 1.4°C in the Ross Sea by 2200, such a significant reduction in growth is likely to affect population structures in nature, for example by delaying sexual maturity and reducing production, with severe impacts on Antarctic fish communities and ecosystems.

Respiration and ingestion rates of calanoid copepod in the northern Benguela Current System measured ex situ during the RSS Discovery D356 cruise in 2010

Respiration rates of 16 calanoid copepod species from the northern Benguela upwelling system were measured on board RRS Discovery in September/October 2010 to determine their energy requirements and assess their significance in the carbon cycle. Copepod species were sampled by different net types. Immediately after the hauls, samples were sorted to species and stages (16 species; females, males and C5 copepodids) according to Bradford-Grieve et al. (1999). Specimens were kept in temperature-controlled refrigerators for at least 12 h before they were used in experiments. Respiration rates of different copepod species were measured onboard by optode respirometry (for details see Köster et al., 2008) with a 10-channel optode respirometer (PreSens Precision Sensing Oxy-10 Mini, Regensburg, Germany) under simulated in situ conditions in temperature-controlled refrigerators. Experiments were run in gas-tight glass bottles (12-13 ml). For each set of experiments, two controls without animals were measured under exactly the same conditions to compensate for potential bias. The number of animals per bottle depended on the copepods size, stage and metabolic activity. Animals were not fed during the experiments but they showed natural species-specific movements. Immediately after the experiments, all specimens were deep-frozen at - 80 °C for later dry mass determination (after lyophilisation for 48 h) in the home lab. The carbon content (% of dry mass) of each species was measured by mass-spectrometry in association with stable isotope analysis and body dry mass was converted to units of carbon. For species without available carbon data, the mean value of all copepod species (44% dry mass) was applied. For the estimation of carbon requirements of copepod species, individual oxygen consumption rates were converted to carbon units, assuming that the expiration of 1 ml oxygen mobilises 0.44 mg of organic carbon by using a respiratory quotient (RQ) of 0.82 for a mixed diet consisting of proteins (RQ = 0.8-1.0), lipids (RQ = 0.7) and carbohydrates (RQ = 1.0) (Auel and Werner, 2003). The carbon ingestion rates were calculated using the energy budget and the potential maximum ingestion rate approach. To allow for physiological comparisons of respiration rates of deep- and shallow-living copepod species without the effects of ambient temperature and different individual body mass, individual respiration rates were temperature- (15°C, Q10=2) and size-adjusted. The scaling coefficient of 0.76 (R2=0.556) is used for the standardisation of body dry mass to 0.3 mg (mean dry mass of all analysed copepods), applying the allometric equation R= (R15°C/M0.76)×0.30.76, where R is respiration and M is individual dry mass in mg.

Combined effects of different CO2 levels and N sources on the diazotrophic cyanobacterium Trichodesmium

To predict effects of climate change and possible feedbacks, it is crucial to understand the mechanisms behind CO2 responses of biogeochemically relevant phytoplankton species. Previous experiments on the abundant N2 fixers Trichodesmium demonstrated strong CO2 responses, which were attributed to an energy reallocation between its carbon (C) and nitrogen (N) acquisition. Pursuing this hypothesis, we manipulated the cellular energy budget by growing Trichodesmium erythraeum IMS101 under different CO2 partial pressure (pCO2) levels (180, 380, 980 and 1400?µatm) and N sources (N2 and NO3-). Subsequently, biomass production and the main energy-generating processes (photosynthesis and respiration) and energy-consuming processes (N2 fixation and C acquisition) were measured. While oxygen fluxes and chlorophyll fluorescence indicated that energy generation and its diurnal cycle was neither affected by pCO2 nor N source, cells differed in production rates and composition. Elevated pCO2 increased N2 fixation and organic C and N contents. The degree of stimulation was higher for nitrogenase activity than for cell contents, indicating a pCO2 effect on the transfer efficiency from N2 to biomass. pCO2-dependent changes in the diurnal cycle of N2 fixation correlated well with C affinities, confirming the interactions between N and C acquisition. Regarding effects of the N source, production rates were enhanced in NO3-grown cells, which we attribute to the higher N retention and lower ATP demand compared with N2 fixation. pCO2 effects on C affinity were less pronounced in NO3- users than N2 fixers. Our study illustrates the necessity to understand energy budgets and fluxes under different environmental conditions for explaining indirect effects of rising pCO2.

Seawater carbonate chemistry, survival rate, shell mass growth, shell dissolution and locomotory speed of Limacina retroversa in a laboratory experiment

Anthropogenic carbon dioxide emissions induce ocean acidification, thereby reducing carbonate ion concentration, which may affect the ability of calcifying organisms to build shells. Pteropods, the main planktonic producers of aragonite in the worlds' oceans, may be particularly vulnerable to changes in sea water chemistry. The negative effects are expected to be most severe at high-latitudes, where natural carbonate ion concentrations are low. In this study we investigated the combined effects of ocean acidification and freshening on Limacina retroversa, the dominant pteropod in sub polar areas. Living L. retroversa, collected in Northern Norwegian Sea, were exposed to four different pH values ranging from the pre-industrial level to the forecasted end of century ocean acidification scenario. Since over the past half-century the Norwegian Sea has experienced a progressive freshening with time, each pH level was combined with a salinity gradient in two factorial, randomized experiments investigating shell degradation, swimming behavior and survival. In addition, to investigate shell degradation without any physiologic influence, one perturbation experiments using only shells of dead pteropods was performed. Lower pH reduced shell mass whereas shell dissolution increased with pCO2. Interestingly, shells of dead organisms had a higher degree of dissolution than shells of living individuals. Mortality of Limacina retroversa was strongly affected only when both pH and salinity reduced simultaneously. The combined effects of lower salinity and lower pH also affected negatively the ability of pteropods to swim upwards. Results suggest that the energy cost of maintaining ion balance and avoiding sinking (in low salinity scenario) combined with the extra energy cost necessary to counteract shell dissolution (in high pCO2 scenario), exceed the available energy budget of this organism causing the pteropods to change swimming behavior and begin to collapse. Since L. retroversa play an important role in the transport of carbonates to the deep oceans these findings have significant implications for the mechanisms influencing the inorganic carbon cycle in the sub-polar area.

Seawater carbonate chemistry and biological processes of Mytilus edulis during experiments, 2011

Progressive ocean acidification due to anthropogenic CO2 emissions will alter marine ecosytem processes. Calcifying organisms might be particularly vulnerable to these alterations in the speciation of the marine carbonate system. While previous research efforts have mainly focused on external dissolution of shells in seawater under saturated with respect to calcium carbonate, the internal shell interface might be more vulnerable to acidification. In the case of the blue mussel Mytilus edulis, high body fluid pCO2 causes low pH and low carbonate concentrations in the extrapallial fluid, which is in direct contact with the inner shell surface. In order to test whether elevated seawater pCO2 impacts calcification and inner shell surface integrity we exposed Baltic M. edulis to four different seawater pCO2 (39, 142, 240, 405 Pa) and two food algae (310-350 cells mL-1 vs. 1600-2000 cells mL-1) concentrations for a period of seven weeks during winter (5°C). We found that low food algae concentrations and high pCO2 values each significantly decreased shell length growth. Internal shell surface corrosion of nacreous ( = aragonite) layers was documented via stereomicroscopy and SEM at the two highest pCO2 treatments in the high food group, while it was found in all treatments in the low food group. Both factors, food and pCO2, significantly influenced the magnitude of inner shell surface dissolution. Our findings illustrate for the first time that integrity of inner shell surfaces is tightly coupled to the animals' energy budget under conditions of CO2 stress. It is likely that under food limited conditions, energy is allocated to more vital processes (e.g. somatic mass maintenance) instead of shell conservation. It is evident from our results that mussels exert significant biological control over the structural integrity of their inner shell surfaces.