Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990–2013: Quantifying the epidemiological transition

Background The Global Burden of Disease Study 2013 (GBD 2013) aims to bring together all available epidemiological data using a coherent measurement framework, standardised estimation methods, and transparent data sources to enable comparisons of health loss over time and across causes, age–sex groups, and countries. The GBD can be used to generate summary measures such as disability-adjusted life-years (DALYs) and healthy life expectancy (HALE) that make possible comparative assessments of broad epidemiological patterns across countries and time. These summary measures can also be used to quantify the component of variation in epidemiology that is related to sociodemographic development. Methods We used the published GBD 2013 data for age-specific mortality, years of life lost due to premature mortality (YLLs), and years lived with disability (YLDs) to calculate DALYs and HALE for 1990, 1995, 2000, 2005, 2010, and 2013 for 188 countries. We calculated HALE using the Sullivan method; 95% uncertainty intervals (UIs) represent uncertainty in age-specific death rates and YLDs per person for each country, age, sex, and year. We estimated DALYs for 306 causes for each country as the sum of YLLs and YLDs; 95% UIs represent uncertainty in YLL and YLD rates. We quantified patterns of the epidemiological transition with a composite indicator of sociodemographic status, which we constructed from income per person, average years of schooling after age 15 years, and the total fertility rate and mean age of the population. We applied hierarchical regression to DALY rates by cause across countries to decompose variance related to the sociodemographic status variable, country, and time. Findings Worldwide, from 1990 to 2013, life expectancy at birth rose by 6·2 years (95% UI 5·6–6·6), from 65·3 years (65·0–65·6) in 1990 to 71·5 years (71·0–71·9) in 2013, HALE at birth rose by 5·4 years (4·9–5·8), from 56·9 years (54·5–59·1) to 62·3 years (59·7–64·8), total DALYs fell by 3·6% (0·3–7·4), and age-standardised DALY rates per 100 000 people fell by 26·7% (24·6–29·1). For communicable, maternal, neonatal, and nutritional disorders, global DALY numbers, crude rates, and age-standardised rates have all declined between 1990 and 2013, whereas for non–communicable diseases, global DALYs have been increasing, DALY rates have remained nearly constant, and age-standardised DALY rates declined during the same period. From 2005 to 2013, the number of DALYs increased for most specific non-communicable diseases, including cardiovascular diseases and neoplasms, in addition to dengue, food-borne trematodes, and leishmaniasis; DALYs decreased for nearly all other causes. By 2013, the five leading causes of DALYs were ischaemic heart disease, lower respiratory infections, cerebrovascular disease, low back and neck pain, and road injuries. Sociodemographic status explained more than 50% of the variance between countries and over time for diarrhoea, lower respiratory infections, and other common infectious diseases; maternal disorders; neonatal disorders; nutritional deficiencies; other communicable, maternal, neonatal, and nutritional diseases; musculoskeletal disorders; and other non-communicable diseases. However, sociodemographic status explained less than 10% of the variance in DALY rates for cardiovascular diseases; chronic respiratory diseases; cirrhosis; diabetes, urogenital, blood, and endocrine diseases; unintentional injuries; and self-harm and interpersonal violence. Predictably, increased sociodemographic status was associated with a shift in burden from YLLs to YLDs, driven by declines in YLLs and increases in YLDs from musculoskeletal disorders, neurological disorders, and mental and substance use disorders. In most country-specific estimates, the increase in life expectancy was greater than that in HALE. Leading causes of DALYs are highly variable across countries. Interpretation Global health is improving. Population growth and ageing have driven up numbers of DALYs, but crude rates have remained relatively constant, showing that progress in health does not mean fewer demands on health systems. The notion of an epidemiological transition—in which increasing sociodemographic status brings structured change in disease burden—is useful, but there is tremendous variation in burden of disease that is not associated with sociodemographic status. This further underscores the need for country-specific assessments of DALYs and HALE to appropriately inform health policy decisions and attendant actions.

Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks in 188 countries, 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013

