Incorporation of a physically based melt pond scheme into the sea ice component of a climate model

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2005 and 2007. These minima have not been predicted in the IPCC AR4 report, suggesting that the sea ice component of climate models should more realistically represent the processes controlling the sea ice mass balance. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds accumulate on the surface of sea ice from snow and sea ice melt and their presence reduces the albedo of the ice cover, leading to further melt. Toward the end of the melt season, melt ponds cover up to 50% of the sea ice surface. We have developed a melt pond evolution theory. Here, we have incorporated this melt pond theory into the Los Alamos CICE sea ice model, which has required us to include the refreezing of melt ponds. We present results showing that the presence, or otherwise, of a representation of melt ponds has a significant effect on the predicted sea ice thickness and extent. We also present a sensitivity study to uncertainty in the sea ice permeability, number of thickness categories in the model representation, meltwater redistribution scheme, and pond albedo. We conclude with a recommendation that our melt pond scheme is included in sea ice models, and the number of thickness categories should be increased and concentrated at lower thicknesses.

A continuum model of melt pond evolution on Arctic sea ice

[1] During the Northern Hemisphere summer, absorbed solar radiation melts snow and the upper surface of Arctic sea ice to generate meltwater that accumulates in ponds. The melt ponds reduce the albedo of the sea ice cover during the melting season, with a significant impact on the heat and mass budget of the sea ice and the upper ocean. We have developed a model, designed to be suitable for inclusion into a global circulation model (GCM), which simulates the formation and evolution of the melt pond cover. In order to be compatible with existing GCM sea ice models, our melt pond model builds upon the existing theory of the evolution of the sea ice thickness distribution. Since this theory does not describe the topography of the ice cover, which is crucial to determining the location, extent, and depth of individual ponds, we have needed to introduce some assumptions. We describe our model, present calculations and a sensitivity analysis, and discuss our results.

A model of melt pond evolution on sea ice

A one-dimensional, thermodynamic, and radiative model of a melt pond on sea ice is presented that explicitly treats the melt pond as an extra phase. A two-stream radiation model, which allows albedo to be determined from bulk optical properties, and a parameterization of the summertime evolution of optical properties, is used. Heat transport within the sea ice is described using an equation describing heat transport in a mushy layer of a binary alloy (salt water). The model is tested by comparison of numerical simulations with SHEBA data and previous modeling. The presence of melt ponds on the sea ice surface is demonstrated to have a significant effect on the heat and mass balance. Sensitivity tests indicate that the maximum melt pond depth is highly sensitive to optical parameters and drainage. INDEX TERMS: 4207 Oceanography: General: Arctic and Antarctic oceanography; 4255 Oceanography: General: Numerical modeling; 4299 Oceanography: General: General or miscellaneous; KEYWORDS: sea ice, melt pond, albedo, Arctic Ocean, radiation model, thermodynamic

September Arctic sea-ice minimum predicted by spring melt-pond fraction

The area of Arctic September sea ice has diminished from about 7 million km2 in the 1990s to less than 5 million km2 in five of the past seven years, with a record minimum of 3.6 million km2 in 2012 (ref. 1). The strength of this decrease is greater than expected by the scientific community, the reasons for this are not fully understood, and its simulation is an on-going challenge for existing climate models2, 3. With growing Arctic marine activity there is an urgent demand for forecasting Arctic summer sea ice4. Previous attempts at seasonal forecasts of ice extent were of limited skill5, 6, 7, 8, 9. However, here we show that the Arctic sea-ice minimum can be accurately forecasted from melt-pond area in spring. We find a strong correlation between the spring pond fraction and September sea-ice extent. This is explained by a positive feedback mechanism: more ponds reduce the albedo; a lower albedo causes more melting; more melting increases pond fraction. Our results help explain the acceleration of Arctic sea-ice decrease during the past decade. The inclusion of our new melt-pond model10 promises to improve the skill of future forecast and climate models in Arctic regions and beyond.

Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice

