On eddy viscosity, energy cascades, and the horizontal resolution of gridded satellite altimeter products

Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 283–300, doi:10.1175/JPO-D-11-0240.1.Motivated by the recent interest in ocean energetics, the widespread use of horizontal eddy viscosity in models, and the promise of high horizontal resolution data from the planned wide-swath satellite altimeter, this paper explores the impacts of horizontal eddy viscosity and horizontal grid resolution on geostrophic turbulence, with a particular focus on spectral kinetic energy fluxes Π(K) computed in the isotropic wavenumber (K) domain. The paper utilizes idealized two-layer quasigeostrophic (QG) models, realistic high-resolution ocean general circulation models, and present-generation gridded satellite altimeter data. Adding horizontal eddy viscosity to the QG model results in a forward cascade at smaller scales, in apparent agreement with results from present-generation altimetry. Eddy viscosity is taken to roughly represent coupling of mesoscale eddies to internal waves or to submesoscale eddies. Filtering the output of either the QG or realistic models before computing Π(K) also greatly increases the forward cascade. Such filtering mimics the smoothing inherent in the construction of present-generation gridded altimeter data. It is therefore difficult to say whether the forward cascades seen in present-generation altimeter data are due to real physics (represented here by eddy viscosity) or to insufficient horizontal resolution. The inverse cascade at larger scales remains in the models even after filtering, suggesting that its existence in the models and in altimeter data is robust. However, the magnitude of the inverse cascade is affected by filtering, suggesting that the wide-swath altimeter will allow a more accurate determination of the inverse cascade at larger scales as well as providing important constraints on smaller-scale dynamics.BKA received support from Office of Naval Research Grant N00014-11-1-0487, National Science Foundation (NSF) Grants OCE-0924481 and OCE- 09607820, and University of Michigan startup funds. KLP acknowledges support from Woods Hole Oceanographic Institution bridge support funds. RBS acknowledges support from NSF grants OCE-0960834 and OCE-0851457, a contract with the National Oceanography Centre, Southampton, and a NASA subcontract to Boston University. JFS and JGR were supported by the projects ‘‘Global and remote littoral forcing in global ocean models’’ and ‘‘Agesotrophic vorticity dynamics of the ocean,’’ respectively, both sponsored by the Office of Naval Research under program element 601153N.2013-08-0

Effects of stencil width on surface ocean geostrophic velocity and vorticity estimation from gridded satellite altimeter data

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/90421/1/2011JC007367.pd

Effects of stencil width on surface ocean geostrophic velocity and vorticity estimation from gridded satellite altimeter data

This paper examines the effect of “stencil width” on surface ocean geostrophic velocity and vorticity estimated from differentiating gridded satellite altimeter sea surface height products. In oceanographic applications, the value of the first derivative at a central grid point is generally obtained by differencing the sea surface heights at adjacent grid points. This is called a “three-point stencil centered difference”. Here the stencil width is increased from three to five, seven, and nine points, using well-known formulae from the numerical analysis literature. The discrepancies between velocities computed with successive stencils decreases with increasing stencil width, suggesting that wide stencil results are more reliable. Significant speed-dependent biases (up to 10–20%) are found between results computed from three-point stencils versus those computed from wider stencils. The geostrophic velocity, and the variance of geostrophic velocity, are underestimated with thin stencils. Similar results are seen in geostrophic velocities computed from high-resolution model output. In contrast to the case when three-point stencils are used, wider stencils yield estimates of the anisotropy of velocity variance that are insensitive to the differences in grid spacing between two widely used altimeter products. Three-point stencils yield incorrect anisotropies on the 1/4° anisotropic AVISO grid; we recommend the use of 7-point stencils. Despite the demonstrated inadequacies of the three-point stencils, the conclusions of earlier studies based on them, that the zonally averaged midlatitude eddy kinetic energy field is nearly isotropic, are found to pertain also with wider stencils. Finally, the paper also examines the strengths and limitations of applying noise-suppressing differentiators, versus classic centered differences, to altimeter data

Toward an internal gravity wave spectrum in global ocean models

Office of Naval Research, National Science FoundationPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/111824/1/grl_internalwavespectrum_2015.pd

Sexual selection on male vocal fundamental frequency in humans and other anthropoids

D.A.P. was supported by a National Institutes of Mental Health T32 MH70343-05 fellowship. J.R.W. was supported by a National Science Foundation predoctoral fellowship.In many primates, including humans, the vocalizations of males and females differ dramatically, with male vocalizations and vocal anatomy often seeming to exaggerate apparent body size. These traits may be favoured by sexual selection because low-frequency male vocalizations intimidate rivals and/or attract females, but this hypothesis has not been systematically tested across primates, nor is it clear why competitors and potential mates should attend to vocalization frequencies. Here we show across anthropoids that sexual dimorphism in fundamental frequency (F0) increased during evolutionary transitions towards polygyny, and decreased during transitions towards monogamy. Surprisingly, humans exhibit greater F0 sexual dimorphism than any other ape. We also show that low-F0 vocalizations predict perceptions of men’s dominance and attractiveness, and predict hormone profiles (low cortisol and high testosterone) related to immune function. These results suggest that low male F0 signals condition to competitors and mates, and evolved in male anthropoids in response to the intensity of mating competition.PostprintPeer reviewe