Background The Global Burden of Disease, Injuries, and Risk Factor study 2013 (GBD 2013) is the first of a series of annual updates of the GBD. Risk factor quantification, particularly of modifiable risk factors, can help to identify emerging threats to population health and opportunities for prevention. The GBD 2013 provides a timely opportunity to update the comparative risk assessment with new data for exposure, relative risks, and evidence on the appropriate counterfactual risk distribution. Methods Attributable deaths, years of life lost, years lived with disability, and disability-adjusted life-years (DALYs) have been estimated for 79 risks or clusters of risks using the GBD 2010 methods. Risk–outcome pairs meeting explicit evidence criteria were assessed for 188 countries for the period 1990–2013 by age and sex using three inputs: risk exposure, relative risks, and the theoretical minimum risk exposure level (TMREL). Risks are organised into a hierarchy with blocks of behavioural, environmental and occupational, and metabolic risks at the first level of the hierarchy. The next level in the hierarchy includes nine clusters of related risks and two individual risks, with more detail provided at levels 3 and 4 of the hierarchy. Compared with GBD 2010, six new risk factors have been added: handwashing practices, occupational exposure to trichloroethylene, childhood wasting, childhood stunting, unsafe sex, and low glomerular filtration rate. For most risks, data for exposure were synthesised with a Bayesian meta-regression method, DisMod-MR 2.0, or spatial-temporal Gaussian process regression. Relative risks were based on meta-regressions of published cohort and intervention studies. Attributable burden for clusters of risks and all risks combined took into account evidence on the mediation of some risks such as high body-mass index (BMI) through other risks such as high systolic blood pressure and high cholesterol. Findings All risks combined account for 57·2% (95% uncertainty interval [UI] 55·8–58·5) of deaths and 41·6% (40·1–43·0) of DALYs. Risks quantified account for 87·9% (86·5–89·3) of cardiovascular disease DALYs, ranging to a low of 0% for neonatal disorders and neglected tropical diseases and malaria. In terms of global DALYs in 2013, six risks or clusters of risks each caused more than 5% of DALYs: dietary risks accounting for 11·3 million deaths and 241·4 million DALYs, high systolic blood pressure for 10·4 million deaths and 208·1 million DALYs, child and maternal malnutrition for 1·7 million deaths and 176·9 million DALYs, tobacco smoke for 6·1 million deaths and 143·5 million DALYs, air pollution for 5·5 million deaths and 141·5 million DALYs, and high BMI for 4·4 million deaths and 134·0 million DALYs. Risk factor patterns vary across regions and countries and with time. In sub-Saharan Africa, the leading risk factors are child and maternal malnutrition, unsafe sex, and unsafe water, sanitation, and handwashing. In women, in nearly all countries in the Americas, north Africa, and the Middle East, and in many other high-income countries, high BMI is the leading risk factor, with high systolic blood pressure as the leading risk in most of Central and Eastern Europe and south and east Asia. For men, high systolic blood pressure or tobacco use are the leading risks in nearly all high-income countries, in north Africa and the Middle East, Europe, and Asia. For men and women, unsafe sex is the leading risk in a corridor from Kenya to South Africa. Interpretation Behavioural, environmental and occupational, and metabolic risks can explain half of global mortality and more than one-third of global DALYs providing many opportunities for prevention. Of the larger risks, the attributable burden of high BMI has increased in the past 23 years. In view of the prominence of behavioural risk factors, behavioural and social science research on interventions for these risks should be strengthened. Many prevention and primary care policy options are available now to act on key risks.

Ferrocene in conjugation with a Fischer carbene: Synthesis, NLO, and electrochemical behavior of a novel organometallic push-pull system

Using an iterative sequence of Wittig olefination, reduction, oxidation, and condensation of an active methylene group to carbonyl, it was possible to prepare a series of organometallic push-pull molecules [(CO)(5)M=C(OCH3)(-CH=CH-)(n)(C5H4)Fe(C5H5), M = W, Cr, n = 1-4] in which ferrocene is the donor element and a Fisher carbene moeity is the acceptor group. The molecular first hyperpolarizability beta was determined by hyper-Rayleigh scattering experiments. The beta values ranged from 110 x 10(-30) to 2420 x 10(-30) esu in acetonitrile, and they are among the highest reported for organometallic molecules so far. Electrochemical measurements are consistent with the push-pull nature of these compounds.

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