Over Arctic sea ice, pressure ridges and floe andmelt pond edges all introduce discrete obstructions to the flow of air or water past the ice and are a source of form drag. In current climate models form drag is only accounted for by tuning the air–ice and ice–ocean drag coefficients, that is, by effectively altering the roughness length in a surface drag parameterization. The existing approach of the skin drag parameter tuning is poorly constrained by observations and fails to describe correctly the physics associated with the air–ice and ocean–ice drag. Here, the authors combine recent theoretical developments to deduce the total neutral form drag coefficients from properties of the ice cover such as ice concentration, vertical extent and area of the ridges, freeboard and floe draft, and the size of floes and melt ponds. The drag coefficients are incorporated into the Los Alamos Sea Ice Model (CICE) and show the influence of the new drag parameterization on the motion and state of the ice cover, with the most noticeable being a depletion of sea ice over the west boundary of the Arctic Ocean and over the Beaufort Sea. The new parameterization allows the drag coefficients to be coupled to the sea ice state and therefore to evolve spatially and temporally. It is found that the range of values predicted for the drag coefficients agree with the range of values measured in several regions of the Arctic. Finally, the implications of the new form drag formulation for the spinup or spindown of the Arctic Ocean are discussed.

Processes controlling surface, bottom and lateral melt of Arctic sea ice in a state of the art sea ice model

We present a modelling study of processes controlling the summer melt of the Arctic sea ice cover. We perform a sensitivity study and focus our interest on the thermodynamics at the ice–atmosphere and ice–ocean interfaces. We use the Los Alamos community sea ice model CICE, and additionally implement and test three new parametrization schemes: (i) a prognostic mixed layer; (ii) a three equation boundary condition for the salt and heat flux at the ice–ocean interface; and (iii) a new lateral melt parametrization. Recent additions to the CICE model are also tested, including explicit melt ponds, a form drag parametrization and a halodynamic brine drainage scheme. The various sea ice parametrizations tested in this sensitivity study introduce a wide spread in the simulated sea ice characteristics. For each simulation, the total melt is decomposed into its surface, bottom and lateral melt components to assess the processes driving melt and how this varies regionally and temporally. Because this study quantifies the relative importance of several processes in driving the summer melt of sea ice, this work can serve as a guide for future research priorities.

Dynamics of Aeromonas species isolated from wastewater treatment system

