@article{9152,
  abstract     = {Previous numerical studies of the dissipation of internal tides in idealized settings suggest the existence of a critical latitude (~29°) where dissipation is enhanced. But observations only indicate a modest enhancement at this latitude. To resolve this difference between observational and numerical results, the authors study the latitudinal dependence of internal tides’ dissipation in more realistic conditions. In particular, the ocean is not a quiescent medium; the presence of large-scale currents or mesoscale eddies can impact the propagation and dissipation of internal tides. This paper investigates the impact of a weak background mean current in numerical simulations. The authors focus on the local dissipation of high spatial mode internal waves near their generation site. The vertical profile of dissipation and its variation with latitude without the mean current are consistent with earlier studies. But adding a weak mean current has a major impact on the latitudinal distribution of dissipation. The peak at the critical latitude disappears, and the dissipation is closer to a constant, albeit with two weak peaks at ~25° and ~35° latitude. This disappearance results from the Doppler shift of the internal tides’ frequency, which hinders the nonlinear transfer of energy to small-scale secondary waves via the parametric subharmonic instability (PSI). The new two weak peaks correspond to the Doppler-shifted critical latitudes of the left- and right-propagating waves. The results are confirmed in simulations with simple sinusoidal topography. Thus, although nonlinear transfers via PSI are efficient at dissipating internal tides, the exact location of the dissipation is sensitive to large-scale oceanic conditions.},
  author       = {Richet, O. and Muller, Caroline J and Chomaz, J.-M.},
  issn         = {0022-3670},
  journal      = {Journal of Physical Oceanography},
  keywords     = {Oceanography},
  number       = {6},
  pages        = {1457--1472},
  publisher    = {American Meteorological Society},
  title        = {{Impact of a mean current on the internal tide energy dissipation at the critical latitude}},
  doi          = {10.1175/jpo-d-16-0197.1},
  volume       = {47},
  year         = {2017},
}