Massive black-hole binary inspirals: results from the LISA parameter estimation taskforce

The LISA Parameter Estimation Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.

Family violence in culturally and linguistically diverse communities: An evaluation of a family relationship centre

Federal legislative changes in Australia have sought to improve how family relationship centres (FRCs) can be more responsive to culturally and linguistically diverse (CALD) communities when addressing family disputes and family violence. Research on the prevalence of family violence against women from CALD backgrounds is sparse. This article seeks to contribute to the understanding of this issue by describing findings of an evaluation of the FRC at Broadmeadows conducted by the authors. The findings focus on family violence and the appropriateness of the services offered to three CALD community groups, namely Iraqi, Lebanese and Turkish. The final section provides reflections on these findings and offers suggestions about how FRC services might better cater to the family violence needs of CALD communities by developing more extensive partnerships with community groups and by expanding the range of processes they offer beyond mediation.

Designs for phase I clinical trials with multiple courses of subjects at different doses

The goal of this article is to provide a new design framework and its corresponding estimation for phase I trials. Existing phase I designs assign each subject to one dose level based on responses from previous subjects. Yet it is possible that subjects with neither toxicity nor efficacy responses can be treated at higher dose levels, and their subsequent responses to higher doses will provide more information. In addition, for some trials, it might be possible to obtain multiple responses (repeated measures) from a subject at different dose levels. In this article, a nonparametric estimation method is developed for such studies. We also explore how the designs of multiple doses per subject can be implemented to improve design efficiency. The gain of efficiency from "single dose per subject" to "multiple doses per subject" is evaluated for several scenarios. Our numerical study shows that using "multiple doses per subject" and the proposed estimation method together increases the efficiency substantially.

Estimating equations for removal data analysis

We consider the problem of estimating a population size from successive catches taken during a removal experiment and propose two estimating functions approaches, the traditional quasi-likelihood (TQL) approach for dependent observations and the conditional quasi-likelihood (CQL) approach using the conditional mean and conditional variance of the catch given previous catches. Asymptotic covariance of the estimates and the relationship between the two methods are derived. Simulation results and application to the catch data from smallmouth bass show that the proposed estimating functions perform better than other existing methods, especially in the presence of overdispersion.

The impact of global positioning systems and plotters on fishing power in the northern prawn fishery, Australia

The impact of global positioning systems (GPS) and plotter systems on the relative fishing power of the northern prawn fishery fleet on tiger prawns (Penaeus esculentus Haswell, 1879, and P. semisulcatus de Haan, 1850) was investigated from commercial catch data. A generalized linear model was used to account for differences in fishing power between boats and changes in prawn abundance. It was found that boats that used a GPS alone had 4% greater fishing power than boats without a CPS. The addition of a plotter raised the power by 7% over boats without the equipment. For each year between the first to third that a fisher has been working with plotters, there is an additional 2 or 3% increase. It appears that when all boats have a GPS and plotter for at least 3 years, the fishing power of the fleet will increase by 12%. Management controls have reduced the efficiency of each boat and lowered the number of days available to fish, but this may not have been sufficient to counteract the increases. Further limits will be needed to maintain the desired levels of mortality.

Fishery biology and management of Protonibea diacanthus (Sciaenidae) aggregations in far Northern Cape York Peninsula waters

