Refine
Year of publication
Document Type
- Article (17)
- Postprint (1)
- Working Paper (1)
Language
- English (19)
Is part of the Bibliography
- yes (19)
Keywords
A catalog of genetic loci associated with kidney function from analyses of a million individuals
(2019)
Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through transancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these,147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.
We report on the gamma-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average Large Area Telescope (LAT) gamma-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 +/- 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of two), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 +/- 0.14, and the softest one is 2.51 +/- 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3 GeV. In this paper, we also present the first results from the 4.5 month long multifrequency campaign (2009 March 15-August 1) on Mrk 501, which included the Very Long Baseline Array (VLBA), Swift, RXTE, MAGIC, and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The extensive radio to TeV data set from this campaign provides us with the most detailed spectral energy distribution yet collected for this source during its relatively low activity. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton (SSC) model. In the framework of this model, we find that the dominant emission region is characterized by a size less than or similar to 0.1 pc (comparable within a factor of few to the size of the partially resolved VLBA core at 15-43 GHz), and that the total jet power (similar or equal to 10(44) erg s(-1)) constitutes only a small fraction (similar to 10(-3)) of the Eddington luminosity. The energy distribution of the freshly accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3 GeV-10 TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20 GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks. We find that the ultrarelativistic electrons and mildly relativistic protons within the blazar zone, if comparable in number, are in approximate energy equipartition, with their energy dominating the jet magnetic field energy by about two orders of magnitude.
Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. However, most of the gamma-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available gamma-ray instrumentation was easier to perform.
Aims. Our goal is to characterize the source gamma-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity.
Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign.
Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density.
Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
The unidentified very-high-energy (VHE; E > 0.1 TeV) gamma -ray source, HESS J1826-130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady gamma -ray flux from HESS J1826-130, which appears extended with a half-width of 0.21 degrees +/- 0.02 <br /> (stat)degrees <br /> stat degrees +/- 0.05 <br /> (sys)degrees sys degrees . The source spectrum is best fit with either a power-law function with a spectral index Gamma = 1.78 +/- 0.10(stat) +/- 0.20(sys) and an exponential cut-off at 15.2 <br /> (+5.5)(-3.2) -3.2+5.5 TeV, or a broken power-law with Gamma (1) = 1.96 +/- 0.06(stat) +/- 0.20(sys), Gamma (2) = 3.59 +/- 0.69(stat) +/- 0.20(sys) for energies below and above E-br = 11.2 +/- 2.7 TeV, respectively. The VHE flux from HESS J1826-130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula, HESS J1825-137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826-130 VHE emission related to PSR J1826-1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826-130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to greater than or similar to 200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants, molecular clouds, and counterparts detected in radio, X-rays, and TeV energies.
Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature. Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays). Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of similar to 0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was similar to 3 CU, and the peak of the high-energy spectral component was found to be at similar to 2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices <2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance <3 sigma) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency-peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone SSC scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behavior seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.
Ligands incorporating a tetraazamacrocycle receptor, a 'click'-derived triazole and a 1,8-naphthalimide fluorophore have proven utility as probes for metal ions. Three new cyclam-based molecular probes are reported, in which a piperidinyl group has been introduced at the 4-position of the naphthalimide fluorophore. These compounds have been synthesized using the copper(I)-catalyzed azide-alkyne Huisgen cycloaddition and their photophysical properties studied in detail. The alkylamino group induces the expected red-shift in absorption and emission spectra relative to the simple naphthalimide derivatives and gives rise to extended fluorescence lifetimes in aqueous buffer. The photophysical properties of these systems are shown to be highly solvent-dependent. Screening the fluorescence responses of the new conjugates to a wide variety of metal ions reveals significant and selective fluorescence quenching in the presence of copper(II), yet no fluorescence enhancement with zinc(II) as observed previously for the simple naphthalimide derivatives. Reasons for this different behaviour are proposed. Cytotoxicity testing shows that these new cyclam-triazole-dye conjugates display little or no toxicity against either DLD-1 colon carcinoma cells or MDA-MB-231 breast carcinoma cells, suggesting a potential role for these and related systems in biological sensing applications.
This study pushes our understanding of research reliability by reproducing and replicating claims from 110 papers in leading economic and political science journals. The analysis involves computational reproducibility checks and robustness assessments. It reveals several patterns. First, we uncover a high rate of fully computationally reproducible results (over 85%). Second, excluding minor issues like missing packages or broken pathways, we uncover coding errors for about 25% of studies, with some studies containing multiple errors. Third, we test the robustness of the results to 5,511 re-analyses. We find a robustness reproducibility of about 70%. Robustness reproducibility rates are relatively higher for re-analyses that introduce new data and lower for re-analyses that change the sample or the definition of the dependent variable. Fourth, 52% of re-analysis effect size estimates are smaller than the original published estimates and the average statistical significance of a re-analysis is 77% of the original. Lastly, we rely on six teams of researchers working independently to answer eight additional research questions on the determinants of robustness reproducibility. Most teams find a negative relationship between replicators' experience and reproducibility, while finding no relationship between reproducibility and the provision of intermediate or even raw data combined with the necessary cleaning codes.
Statistical properties of a free-electron laser revealed by Hanbury Brown-Twiss interferometry
(2017)
We present a comprehensive experimental analysis of statistical properties of the self-amplified spontaneous emission free-electron laser (FEL) FLASH by means of Hanbury Brown and Twiss interferometry. The experiments were performed at FEL wavelengths of 5.5, 13.4, and 20.8 nm. We determined the second-order intensity correlation function for all wavelengths and different operation conditions of FLASH. In all experiments a high degree of spatial coherence (above 50%) was obtained. Our analysis performed in spatial and spectral domains provided us with the independent measurements of an average pulse duration of the FEL that were below 60 fs. To explain the complicated behavior of the second-order intensity correlation function we developed an advanced theoretical model that includes the presence of multiple beams and external positional jitter of the FEL pulses. By this analysis we determined that in one of the experiments external positional jitter was about 25% of the beam size. We envision that methods developed in our study will be used widely for analysis and diagnostics of FEL radiation.
Glacial-interglacial variations in CO2 and methane in polar ice cores have been attributed, in part, to changes in global wetland extent, but the wetland distribution before the Last Glacial Maximum (LGM, 21 ka to 18 ka) remains virtually unknown. We present a study of global peatland extent and carbon (C) stocks through the last glacial cycle (130 ka to present) using a newly compiled database of 1,063 detailed stratigraphic records of peat deposits buried by mineral sediments, as well as a global peatland model. Quantitative agreement between modeling and observations shows extensive peat accumulation before the LGM in northern latitudes (> 40 degrees N), particularly during warmer periods including the last interglacial (130 ka to 116 ka, MIS 5e) and the interstadial (57 ka to 29 ka, MIS 3). During cooling periods of glacial advance and permafrost formation, the burial of northern peatlands by glaciers and mineral sediments decreased active peatland extent, thickness, and modeled C stocks by 70 to 90% from warmer times. Tropical peatland extent and C stocks show little temporal variation throughout the study period. While the increased burial of northern peats was correlated with cooling periods, the burial of tropical peat was predominately driven by changes in sea level and regional hydrology. Peat burial by mineral sediments represents a mechanism for long-term terrestrial C storage in the Earth system. These results show that northern peatlands accumulate significant C stocks during warmer times, indicating their potential for C sequestration during the warming Anthropocene.
A novel surface coating with durable broad-spectrum antibacterial ability was prepared based on mussel inspired dendritic polyglycerol (MI-dPG) embedded with copper nanoparticles (Cu NPs). The functional surface coating is fabricated via a facile dip-coating process followed by in situ reduction of copper ions with a MI-dPG coating to introduce Cu NPs into the coating matrix. This coating has been demonstrated to possess efficient long-term antibacterial properties against Escherichia coli (E. coli), Staphylococcus aureus (S. aureus), and kanamycin-resistant E. coli through an "attract-kill-release" strategy. The synergistic antibacterial activity of the coating was shown by the combination of two functions of the contact killing, reactive oxygen species production and Cu ions released from the coating. Furthermore, this coating inhibited biofilm formation and showed good compatibility to eukaryotic cells. Thus, this newly developed Cu NP-incorporated MI-dPG surface coating may find potential application in the design of antimicrobial coating, such as implantable devices.
Bacterial pore-forming toxins compromise plasmalemmal integrity, leading to Ca2+ influx, leakage of the cytoplasm, and cell death. Such lesions can be repaired by microvesicular shedding or by the endocytic uptake of the injured membrane sites. Cells have at their disposal an entire toolbox of repair proteins for the identification and elimination of membrane lesions. Sphingomyelinases catalyze the breakdown of sphingomyelin into ceramide and phosphocholine. Sphingomyelin is predominantly localized in the outer leaflet, where it is hydrolyzed by acid sphingomyelinase (ASM) after lysosomal fusion with the plasma membrane. The magnesium-dependent neutral sphingomyelinase (NSM)-2 is found at the inner leaflet of the plasmalemma. Because either sphingomyelinase has been ascribed a role in the cellular stress response, we investigated their role in plasma membrane repair and cellular survival after treatment with the pore-forming toxins listeriolysin O (LLO) or pneumolysin (PLY). Jurkat T cells, in which ASM or NSM-2 was down-regulated [ASM knockdown (KD) or NSM-2 KD cells], showed inverse reactions to toxin-induced membrane damage: ASM KD cells displayed reduced toxin resistance, decreased viability, and defects in membrane repair. In contrast, the down-regulation of NSM-2 led to an increase in viability and enhanced plasmalemmal repair. Yet, in addition to the increased plasmalemmal repair, the enhanced toxin resistance of NSM-2 KD cells also appeared to be dependent on the activation of p38/MAPK, which was constitutively activated, whereas in ASM KD cells, the p38/MAPK activation was constitutively blunted.Schoenauer, R., Larpin, Y., Babiychuk, E. B., Drucker, P., Babiychuk, V. S., Avota, E., Schneider-Schaulies, S., Schumacher, F., Kleuser, B., Koffel, R., Draeger, A. Down-regulation of acid sphingomyelinase and neutral sphingomyelinase-2 inversely determines the cellular resistance to plasmalemmal injury by pore-forming toxins.
Landslides arc common in high-latitude forest ecosystems that have developed on permafrost. The most vulnerable areas in the permafrost territories of Siberia occur on the south-facing slopes of northern rivers, where they arc observed on about 20% of the total area of river slopes. Landslide disturbances will likely increase with climate change especially due to increasing summer-autumn precipitation. These processes are the most destructive natural disturbance agent and lead to the complete removal of pre-slide forest ecosystems (vegetation cover and soil). To evaluate postsliding ecosystem succession, we undertook integrated ecological research at landslides of different age classes along the Nizhnyaya Tunguska River and the Kochechum River (Tura, Krasnoyarsk region, Russia). Just after the event (at the one-year-old site), we registered a drop in soil respiration, a threefold lower microbial respiration rate, and a fourfold smaller mineral soil carbon and nitrogen stock at bare soil (melkozem) plots at the middle location of the site as compared with the non affected control site. The recovery of disturbed areas began with the re-establishment of plant cover and the following accumulation of an organic soil layer. During the 35-year succession (L1972), the accumulated layer (0 layer)at the oldest site contained similar C- and N stocks to those found at the control sites. However, the mineral soil C- and N stocks and the microbial biomass even of the oldest landslide area- did not reach the value of these parameters in control plots. Later, the soil respiration level and the eco-physiological status of soil microbiota also recovered due to these changes. This study demonstrates that the recovery after landslides in permafrost forests takes several decades. In addition, the degradation of permafrost due to landslides clearly hinders the accumulation of soil organic matter in the mineral soil. (C) 2018 Elsevier B.v. All rights reserved.
Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m(2)/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m(2) at follow-up among those with eGFRcrea 60 mL/min/1.73m(2) or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or (LARP4B). Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.
Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m(2)/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m(2) at follow-up among those with eGFRcrea 60 mL/min/1.73m(2) or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or (LARP4B). Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.
Geomagnetic indices are convenient quantities that distill the complicated physics of some region or aspect of near-Earth space into a single parameter. Most of the best-known indices are calculated from ground-based magnetometer data sets, such as Dst, SYM-H, Kp, AE, AL, and PC. Many models have been created that predict the values of these indices, often using solar wind measurements upstream from Earth as the input variables to the calculation. This document reviews the current state of models that predict geomagnetic indices and the methods used to assess their ability to reproduce the target index time series. These existing methods are synthesized into a baseline collection of metrics for benchmarking a new or updated geomagnetic index prediction model. These methods fall into two categories: (1) fit performance metrics such as root-mean-square error and mean absolute error that are applied to a time series comparison of model output and observations and (2) event detection performance metrics such as Heidke Skill Score and probability of detection that are derived from a contingency table that compares model and observation values exceeding (or not) a threshold value. A few examples of codes being used with this set of metrics are presented, and other aspects of metrics assessment best practices, limitations, and uncertainties are discussed, including several caveats to consider when using geomagnetic indices. Plain Language Summary One aspect of space weather is a magnetic signature across the surface of the Earth. The creation of this signal involves nonlinear interactions of electromagnetic forces on charged particles and can therefore be difficult to predict. The perturbations that space storms and other activity causes in some observation sets, however, are fairly regular in their pattern. Some of these measurements have been compiled together into a single value, a geomagnetic index. Several such indices exist, providing a global estimate of the activity in different parts of geospace. Models have been developed to predict the time series of these indices, and various statistical methods are used to assess their performance at reproducing the original index. Existing studies of geomagnetic indices, however, use different approaches to quantify the performance of the model. This document defines a standardized set of statistical analyses as a baseline set of comparison tools that are recommended to assess geomagnetic index prediction models. It also discusses best practices, limitations, uncertainties, and caveats to consider when conducting a model assessment.