Non-linear regression models for Approximate Bayesian Computation

Approximate Bayesian inference on the basis of summary statistics is well-suited to complex problems for which the likelihood is either mathematically or computationally intractable. However the methods that use rejection suffer from the curse of dimensionality when the number of summary statistics is increased. Here we propose a machine-learning approach to the estimation of the posterior density by introducing two innovations. The new method fits a nonlinear conditional heteroscedastic regression of the parameter on the summary statistics, and then adaptively improves estimation using importance sampling. The new algorithm is compared to the state-of-the-art approximate Bayesian methods, and achieves considerable reduction of the computational burden in two examples of inference in statistical genetics and in a queueing model.Comment: 4 figures; version 3 minor changes; to appear in Statistics and Computin

A critical comparison of integral projection and matrix projection models for demographic analysis

Structured demographic models are among the most common and useful tools in population biology. However, the introduction of integral projection models (IPMs) has caused a profound shift in the way many demographic models are conceptualized. Some researchers have argued that IPMs, by explicitly representing demographic processes as continuous functions of state variables such as size, are more statistically efficient, biologically realistic, and accurate than classic matrix projection models, calling into question the usefulness of the many studies based on matrix models. Here, we evaluate how IPMs and matrix models differ, as well as the extent to which these differences matter for estimation of key model outputs, including population growth rates, sensitivity patterns, and life spans. First, we detail the steps in constructing and using each type of model. Second, we present a review of published demographic models, concentrating on size-based studies, which shows significant overlap in the way IPMs and matrix models are constructed and analyzed. Third, to assess the impact of various modeling decisions on demographic predictions, we ran a series of simulations based on size-based demographic data sets for five biologically diverse species. We found little evidence that discrete vital rate estimation is less accurate than continuous functions across a wide range of sample sizes or size classes (equivalently bin numbers or mesh points). Most model outputs quickly converged with modest class numbers (≥10), regardless of most other modeling decisions. Another surprising result was that the most commonly used method to discretize growth rates for IPM analyses can introduce substantial error into model outputs. Finally, we show that empirical sample sizes generally matter more than modeling approach for the accuracy of demographic outputs. Based on these results, we provide specific recommendations to those constructing and evaluating structured population models. Both our literature review and simulations question the treatment of IPMs as a clearly distinct modeling approach or one that is inherently more accurate than classic matrix models. Importantly, this suggests that matrix models, representing the vast majority of past demographic analyses available for comparative and conservation work, continue to be useful and important sources of demographic information.Support for this work was provided by NSF awards 1146489, 1242558, 1242355, 1353781, 1340024, 1753980, and 1753954, 1144807, 0841423, and 1144083. Support also came from USDA NIFA Postdoctoral Fellowship (award no. 2019-67012-29726/project accession no. 1019364) for R. K. Shriver; the Swiss Polar Institute of Food and Agriculture for N. I. Chardon; the ICREA under the ICREA Academia Programme for C. Linares; and SERDP contract RC-2512 and USDA National Institute of Food and Agriculture, Hatch project 1016746 for A .M. Louthan. This is Contribution no. 21-177-J from the Kansas Agricultural Experiment Station

Pan-Cancer Analysis of lncRNA Regulation Supports Their Targeting of Cancer Genes in Each Tumor Context

Long noncoding RNAs (lncRNAs) are commonly dys-regulated in tumors, but only a handful are known toplay pathophysiological roles in cancer. We inferredlncRNAs that dysregulate cancer pathways, onco-genes, and tumor suppressors (cancer genes) bymodeling their effects on the activity of transcriptionfactors, RNA-binding proteins, and microRNAs in5,185 TCGA tumors and 1,019 ENCODE assays.Our predictions included hundreds of candidateonco- and tumor-suppressor lncRNAs (cancerlncRNAs) whose somatic alterations account for thedysregulation of dozens of cancer genes and path-ways in each of 14 tumor contexts. To demonstrateproof of concept, we showed that perturbations tar-geting OIP5-AS1 (an inferred tumor suppressor) andTUG1 and WT1-AS (inferred onco-lncRNAs) dysre-gulated cancer genes and altered proliferation ofbreast and gynecologic cancer cells. Our analysis in-dicates that, although most lncRNAs are dysregu-lated in a tumor-specific manner, some, includingOIP5-AS1, TUG1, NEAT1, MEG3, and TSIX, synergis-tically dysregulate cancer pathways in multiple tumorcontexts

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smokin