The sciaenid Protonibea diacanthus is a large, long-lived predatory fish of inshore northern Australian waters, which forms annual aggregations that are fished extensively by traditional (subsistence) and recreational fishers. There are now widespread concerns that the resource is being overexploited. Indigenous fishers of the Cape York Northern Peninsula Area (NPA) relate that large adult fish (up to 1500 mm total length (TL)) made up the bulk of the catch from the sciaenid aggregations until about 1994. In contrast, sexually mature P. diacanthus comprised only a small component (12 fish out of 270=4.4%) examined in a 1999–2000 sampling programme that was biased towards the largest individuals available. At 790 mm TL, the minimum size at first maturity for female P. diacanthus in this study is much smaller than the 920 mm TL reported previously in Queensland waters. Developing ovaries were observed in specimens sampled from sciaenid aggregations which formed in NPA waters between May and September 2000. However, no fish with ripe or spent gonads were found in the study, so the current timing and location of the spawning season for P. diacanthus in the region remain unknown. Food items observed in the analysis of the diet of P. diacanthus from the NPA included a variety of teleosts and invertebrates. The range of animal taxa represented in the prey items support the description of an ‘opportunistic predator’ attributed to the species. In our sampling, the stomach contents of fish caught during the time of the aggregation events did not differ from those observed at other times of the year. A total of 114 P. diacanthus were tagged and released at aggregation sites during the study period, and 3 fish (2.6%) were subsequently recaptured. The low rate of tag returns from the wild stock tagging programme, both in this study (2.6%) and from recreational fisher tag/release programmes for the sciaenid elsewhere in Queensland (6.5%), were not explained by tag loss nor mortality, given the high retention rate of tags and the zero mortality seen in tank trials. In response to the biological findings from this study, indigenous community councils of the NPA imposed a 2-year fishing moratorium for P. diacanthus. Surveys at aggregation sites in 2002 and 2003 established that much larger fish (mean size 103.5 cm TL) were again present on the grounds, albeit in very low numbers. These recent preliminary results highlight the critical need for continued monitoring and management of the P. diacanthus fishery in the NPA, if prospects for resource recovery are to be realised. The NPA initiative has provided a rare opportunity to negotiate a co-management strategy, based on scientific data and traditional knowledge, for the recovery of a cultural and economically significant fished resource.

Rapid Global Expansion of the Fungal Disease Chytridiomycosis into Declining and Healthy Amphibian Populations.

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis, is enigmatic because it occurs globally in both declining and apparently healthy (non-declining) amphibian populations. This distribution has fueled debate concerning whether, in sites where it has recently been found, the pathogen was introduced or is endemic. In this study, we addressed the molecular population genetics of a global collection of fungal strains from both declining and healthy amphibian populations using DNA sequence variation from 17 nuclear loci and a large fragment from the mitochondrial genome. We found a low rate of DNA polymorphism, with only two sequence alleles detected at each locus, but a high diversity of diploid genotypes. Half of the loci displayed an excess of heterozygous genotypes, consistent with a primarily clonal mode of reproduction. Despite the absence of obvious sex, genotypic diversity was high (44 unique genotypes out of 59 strains). We provide evidence that the observed genotypic variation can be generated by loss of heterozygosity through mitotic recombination. One strain isolated from a bullfrog possessed as much allelic diversity as the entire global sample, suggesting the current epidemic can be traced back to the outbreak of a single clonal lineage. These data are consistent with the current chytridiomycosis epidemic resulting from a novel pathogen undergoing a rapid and recent range expansion. The widespread occurrence of the same lineage in both healthy and declining populations suggests that the outcome of the disease is contingent on environmental factors and host resistance.

Great Barrier Reef Fish Biological Data

Biological measurements on fish sampled during the course of FRDC funded project Growth, Reproduction and Recruitment of Great Barrier Reef Food Fish Stocks (FRDC 90/18). The comma-delimited ascii file comprises the following fields: 1. Cruise number 2. Date (d-m-y) 3, Region (descriptor of part of Queensland coast or Great Barrier Reef system) 4. Reef (name or number) 5. Data source (Res=research, Rec=recreational fisher, Com=commercial fisher) 6. Capture method 7. Trap number (where appropriate) 8. Species name 9. LthStd (standard length, cm) 10. LthFrk (fork length, cm) 11. LthTot (total length, cm) 12. WtTot (approx total weight, g; weighed at sea) 13. FrameWt (weight of frame [after filleting, with viscera], g; weighed in lab) 14. Sex (macroscopic examination only) 15. GonadWt (g) Data obtained by the Department Employment, Economic Development and Innovation (formerly Primary Industries and Fisheries) between 1988 and 1993, primarily in the southern Great Barrier Reef (Capricorn-Bunker and Swain Groups), with fish traps and handlining.