Aeromonas are widely distributed in the aquatic environment, and are considered to be emerging organisms that can produce a series of virulence factors. The present study was carried out in a sanitary sewage stabilization pond treatment system, located in Lins, State of Sao Paulo, Brazil. Most probable number was applied for estimation of the genus Aeromonas. Colony isolation was carried out on blood agar ampicillin and confirmed by biochemical characterization. Aeromonas species were isolated in 72.4% of influent samples, and in 55.2 and 48.3% of effluent from anaerobic and facultative lagoons, respectively. Thirteen Aeromonas species were isolated, representing most of the recognized species of these organisms. Even though it was possible to observe a tendency of decrease, total elimination of these organisms from the studied system was not achieved. Understanding of the pathogenic organism`s dynamics in wastewater treatment systems with a reuse potential is especially important because of the risk it represents.

The effects of stocking density in physiological parameters and growth of the endangered teleost species piabanha, Brycon insignis (Steindachner, 1877)

This study investigated the effects of stocking density on the growth and fatty acid (FA) of Brycon insignis metabolism. Fingerlings (360) were distributed into eight ponds at two stocking densities (105 and 210 g/m(3)). The analysis of growth showed that the condition factor (K) and the coefficient of variation (CV) for body mass were not affected by stocking density. However, final body mass and length, specific growth rate (SGR), and weight gain (WG) were higher in the low stocking density group, which also presented a higher feed efficiency (FE) and survival (S). By contrast, muscle protein levels were higher in the high stocking density group. The plasma and muscle lipid content were not affected by stocking density, but fish reared at lower stocking density presented higher lipid concentration in the liver, with no differences in hepatosomatic index values. Even with the differences observed in metabolic and growth parameters, plasma cortisol was not affected by stocking density. The FA profile in the muscle and liver neutral fraction were not affected by stocking density, but the FA in the polar fractions differed between the two stocking densities. In the liver, total polyunsaturated fatty acids (PUFA) and PUFA n - 3 increased in higher stocking density, mainly due to an increase in docosahexaenoic acid (DHA). In addition, PUFA n - 6 were also increased in the higher stocking density group, mainly due to an increase in arachidonic acid (AA) and docosadienoic acid (22:2n - 6). In the muscle polar fraction, the saturated fatty acids (SFA) and monounsaturated fatty acids (MUFA) decreased in the animals from the higher stocking density group, and this reduction was compensated by an increase in PUFA n - 3 and PUFA n - 6, mainly the FA with 20-22 carbons (20:4n - 6: 22:4:n - 6; 22:5n - 6, 22:5n - 3, and 22:6n - 3). A different profile was observed for the C18 PUFAs, mainly 18:2n - 6 and 18:4n - 6, which were higher in the lower density stocking group. The data suggest that when living in high stocking density, B. insignis differentially utilizes the hepatic lipids as energy source and remodels the membrane fatty acids, with higher amounts of DHA in the polar muscle fraction compensated for by a decrease in MUFA. The zootechnical and physiological indices reveal that the lower stocking density group achieve overall better performance. (C) 2010 Elsevier B.V. All rights reserved.

The importance of dietary calcium consumption in two species of semi-terrestrial grapsoid crabs

Calcium (Ca) is essential for crustaceans, due to calcium carbonate (CaCO(3)) deposition in the new exoskeleton to harden it. The purpose of this work was to study short term Ca balance in terms of dietary Ca ingestion in two phylogenetically related crabs (Superfamily Grapsoidea) showing different degrees of terrestrial adaptations: Sesarma rectum Randall, 1840 and Neohelice granulata (Dana, 1851). Dietary Ca ingestion was studied using purified diets with different Ca concentrations (0, 2.2 and 6.66% Ca), together with measurements of Ca excretion and Ca hemolymph levels. The results showed that both crabs had the same response to foods containing different levels of Ca, with both species eating more of the high Ca diet. However, S. rectum consumed more per mg body mass at all Ca concentrations (6 mg.g(-1) for S. rectum against 3 mg.g(-1) for N. granulata). Both species excreted/egested Ca differently: S. rectum excreted Ca proportionally to ingestion, whereas N. granulata maintained constant faecal Ca output at all dietary Ca levels. Moreover, Ca hemolymph levels for crabs fed the different diets were independent of dietary Ca. In conclusion, both S. rectum and N. granulata seem to regulate the consumption of diets containing more Ca, which suggests a fine balance for Ca intake.

Periodic window arising in the parameter space of an impact oscillator

In the bi-dimensional parameter space of an impact-pair system, shrimp-shaped periodic windows are embedded in chaotic regions. We show that a weak periodic forcing generates new periodic windows near the unperturbed one with its shape and periodicity. Thus, the new periodic windows are parameter range extensions for which the controlled periodic oscillations substitute the chaotic oscillations. We identify periodic and chaotic attractors by their largest Lyapunov exponents. (C) 2010 Elsevier B.V. All rights reserved.

The Itajai foreland basin: a tectono-sedimentary record of the Ediacaran period, Southern Brazil

The Itajai Basin located in the southern border of the Luis Alves Microplate is considered as a peripheral foreland basin related to the Dom Feliciano Belt. It presents an excellent record of the Ediacaran period, and its upper parts display the best Brazilian example of Precambrian turbiditic deposits. The basal succession of Itajai Group is represented by sandstones and conglomerates (BaA(0) Formation) deposited in alluvial and deltaic-fan systems. The marine upper sequences correspond to the Ribeiro Carvalho (channelized and non-channelized proximal silty-argillaceous rhythmic turbidites), Ribeiro Neisse (arkosic sandstones and siltites), and Ribeiro do Bode (distal silty turbidites) formations. The ApiA(0)na Formation felsic volcanic rocks crosscut the sedimentary succession. The Cambrian Subida leucosyenogranite represents the last felsic magmatic activity to affect the Itajai Basin. The Brusque Group and the Florianpolis Batholith are proposed as source areas for the sediments of the upper sequence. For the lower continental units the source areas are the Santa Catarina, So Miguel and CamboriA(0) complexes. The lack of any oceanic crust in the Itajai Basin suggests that the marine units were deposited in a restricted, internal sea. The sedimentation started around 600 Ma and ended before 560 Ma as indicated by the emplacement of rhyolitic domes. The Itajai Basin is temporally and tectonically correlated with the Camaqu Basin in Rio Grande do Sul and the Arroyo del Soldado/Piriapolis Basin in Uruguay. It also has several tectono-sedimentary characteristics in common with the African-equivalent Nama Basin.

Chemostratigraphy of the Tamengo Formation (Corumba Group, Brazil) A contribution to the calibration of the Ediacaran carbon-isotope curve

The Corumba Group cropping out in the southern Paraguay Belt in Brazil is one of the most complete Ediacaran sedimentary archives of palaeogeographic climatic biogeochemical and biotic evolution in southwestern Gondwana The unit hosts a rich fossil record including acritarchs vendotaenids (Vendo taenia Eoholynia) soft-bodied metazoans (Corumbella) and skeletal fossils (Cloudina Titanotheca) The Tamengo Formation made up mainly of limestones and marls provides a rich bio- and chemostratigraphic record Several outcrops formerly assigned to the Cuiaba Group are here included in the Tamengo Formation on the basis of lithological and chemostratigraphical criteria High-resolution carbon isotopic analyses are reported for the Tamengo Formation showing (from base to top) (1) a positive delta(13)C excursion to +4 parts per thousand PDB above post-glacial negative values (2) a negative excursion to -3 5 parts per thousand associated with a marked regression and subsequent transgression (3) a positive excursion to +5 5 parts per thousand and (4) a plateau characterized by delta(13)C around +3 parts per thousand A U-Pb SHRIMP zircon age of an ash bed Interbedded in the upper part of the delta(13)C positive plateau yielded 543 +/- 3 Ma which is considered as the depositional age (Babinski et al 2008a) The positive plateau in the upper Tamengo Formation and the preceding positive excursion are ubiquitous features in several successions worldwide including the Nama Group (Namibia) the Dengying Formation (South China) and the Nafun and Ara groups (Oman) This plateau is constrained between 542 and 551 Ma thus consistent with the age of the upper Tamengo Formation The negative excursion of the lower Tamengo Formation may be correlated to the Shuram-Wonoka negative anomaly although delta(13)C values do not fall beyond -3 5 parts per thousand in the Brazilian sections Sedimentary breccias occur just beneath this negative excursion in the lower Tamengo Formation One possible interpretation of the origin of these breccias is a glacioeustatic sea-level fall but a tectonic interpretation cannot be completely ruled out Published by Elsevier B V

Cambrian mafic to felsic magmatism and its connections with transcurrent shear zones of the Borborema Province (NE Brazil): Implications for the late assembly of the West Gondwana

New U-Pb (SHRIMP) and (40)Ar/(39)Ar isotopic data of igneous rocks and mylonites of the Borborema Province (NE Brazil) show that a wide range of tectonothermal events affected the province during the transition from the Precambrian to the Cambrian. Concordant zircon U-Pb data constrained the crystallization age of mafic stocks, mafic to felsic dikes and granite batholiths between 548 and 533 Ma. These bodies were emplaced in a regional strain field combining extension and dextral shearing. The ductile shear deformation overprinted an older basement fabric to develop a low- to medium metamorphic grade vertical mylonite belt that cut the province in the E-W direction. Magnetic fabrics of the Cambrian batholiths determined by anisotropy of magnetic susceptibility are consistent with syntectonic emplacement. The magmatic pulses and shear deformation would have supplied enough heat to reset the synkinematic micas of mylonites to yield (40)Ar/(39)Ar plateau cooling ages between ca. 550 and 510 Ma. These results provide evidence that emplacement of Early Cambrian mafic and felsic magmas were accompanied by regional-scale shear deformations, probably in the consequence of late collisions along the West Gondwana margin. (C) 2010 Published by Elsevier B.V.

The Cariris Velhos tectonic event in Northeast Brazil

The Borborema Province in northeastern South America is a typical Brasiliano-Pan-African branching system of Neoproterozoic orogens that forms part of the Western Gondwana assembly. The province is positioned between the Sao Luis-West Africa craton to the north and the Sao Francisco (Congo-Kasai) craton to the south. For this province the main characteristics are (a) its subdivision into five major tectonic domains, bounded mostly by long shear zones, as follows: Medio Coreau, Ceara Central, Rio Grande do Norte, Transversal, and Southern; (b) the alternation of supracrustal belts with reworked basement inliers (Archean nuclei + Paleoproterozoic belts); and (c) the diversity of granitic plutonism, from Neoproterozoic to Early Cambrian ages, that affect supracrustal rocks as well as basement inliers. Recently, orogenic rock assemblages of early Tonian (1000-920 Ma) orogenic evolution have been recognized, which are restricted to the Transversal and Southern domains of the Province. Within the Transversal Zone, the Alto Pajeu terrane locally includes some remnants of oceanic crust along with island arc and continental arc rock assemblages, but the dominant supracrustal rocks are mature and immature pelitic metasedimentary and metavolcaniclastic rocks. Contiguous and parallel to the Alto Pajeu terrane, the Riacho Gravata subterrane consists mainly of low-grade metamorphic successions of metarhythmites, some of which are clearly turbiditic in origin, metaconglomerates, and sporadic marbles, along with interbedded metarhyolitic and metadacitic volcanic or metavolcaniclastic rocks. Both terrane and subterrane are cut by syn-contractional intrusive sheets of dominantly peraluminous high-K calc-alkaline, granititic to granodioritic metaplutonic rocks. The geochemical patterns of both supracrustal and intrusive rocks show similarities with associations of mature continental arc volcano-sedimentary sequences, but some subordinate intra-plate characteristics are also found. In both the Alto Pajeu and Riacho Gravata terranes, TIMS and SHRIMP U-Pb isotopic data from zircons from both metavolcanic and metaplutonic rocks yield ages between 1.0 and 0.92 Ga, which define the time span for an event of orogenic character, the Cariris Velhos event. Less extensive occurrences of rocks of Cariris Velhos age are recognized mainly in the southernmost domains of the Province, as for example in the Polo Redondo-Maranco terrane, where arc-affinity migmatite-granitic and meta-volcano-sedimentary rocks show U-Pb ages (SHRIMP data) around 0.98-0.97 Ga. For all these domains, Sm-Nd data exhibit Tom model ages between 1.9 and 1.1 Ga with corresponding slightly negative to slightly positive epsilon(Nd)(t) values. These domains, along with the Borborema Province as a whole, were significantly affected by tectonic and magmatic events of the Brasiliano Cycle (0.7-0.5 Ga), so that it is possible that there are some other early Tonian rock assemblages which were completely masked and hidden by these later Brasiliano events. Cariris Velhos processes are younger than the majority of orogenic systems at the end of Mesoproterozoic Era and beginning of Neoproterozoic throughout the world, e.g. Irumide belt, Kibaride belt and Namaqua-Natal belt, and considerably younger than those of the youngest orogenic process (Ottawan) in the Grenvillian System. Therefore, they were probably not associated with the proposed assembly of Rodinia. We suggest, instead, that Cariris Velhos magmatism and tectonism could have been related to a continental margin magmatic arc, with possible back-arc associations, and that this margin may have been a short-lived (<100 m.y.) leading edge of the newly assembled Rodinia supercontinent. (C) 2009 Elsevier Ltd. All rights reserved.

Late Paleozoic-Early Triassic magmatism on the western margin of Gondwana: Collahuasi area, Northern Chile

The basement in the `Altiplano` high plateau of the Andes of northern Chile mostly consists of late Paleozoic to Early Triassic felsic igneous rocks (Collahuasi Group) that were emplaced and extruded along the western margin of the Gondwana supercontinent. This igneous Suite crops out in the Collalluasi area and forms the backbone of most of the high Andes from latitude 20 degrees to 22 degrees S. Rocks of the Collahuasi Group and correlative formations form art extensive belt of volcanic and subvolcanic rocks throughout the main Andes of Chile, the Frontal Cordillera of Argentina (Choiyoi Group or Choiyoi Granite-Rhyolite Province), and the Eastern Cordillera of Peru. Thirteen new SHRIMP U-Pb zircon ages from the Collahuasi area document a bimodal timing for magnatism, with a dominant peak at about 300 Ma and a less significant one at 244 Ma. Copper-Mo porphyry mineralization is related to the younger igneous event. Initial Hf isotopic ratios for the similar to 300 Ma zircons range from about -2 to +6 indicating that the magmas incorporated components with a significant crustal residence time. The 244 Ma magmas were derived from a less enriched source, with the initial HT values ranging from +2 to +6, suggestive of a mixture with a more depleted component. Limited whole rock (144)Nd/(143)Nd and (87)Sr/(86)Sr isotopic ratios further support the likelihood that the Collahuasi Group magmatism incorporated significant older crustal components, or at least a mixture of crustal sources with more and less evolved isotopic signatures. (C) 2007 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.