@article{12648,
  abstract     = {Distributed glacier melt models generally assume that the glacier surface consists of bare exposed ice and snow. In reality, many glaciers are wholly or partially covered in layers of debris that tend to suppress ablation rates. In this paper, an existing physically based point model for the ablation of debris-covered ice is incorporated in a distributed melt model and applied to Haut Glacier d'Arolla, Switzerland, which has three large patches of debris cover on its surface. The model is based on a 10 m resolution digital elevation model (DEM) of the area; each glacier pixel in the DEM is defined as either bare or debris-covered ice, and may be covered in snow that must be melted off before ice ablation is assumed to occur. Each debris-covered pixel is assigned a debris thickness value using probability distributions based on over 1000 manual thickness measurements. Locally observed meteorological data are used to run energy balance calculations in every pixel, using an approach suitable for snow, bare ice or debris-covered ice as appropriate. The use of the debris model significantly reduces the total ablation in the debris-covered areas, however the precise reduction is sensitive to the temperature extrapolation used in the model distribution because air near the debris surface tends to be slightly warmer than over bare ice. Overall results suggest that the debris patches, which cover 10% of the glacierized area, reduce total runoff from the glacierized part of the basin by up to 7%.},
  author       = {Reid, T. D. and Carenzo, M. and Pellicciotti, Francesca and Brock, B. W.},
  issn         = {0148-0227},
  journal      = {Journal of Geophysical Research: Atmospheres},
  keywords     = {Paleontology, Space and Planetary Science, Earth and Planetary Sciences (miscellaneous), Atmospheric Science, Earth-Surface Processes, Geochemistry and Petrology, Soil Science, Water Science and Technology, Ecology, Aquatic Science, Forestry, Oceanography, Geophysics},
  number       = {D18},
  publisher    = {American Geophysical Union},
  title        = {{Including debris cover effects in a distributed model of glacier ablation}},
  doi          = {10.1029/2012jd017795},
  volume       = {117},
  year         = {2012},
}

@article{12651,
  abstract     = {Temperature data from three Automatic Weather Stations and twelve Temperature Loggers are used to investigate the spatiotemporal variability of temperature over a glacier, its main atmospheric controls, the suitability of extrapolation techniques and their effect on melt modeling. We use data collected on Juncal Norte Glacier, central Chile, during one ablation season. We examine temporal and spatial variability in lapse rates (LRs), together with alternative statistical interpolation methods. The main control over the glacier thermal regime is the development of a katabatic boundary layer (KBL). Katabatic wind occurs at night and in the morning and is eroded in the afternoon. LRs reveal strong diurnal variability, with steeper LRs during the day when the katabatic wind weakens and shallower LRs during the night and morning. We suggest that temporally variable LRs should be used to account for the observed change. They tend to be steeper than equivalent constant LRs, and therefore result in a reduction in simulated melt compared to use of constant LRs when extrapolating from lower to higher elevations. In addition to the temporal variability, the temperature-elevation relationship varies also in space. Differences are evident between local LRs and including such variability in melt modeling affects melt simulations. Extrapolation methods based on the spatial variability of the observations after removal of the elevation trend, such as Inverse Distance Weighting or Kriging, do not seem necessary for simulations of gridded temperature data over a glacier.},
  author       = {Petersen, L. and Pellicciotti, Francesca},
  issn         = {0148-0227},
  journal      = {Journal of Geophysical Research: Atmospheres},
  keywords     = {Paleontology, Space and Planetary Science, Earth and Planetary Sciences (miscellaneous), Atmospheric Science, Earth-Surface Processes, Geochemistry and Petrology, Soil Science, Water Science and Technology, Ecology, Aquatic Science, Forestry, Oceanography, Geophysics},
  number       = {D23},
  publisher    = {American Geophysical Union},
  title        = {{Spatial and temporal variability of air temperature on a melting glacier: Atmospheric controls, extrapolation methods and their effect on melt modeling, Juncal Norte Glacier, Chile}},
  doi          = {10.1029/2011jd015842},
  volume       = {116},
  year         = {2011},
}

@article{9145,
  abstract     = {We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Γα (α stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation ε. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs.
As for the RSM results, we show that the Γα’s are functions of both Ri and Rρ (Richardson number and density ratio representing double diffusion, DD); the Γα are different for heat, salt and momentum; in the case of heat, the traditional value Γh = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Γh can be three times as large, a result that we reproduce. The salt Γs is given in terms of Γh. The momentum Γm has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Γm(Ri, Rρ). Having tested Γh, we then test the momentum Γm by showing that the turbulent Prandtl number Γm/Γh vs. Ri reproduces the available data quite well.

As for the dissipation ε, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean’s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable εN−2 from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Γα = 0.2 with the RSM result that brings into the problem the Ri, Rρ dependence of the Γα; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies.
We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Γh(Ri, Rρ). The RSM model reproduces well the NATRE data. Bimodal ε-distribution. NATRE data show that ε(Ri < 1) ≈ 10ε(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri ≫ 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.},
  author       = {Canuto, V.M. and Howard, A.M. and Cheng, Y. and Muller, Caroline J and Leboissetier, A. and Jayne, S.R.},
  issn         = {1463-5003},
  journal      = {Ocean Modelling},
  keywords     = {Computer Science (miscellaneous), Geotechnical Engineering and Engineering Geology, Atmospheric Science, Oceanography},
  number       = {3-4},
  pages        = {70--91},
  publisher    = {Elsevier},
  title        = {{Ocean turbulence, III: New GISS vertical mixing scheme}},
  doi          = {10.1016/j.ocemod.2010.04.006},
  volume       = {34},
  year         = {2010},
}

@article{9147,
  abstract     = {As part of an ongoing effort to develop a parameterization of wave-induced abyssal mixing, the authors derive an heuristic model for nonlinear wave breaking and energy dissipation associated with internal tides. Then the saturation and dissipation of internal tides for idealized and observed topography samples are investigated. One of the main results is that the wave-induced mixing could be more intense and more confined to the bottom than previously assumed in numerical models. Furthermore, in this model wave breaking and mixing clearly depend on the small scales of the topography below 10 km or so, which is below the current resolution of global bathymetry. This motivates the use of a statistical approach to represent the unresolved topography when addressing the role of internal tides in mixing the deep ocean.},
  author       = {Muller, Caroline J and Bühler, Oliver},
  issn         = {1520-0485},
  journal      = {Journal of Physical Oceanography},
  keywords     = {Oceanography},
  number       = {9},
  pages        = {2077--2096},
  publisher    = {American Meteorological Society},
  title        = {{Saturation of the internal tides and induced mixing in the abyssal ocean}},
  doi          = {10.1175/2009jpo4141.1},
  volume       = {39},
  year         = {2009},
}

@article{12658,
  abstract     = {[1] During the ablation period 2001 a glaciometeorological experiment was carried out on Haut Glacier d'Arolla, Switzerland. Five meteorological stations were installed on the glacier, and one permanent automatic weather station in the glacier foreland. The altitudes of the stations ranged between 2500 and 3000 m a.s.l., and they were in operation from end of May to beginning of September 2001. The spatial arrangement of the stations and temporal duration of the measurements generated a unique data set enabling the analysis of the spatial and temporal variability of the meteorological variables across an alpine glacier. All measurements were taken at a nominal height of 2 m, and hourly averages were derived for the analysis. The wind regime was dominated by the glacier wind (mean value 2.8 m s−1) but due to erosion by the synoptic gradient wind, occasionally the wind would blow up the valley. A slight decrease in mean 2 m air temperatures with altitude was found, however the 2 m air temperature gradient varied greatly and frequently changed its sign. Mean relative humidity was 71% and exhibited limited spatial variation. Mean incoming shortwave radiation and albedo both generally increased with elevation. The different components of shortwave radiation are quantified with a parameterization scheme. Resulting spatial variations are mainly due to horizon obstruction and reflections from surrounding slopes, i.e., topography. The effect of clouds accounts for a loss of 30% of the extraterrestrial flux. Albedos derived from a Landsat TM image of 30 July show remarkably constant values, in the range 0.49 to 0.50, across snow covered parts of the glacier, while albedo is highly spatially variable below the zone of continuous snow cover. These results are verified with ground measurements and compared with parameterized albedo. Mean longwave radiative fluxes decreased with elevation due to lower air temperatures and the effect of upper hemisphere slopes. It is shown through parameterization that this effect would even be more pronounced without the effect of clouds. Results are discussed with respect to a similar study which has been carried out on Pasterze Glacier (Austria). The presented algorithms for interpolating, parameterizing and simulating variables and parameters in alpine regions are integrated in the software package AMUNDSEN which is freely available to be adapted and further developed by the community.},
  author       = {Strasser, Ulrich and Corripio, Javier and Pellicciotti, Francesca and Burlando, Paolo and Brock, Ben and Funk, Martin},
  issn         = {0148-0227},
  journal      = {Journal of Geophysical Research: Atmospheres},
  keywords     = {Paleontology, Space and Planetary Science, Earth and Planetary Sciences (miscellaneous), Atmospheric Science, Earth-Surface Processes, Geochemistry and Petrology, Soil Science, Water Science and Technology, Ecology, Aquatic Science, Forestry, Oceanography, Geophysics},
  number       = {D3},
  publisher    = {American Geophysical Union},
  title        = {{Spatial and temporal variability of meteorological variables at Haut Glacier d'Arolla (Switzerland) during the ablation season 2001: Measurements and simulations}},
  doi          = {10.1029/2003jd003973},
  volume       = {109},
  year         = {2004},
}