Probabilistic load flow for distribution systems with uncertain PV generation

Large integration of solar Photo Voltaic (PV) in distribution network has resulted in over-voltage problems. Several control techniques are developed to address over-voltage problem using Deterministic Load Flow (DLF). However, intermittent characteristics of PV generation require Probabilistic Load Flow (PLF) to introduce variability in analysis that is ignored in DLF. The traditional PLF techniques are not suitable for distribution systems and suffer from several drawbacks such as computational burden (Monte Carlo, Conventional convolution), sensitive accuracy with the complexity of system (point estimation method), requirement of necessary linearization (multi-linear simulation) and convergence problem (Gram–Charlier expansion, Cornish Fisher expansion). In this research, Latin Hypercube Sampling with Cholesky Decomposition (LHS-CD) is used to quantify the over-voltage issues with and without the voltage control algorithm in the distribution network with active generation. LHS technique is verified with a test network and real system from an Australian distribution network service provider. Accuracy and computational burden of simulated results are also compared with Monte Carlo simulations.

Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection.

In order to investigate the effect of long term recurrent selection on the pattern of gene diversity, thirty randomly-selected individuals from the progenitors (p) and four selection cycles (C0, C3, C6 and C11) were sampled for DNA analysis from the tropical maize (Zea mays L.) breeding populations, Atherton 1 (AT1) and Atherton 2 (AT2). Fifteen polymorphic Simple Sequence Repeat markers amplified a total of 284 and 257 alleles in AT1 and AT2 populations, respectively. Reductions in the number of alleles were observed at advanced selection cycles. About 11 and 12% of the alleles in AT1 and AT2 populations respectively, were near to fixation. However, a higher number of alleles (37% in AT1 and 33% in AT2) were close to extinction. Fisher's exact test and analysis of molecular variance (AMOVA) showed significant population differentiations. Gene diversity estimates and AMOVA revealed increased genetic differentiations at the expense of loss of heterozygosity. Population differentiations were mainly due to fixation of complementary alleles at a locus in the two breeding populations. The estimates of effective population at an advanced selection cycle were close to the population size predicted by the breeding method.

Determination of management units for grey mackerel fisheries in northern Australia.

The requirement for Queensland, Northern Territory and Western Australian jurisdictions to ensure sustainable harvest of fish resources and their optimal use relies on robust information on the resource status. For grey mackerel (Scomberomorus semifasciatus) fisheries, each of these jurisdictions has their own management regime in their corresponding waters. The lack of information on stock structure of grey mackerel, however, means that the appropriate spatial scale of management is not known. As well, fishers require assurance of future sustainability to encourage investment and long-term involvement in a fishery that supplies lucrative overseas markets. These management and fisher-unfriendly circumstances must be viewed in the context of recent 3-fold increases in catches of grey mackerel along the Queensland east coast, combined with significant and increasing catches in other parts of the species' northern Australian range. Establishing the stock structure of grey mackerel would also immensely improve the relevance of resource assessments for fishery management of grey mackerel across northern Australia. This highlighted the urgent need for stock structure information for this species. The impetus for this project came from the strategic recommendations of the FRDC review by Ward and Rogers (2003), "Northern mackerel (Scombridae: Scomberomorus): current and future research needs" (Project No. 2002/096), which promoted the urgency for information on the stock structure of grey mackerel. In following these recommendations this project adopted a multi-technique and phased sampling approach as carried out by Buckworth et al (2007), who examined the stock structure of Spanish mackerel, Scomberomorus commerson, across northern Australia. The project objectives were to determine the stock structure of grey mackerel across their northern Australian range, and use this information to define management units and their appropriate spatial scales. We used multiple techniques concurrently to determine the stock structure of grey mackerel. These techniques were: genetic analyses (mitochondrial DNA and microsatellite DNA), otolith (ear bones) isotope ratios, parasite abundances, and growth parameters. The advantage of using this type of multi-technique approach was that each of the different methods is informative about the fish’s life history at different spatial and temporal scales. Genetics can inform about the evolutionary patterns as well as rates of mixing of fish from adjacent areas, while parasites and otolith microchemistry are directly influenced by the environment and so will inform about the patterns of movement during the fishes lifetime. Growth patterns are influenced by both genetic and environmental factors. Due to these differences the use of these techniques concurrently increases the likelihood of detecting different stocks where they exist. We adopted a phased sampling approach whereby sampling was carried out at broad spatial scales in the first year: east coast, eastern Gulf of Carpentaria (GoC), western GoC, and the NW Northern Territory (NW NT). By comparing the fish samples from each of these locations, and using each of the techniques, we tested the null hypothesis that grey mackerel were comprised of a single homogeneous population across northern Australia. Having rejected the null hypothesis we re-sampled the 1st year locations to test for temporal stability in stock structure, and to assess stock structure at finer spatial scales. This included increased spatial coverage on the east coast, the GoC, and WA. From genetic approaches we determined that there at least four genetic stocks of grey mackerel across northern Australia: WA, NW NT (Timor/Arafura), the GoC and the east Grey mackerel management units in northern Australia ix coast. All markers revealed concordant patterns showing WA and NW NT to be clearly divergent stocks. The mtDNA D-loop fragment appeared to have more power to resolve stock boundaries because it was able to show that the GoC and east coast QLD stocks were genetically differentiated. Patterns of stock structure on a finer scale, or where stock boundaries are located, were less clear. From otolith stable isotope analyses four major groups of S. semifasciatus were identified: WA, NT/GoC, northern east coast and central east coast. Differences in the isotopic composition of whole otoliths indicate that these groups must have spent their life history in different locations. The magnitude of the difference between the groups suggests a prolonged separation period at least equal to the fish’s life span. The parasite abundance analyses, although did not include samples from WA, suggest the existence of at least four stocks of grey mackerel in northern Australia: NW NT, the GoC, northern east coast and central east coast. Grey mackerel parasite fauna on the east coast suggests a separation somewhere between Townsville and Mackay. The NW NT region also appears to comprise a separate stock while within the GoC there exists a high degree of variability in parasite faunas among the regions sampled. This may be due to 1. natural variation within the GoC and there is one grey mackerel stock, or 2. the existence of multiple localised adult sub-stocks (metapopulations) within the GoC. Growth parameter comparisons were only possible from four major locations and identified the NW NT, the GoC, and the east coast as having different population growth characteristics. Through the use of multiple techniques, and by integrating the results from each, we were able to determine that there exist at least five stocks of grey mackerel across northern Australia, with some likelihood of additional stock structuring within the GoC. The major management units determined from this study therefore were Western Australia, NW Northern Territory (Timor/Arafura), the Gulf of Carpentaria, northern east Queensland coast and central east Queensland coast. The management implications of these results indicate the possible need for management of grey mackerel fisheries in Australia to be carried out on regional scales finer than are currently in place. In some regions the spatial scales of management might continue as is currently (e.g. WA), while in other regions, such as the GoC and the east coast, managers should at least monitor fisheries on a more local scale dictated by fishing effort and assess accordingly. Stock assessments should also consider the stock divisions identified, particularly on the east coast and for the GoC, and use life history parameters particular to each stock. We also emphasise that where we have not identified different stocks does not preclude the possibility of the occurrence of further stock division. Further, this study did not, nor did it set out to, assess the status of each of the stocks identified. This we identify as a high priority action for research and development of grey mackerel fisheries, as well as a management strategy evaluation that incorporates the conclusions of this work. Until such time that these priorities are addressed, management of grey mackerel fisheries should be cognisant of these uncertainties, particularly for the GoC and the Queensland east coast.