Aims: To improve knowledge of the epidemiology of melanoma by comparing occupational risks of cutaneous melanoma (CM) by anatomic site in Swedish workers.
Methods: Male workers employed in 1970 and living in the country in 1960 were followed up from 1971 to 1989 using the Swedish Registers of Death and Cancer. A more specifically exposed subcohort included men reporting the same occupation in 1960 and 1970. For each location, occupational risk ratios (RRs) were extracted from Poisson regression models adjusted by age, period, town size, and geographical area. To diminish the influence of socioeconomic factors, intrasector analyses, comparing only jobs belonging to the same occupational sector, were performed. Risk patterns for different locations were compared.
Results: High RRs for different sites were found among workers exposed to UV sources (dentists, physiotherapists, and lithographers), and sun exposed workers (harbour masters, and lighthouse/related work). Risk excesses were seen in fur tailors, tanners/fur dressers, patternmakers/cutters, electrical fitters/wiremen, telephone/telegraph installers/repairmen, and some glass/pottery/tile workers. Results for lower and upper limbs were significantly correlated but somewhat independent of those found in thorax, the most frequent location. Correlation between head/neck and thorax was moderate. Specific risk excesses were found for rolling mill workers in head/neck, for chimney sweeps in upper limbs, and for aircraft pilots/navigators/flight engineers in lower limbs.
Conclusions: High RRs in the trunk among occupations with UV exposure from artificial sources suggest an effect not restricted to exposed sites. An unusual distribution of cases and RRs in chimney sweeps, rolling-mill, or glass/pottery/tile workers suggests local effects of exposures. The not previously reported risk excess in this job and in fur related processes, and the RR in electrical fitters and telephone/telegraph installers deserve further investigation. Disparities between locations, as RRs in thorax and limbs, may reflect differences in aetiological mechanisms.
- CM, cutaneous melanoma
- RR, risk ratio
- UV, ultraviolet
- PCBs, polychlorinated biphenyls
Statistics from Altmetric.com
Although melanoma is one of the neoplasms with higher increases of incidence in the past 20 years, especially in fair skinned populations,1 its aetiology remains poorly understood. Familial cases represent only 4–15%,2 reflecting the fact that environmental factors are the most probable aetiological agents. Even among subjects with germline mutations in the CDKN2A gene, which are the most common cause of inherited susceptibility, the penetrance presents a geographical variation in agreement with rates observed in the general population.3 Thus, environmental factors might also have an important role in familial melanoma. Nowadays, sun exposure is the only generally acknowledged aetiological agent,4 complemented with the increased susceptibility that some cutaneous and pigmentary characteristics confer to sunburn,5 which have been recently associated with certain melanocortin-1 receptor polymorphisms.6 Sunlamps and sunbeds have also been related to the occurrence of this neoplasm.7,8
Differences in exposure or skin contact with harmful agents, including UV sources, might be linked, not only with the risk of the disease, but also with the anatomic site where it appears. Melanoma is more frequently found in covered parts of the body, not usually exposed to sunlight. Both sexes have a different site distribution: while the trunk is the most common location in men, legs are the primary site in women, and this different pattern is consistent in countries with large variation in latitude.9 Differences in trends by site have also been described.10–13 Some anatomic locations depict specific epidemiological features, such as a predominance of lentigo malignant melanoma and an older age at diagnosis in the head and neck.14 These differences have been related to patterns of sun exposure,15–17 or with a site dependent susceptibility of melanocytes.18 Others have suggested the existence of different pathways that could be related with anatomic site.19,20
Occupational settings present a large variation in environmental exposures by anatomic site; consequently, comparison of anatomic risk distributions might be a source of hypothesis about the aetiology of this neoplasm. This approach led some authors to suggest that the increased risk of melanoma in white collar workers, particularly in office jobs,21,22 and its preference for normally covered sites, could be due to the effect of intermittent sun exposure.4,23,24 Recently the Association of Swiss Cancer Registries25 has included melanoma by site in its global analysis of cancer per occupation. However, their data source is a pool of several registries with different occupational classifications and they use as reference all registered cases, which could underestimate effects in jobs highly exposed to carcinogens.
Our main objective is to study in a more comprehensive way the relation between occupation and melanoma in men by giving a comparative view of occupational risks of this neoplasm by sites within a whole country. The high rate of incidence of melanoma in Sweden26 and the availability of a large historical cohort27 make this country an ideal candidate for this analysis. This is a more homogeneous population, and the restriction imposed on all members of the cohort of being present in both the 1960 and 1970 censuses assures their residence in the country for at least 10 years.
Some occupations, such as electrical fitters, show a general risk excess of cutaneous melanoma, suggesting the existence of general promoting factors, while others have only site specific excesses, suggesting local effects of exposures.
If occupational exposure to artificial UV devices increases the risk of cutaneous melanoma, its effect is not limited to exposed sites.
Disparities between locations, such as RRs in thorax and limbs, may reflect differences in aetiological mechanisms.
The relation between total melanoma and occupation has already been studied in Sweden with a related cohort28 and with our material, although with a shorter follow up period29 or in the context of a general study relating occupation and all types of cancer.30,31 None of these has focused on site specific analysis, nor has taken into account some of the main confounders already described for this country,32 such as urban/rural distribution (probably a proxy of travel habits to sunny countries),33 and geographical distribution (reflecting environmental UV exposure). Socioeconomic class has not been considered in these previous studies, despite its well known relation with melanoma,22,34,35 which is thought to reflect lifestyle differences.
Our study presents site specific job title relative risks adjusted by town size and geographical distribution, and is based on comparisons within occupational sectors, to compare people with a more homogeneous socioeconomic status. We have also compared site risk distributions in order to determine any disparities or similarities between locations that could point to possible aetiological differences.
The base population for this historical cohort was made up of Swedish men gainfully employed at the time of the 1970 census, present in the 1960 census, and still alive and aged 25–64 years on 1 January 1971. The cohort included 1 890 497 men followed up until the end of 1989, and, inside it, a subcohort including only those with the same occupation in both the 1960 and 1970 censuses, comprising 755 728 men.
Information was drawn from two linked data sets: (1) the Swedish Cancer Environment register, providing information on incident cases and including a number of demographic variables from the 1960 and 1970 censuses, which was used to compute specific rate numerators; and (2) a background population register comprising all individuals in the 1970 census, with information on occupation and residence in 1970, occupation in 1960, and, if applicable, date of death. This register was used to calculate specific rate denominators. The record linkage between these two registers has been described in detail elsewhere.36 Melanoma was coded 190 in the International Classification of Diseases (7th revision). The fourth digit specifies anatomic site. All head and neck melanomas have been analysed together. It must be taken into account that cutaneous melanomas with multiple sites or non-specified location are included in the global melanoma data but not in the site analysis. They made up 0.7% and 10.6% of all registered melanomas respectively.
The not previously reported excess risk found in fur related processes, and in glass/pottery/tile works, should be further investigated
Person-years in each of 278 occupations were accumulated from 1971 until the date of death or up to the end of 1989. The overall person-time that each worker contributed to the study was allocated to the corresponding cells of the variables of stratification. These were occupation, county and size of town of residence in 1970, five year age classes, and the following calendar periods: 1971–75, 1976–80, 1981–85, and 1986–89. Occupation, town size, and county, drawn from the 1970 census, were regarded as fixed variables, while age and period were time dependent. Clayton’s algorithm37 was used in calculating the exact number of person-years.
Occupation was classified according to the Nordic Classification of Occupations.36 Every occupation is represented by a three digit number. The first digit refers to one of 10 major occupational sectors (0–9), where higher numbers indicate manual occupations and lower numbers non-manual occupations, often requiring longer education and associated with a higher socioeconomic status. Expected cases were calculated for each occupation using the age-period specific rates for the whole cohort as referent. They were initially used to calculate age-period standardised incidence ratios (SIR).
On the assumption that the observed number of cases was distributed in each stratum as a Poisson variable, log-linear Poisson models were fitted to estimate risk ratios (RRs). First, the classification of Swedish counties into low, medium, or high risk for melanoma was carried out using a Poisson model adjusted by sex, age, period, occupational sector, and town size, and taking as reference the incidence for the whole country. Counties were grouped into three intervals according with their RRs: <0.85, 0.85–1.15, or >1.15 (see fig 1).
For each sector and occupation, RRs by site, adjusted by town size and geographical area, were calculated using log-linear Poisson models. Town size was included as a confounder since it showed a positive correlation with melanoma risk. In these models, the number of expected cases was introduced as an offset.37 As the expected numbers were computed based on the age and period specific reference rates, RRs were likewise age and period adjusted.
Occupations with at least five observed cases of total cutaneous melanoma were analysed. Each occupation was compared only with others within the same occupational sector. In general, RRs for different sectors in the cohort and subcohort showed a socioeconomic gradation, with risk excesses in sectors 0–III (white collar workers), reduced risks in production sectors, mining, and agriculture, and an intermediate position for sectors VI (transport & communication) and IX (services), though in head/neck this gradation was not so evident. For this reason, we preferred this intrasectorial approach in order to contrast people with a more homogeneous socioeconomic status, though it might underestimate risks if there is any relevant exposure shared by workers in the whole sector (that is, sun exposure in agriculture). Production sectors (sectors VII and VIII) were treated as a single category.
These analyses were also repeated for all cases with the subcohort reporting the same occupation in the 1970 and in 1960 censuses, a more specifically exposed group which has been used mainly to check the consistency of the results found in the general cohort; however, the uneven reduction of the cohort that this restriction implies has several consequences: (1) a lower number of cases by job title, which reduces the statistical power or even makes it impossible to compute the risk in some jobs; (2) the absence of people in jobs which did not exist or were uncommon in men in 1960 (that is, systems analysts and programmers, or nurses), or in occupations with a lower stability (that is, those implying a promotion, such as staff officers); and (3) a slightly different composition of sectors between the subcohort and the cohort.
Logarithms of occupational RRs obtained with the whole cohort were subsequently used to asses similarities and differences in the risk pattern across different locations, using for each comparison all occupations with two or more cases in both locations considered. Results are presented graphically, highlighting occupations with high RRs and those with discrepant RRs for different sites. Pairwise site agreement was tested using Spearman’s correlation coefficient.
During follow up, 6187 cases of cutaneous melanoma were found in the whole cohort and 2816 in the subcohort. Thorax and upper limbs accounted for 51% and 12% of the cases respectively in both. Head and neck melanomas were slightly more frequent in the subcohort (12% v 14%) and lower limbs in the general cohort (11% v 9%).
Tables 1 and 2 show age, period, town size and geographically adjusted intrasectorial RR by job title for all melanomas in the cohort and subcohort, and by sites in the cohort respectively. All occupations with a minimum of five cases in all melanomas and (a) RR ⩾1.5 for global cutaneous melanoma, or (b) a minimum of three cases and RR ⩾2 in any site, regardless of its statistical significance, are reported. For comparison purposes, all job titles fulfilling the established criteria were included in both tables. However, those jobs with less than three cases in all the specific sites are only reported in table 1. A higher cut-off has been used for specific sites as higher risks were also expected given the lower number of cases. Due to the complementary character of the information that cohort and subcohort analyses provide, the results of both approaches will be discussed together.
Risks in job titles not included in any of the tables are available, as well as RR using as reference the whole male working population, and can be provided on request. As expected, in this last analysis, a substantial proportion of occupations with increased RR belonged to sectors 0–III.
All cutaneous melanoma (table 1)
In general, results in the whole cohort and in the 60–70 subcohort were quite consistent. Significant increases of risks over 50% in the cohort (13 jobs) and/or in the subcohort (13 jobs) were found in 20 job titles if all cases are considered. In sectors 0–III it should be noted the risk excess found in certain health professions (dentists and physiotherapists/occupational therapists). Librarians/archivists and curators, composers/musicians, and insurance agents had significant RRs only in the subcohort analysis, while an excess of risk in university teachers was found in the cohort analysis alone. Well driller was the only job with increased RR in sector V. In several jobs comprising harbour workers, mostly from sector VI, there was also a significant RR (forwarding/shipping agents, harbour masters, and lighthouse/lock operators/ferry and harbour service assistants), while post office clerks had a consistent but lower RR. Of special interest within production sectors (VII and VIII) can be considered some occupations within textile industry, such as fur tailors and tanners/fur dressers, patternmakers/cutters, and leather good makers. It is also noteworthy that electrical fitters and wiremen, and telephone-telegraph installers and repairmen depicted a significant increased incidence. Another interesting result was the risk excess within the glass/pottery/tile industry, with an important, consistent, and significant RR in non-specified workers, and high and nearly significant RRs for glass formers/cutters. Additionally, truck/conveyor operators and typographers/lithographers had consistent high risks in both the cohort and the subcohort. In sector IX, catering supervisors had significant RRs only in the subcohort analysis, while risk excess among prison/reformatory officials was found in the whole cohort.
Head and neck (table 2)
A significant RR >2 was found in 10 job titles in the cohort analysis. Despite socioeconomic adjustment, most of the high risk jobs belonged to sectors 0–III (chemists, university teachers, ministers/priests, composers/musicians, social workers, bank employees, forwarding/shipping agents, and commercial travellers). Within blue-collar workers, only patternmakers and cutters, rolling-mill workers, non-specified glass/pottery/tile workers, and truck/conveyor operators had a statistically significant RR.
Thorax (table 2)
In six job titles, there was a significant RR >2 in the cohort, and in other three occupations a significant but lower RR was also found. Six of these nine jobs belonged to sectors V–VIII, thus showing a predominance of blue-collar occupations in this site after intrasectorial adjustment.
Veterinarians had a marked increase of risk, and a significant risk excess was also found in dentists, well drillers, harbour masters, lighthouse/lock operators/harbour service assistants, fur tailors, electrical fitters/wiremen, typographers/lithographers, and prison/reformatory officials.
Upper limbs (table 2)
In six job titles there was a statistically significant RR >2 in the cohort (teachers in theoretical subjects, lawyers in private practices, insurance agents, other brick/concrete work, typographers/lithographers, and chimney sweeps). Electrical fitters and wiremen had a slightly lower excess of risk, which also achieved statistical significance.
Lower limbs (table 2)
Only four job titles had significant increases of risk in this location, and all of them had very high RRs: forestry manager and supervisors, aircraft pilots/navigators/flight engineers, pattern makers/cutters, and non-specified glass/pottery/tile work. Almost significant high RRs were also found in precision toolmakers, packers, and launderers/dry cleaners.
Estimation of error due to chance
Table 3 shows an estimation of expected significant associations with RR >1 due to chance by site and analysis, assuming two tailed α = 0.05. We also report the number of observed significant associations and those over the established threshold (1.5 for all cases and 2 for the sites). In general, between a third and a quarter of the statistically significant results might be due to chance. The use of a cut-off point should increase specificity of our results, even though RR also critically depends on the number of exposed workers. The importance of our findings should, however, be judged in the light of the available knowledge.
Correlations of results by sites
Figure 2 shows pairwise comparisons between occupational log(RR) per site. In each graph, those occupations with discordant risks (logRR <percentile 15 in one site and logRR >percentile 85 in the other) as well as jobs with increased risk (logRR >percentile 85 for both locations) are signalled. The absence of significant correlation of thorax with either lower or upper limbs is noteworthy. A significant correlation of 0.48 was found between upper and lower limbs, which increased to 0.54 when the analysis was restricted to occupations fulfilling the criteria used for inclusion in table 2. Head/neck had low but significant correlation with thorax and lower limbs, but not with upper limbs.
For many authors, a progressive change from predominantly occupational to predominantly recreational patterns of sun exposure is the most likely cause of the increasing incidence of melanoma in Caucasian populations.9 However, and in accordance with our results, some workers are at increased risk of cutaneous melanoma, and the occupational pattern is not exactly the same in different sites. Holiday travel and leisure time activities cannot account for all the variability of risks found between occupations and their distribution per site. Westerdahl and colleagues33 did not find differences in travel habits to south Europe between Swedish patients with different site locations. Furthermore, the use as reference group of people working in its own professional sector for each job title should at least partially control confounding due to those factors, which are highly related with socioeducational status.
Population occupational studies such as ours imply multiple comparisons, which might produce a certain number of spurious significant associations, the so-called mass significance phenomenon. We have tried to quantify this source of error, and in table 3 we provide an estimation of expected associations due to chance by site and analysis, together with the real number of associations obtained. However, we have also taken into account several considerations to weight the plausibility of the associations found. First, consistency between cohort and subcohort analysis results has been valued, because, though highly related, both analyses have been done separately. Second, more credit has been given to high RRs found in jobs with quite similar occupational exposures. Finally, those results already highlighted in the literature are specifically discussed.
For all cutaneous melanoma, our results show, in the first place, a significant risk excess in the cohort and subcohort analyses for all cases in dentists. Higher incidences than expected by chance for dentistry have been found in pooled SIR of Nordic countries and individually in Norway, Sweden, and Denmark,30 Switzerland,25 England,38 and the USA.39 Typographers/lithographers had also consistent and significant excess. Printers have been associated with melanoma in the combined analysis of incidence in four northern European countries—including Sweden.30 Swiss cancer registers25 described a significant OR of 1.6 for all printers, and 1.7 for the subcategory including lithographers. In a more specific way, Nielsen and colleagues40 found a significant RR of 3.4 in their study of incidence in a cohort of Danish lithographers. Increases in melanoma mortality in printers have been reported in the USA and France.41,42 It could be hypothesised that these associations might be attributed to the use of artificial sources of UV radiation in these two jobs,43 which could also account for the high RR for physiotherapists in the subcohort. However, this possible explanation cannot justify the negative results in other jobs with known exposure to UV radiation,43 such as welders. Site distribution of risk does not clarify this hypothesis. For sun exposure, Chen and colleagues44 found that a history of sunburn at an anatomic site was specifically related to the development of malignant melanoma at that particular site. The most commonly accepted intermittent sun exposure theory also points towards this site specificity, more than to a general “solar circulating factor” which would increase global melanoma risk.45 If an analogous effect were expected for other UV sources, it must be noted that workers exposed are usually dressed; thus exposure occurs mainly in head and neck, or upper limbs if short sleeves are worn. Our data show a non-significant high risk in the head and neck for dentists and a significant excess in upper limbs in lithographers/typographers, but in both jobs there is also an unpredicted significant increase in the thorax. For Swiss printers including lithographers,25 a significant OR for lower limbs, and a high, though non-significant OR in the head/neck and thorax has been reported, while four out of five of Danish lithographers’ melanomas were situated in the trunk. However, as the relation between sun and melanoma is not yet clarified, it is possible that UV effects from solar or non-solar sources could be both local and general.
Harbour masters, and lighthouse and related workers, who could be intermittently sun exposed, had a significant global excess of risk, with high RR in the thorax and head/neck. Similar results were observed in truck/conveyor operators.
Of special interest are the non-previously reported increased risks in fur tailors (both in the cohort and in the subcohort for all cutaneous melanomas, as well as in thorax, and upper limbs), and in tanners and fur dressers. Fur tailors are exposed to trichloroethylene, which was related to an excess of melanoma in a Canadian case-control study;46 an increased incidence was also found in a Californian community where drinking water was contaminated with ammonium perchlorate and trichloroethylene.47 However, the fact that both steps in the manufacture of furs have significant excesses of incidence points to common substances used by both jobs, and rests of dyes or chemicals used by fur dressers, still present in the fur at its tailoring process, could be mobilised with the high temperatures of ironing. We also found a higher risk of melanoma in patternmakers/cutters, with a significant RR >2 for all cases, head/neck, and lower limbs.
It is interesting to point out that electrical fitters/wiremen had a consistent increased incidence in each site; though only in the thorax and upper limb does it achieve statistical significance. High risks in electrical workers, with48 or without socioeconomic adjustments,46,48–50 have been reported. This risk excess has been related to polychlorinated biphenyl (PCB) exposure through skin contact, and a dose-response effect has been described,51 although others have failed to find a clear trend. The generalised effect in all sites found in our data could point towards an exposure affecting the whole body, for example, cloth contamination by PCBs, but other explanations have also been proposed. Telephone/telegraph installers and repairman presented a high global RR, with a non-significant small increase in thorax. Some authors have found a relation between electromagnetic fields (EMF) and melanoma,52,53 which could explain the results in these latter occupations. However, there is not yet enough scientific evidence to support this association.54
A very high RR in non-specified glass/pottery/tile workers has been found in global data, head/neck, and lower limbs. This job title might include crystal glassmakers and some glass blowers, and decorating glass and artistic glaziers. Interestingly glass formers/cutters had nearly significant high risks if all melanomas are considered. Furthermore, another glass industry occupation (glass/ceramic kilnmen), not included in the tables as it only had four cases of melanoma, also showed non-significant risk excess for global melanoma, accompanied by a significant RR of 3.18 in the thorax. These results have not been previously reported, not even in a former study of mortality in the Swedish glass industry.55 Their work environment is highly complex, and workers come into close contact with a variety of potentially hazardous agents, including metals such as arsenic and lead in art glass.56
The head/neck is generally considered to be a chronically sun exposed site, and presents a higher proportion of lentigo malignant melanoma,14,57 whose most widely accepted risk factor is long term cumulative UV radiation.58 Thus, higher risks in this site could be expected for outdoor workers; in fact, a risk excess in such occupations as construction workers,59 or farmers25 has been reported. Our data, in contrast, show a clear majority of indoor jobs, mostly belonging to high socioeconomic sectors. We do not have a specific hypothesis to explain the high RR in university teachers, priests, composers/musicians, or social workers. Intrasectorial variation in the frequency of travel to foreign countries would perhaps justify these results; however, this would not explain why they failed to show clearer risk excess for other locations. High risk of all cutaneous melanoma has been previously reported in teachers and priests or other religious workers in Norway and Finland30 or in Switzerland,25 but socioeconomically adjusted estimators, only provided in the last study, did not confirm this excess. Several studies have suggested a higher risk of melanoma for chemical workers,30,60 though in our study it was mainly limited to head and neck. Finally, raised risk of melanoma in rolling mill workers has been previously found.61
Chimney sweeps had increased incidence in upper limbs, and a high, though not significant RR in lower limbs. In this occupation five of six cases are not located in the thorax, which might suggest skin contact with a harmful agent. The Swiss study25 did not find any melanoma cases in this occupation.
In lower limbs, it is noteworthy that a high RR is found in aircraft pilots/navigators/flight engineers. Flight crews have been repeatedly associated with an excess of melanoma, and usually with high risk estimators.38,62–64 A recent meta-analysis reported a combined, socioeconomically adjusted RR among male pilots for melanoma mortality of 1.97.65 Despite this, no estimation of risk by site has been published, and usually information about anatomic location is scarce. This risk excess has been mainly attributed to sunbathing in sunny destinies, as UV exposure on the flight deck is negligible.66 Other alternative occupational factors, such as melatonin disturbances or radiation dose have been proposed, since risk increases with hours of flight.62,63 However, if their risk was mainly focused on legs, it could be related to clothing habits, as they might be more prone to wear short trousers than other workers, favouring intermittent sun exposure. Forestry managers/supervisors also had a non-previously reported significant RR of 3 in the lower limbs. The use of the global cohort as a reference did not substantially modify these results, signalling that they are not probably due to a privileged socioeconomic status within their sector.
Being the most frequent location, results in the thorax mostly mirrored those found for all sites combined, with some exceptions, such as the significant RR observed in veterinarians. As in our study, Goodman and colleagues39 found a significant risk excess in this site for indoor/outdoor and outdoor workers after adjusting for training level and birth place. In general, it is difficult to think of local occupational exposures occurring preferentially in the thorax, which is usually a cloth covered site, but those factors with a global promoting effect for melanoma should also show their effect at this location, which represents 32% of the whole skin surface.67
Correlation coefficients showed that, globally, the pattern of risk was different in thorax and limbs, while upper and lower limbs results were significantly correlated. This curious finding could perhaps be explained by more similar patterns of occupational—or non-occupational—exposure to agents increasing risk of melanoma, including UV radiation. An alternative explanation could be the site specific differences in susceptibility to melanoma18 or the existence of two site related different pathways for this neoplasm19,20 that some authors have proposed, with the skin of both limbs in men being considered more alike than that of the trunk.
Of interest is the fact that the head and neck failed to correlate with the upper limbs, but showed a significant small correlation with the thorax and lower limbs.
In summary, a high incidence among dentists and printers/lithographers in the trunk only supports the aetiological role of occupational artificial sources of UV radiation if their effect is presumed not to be restricted to exposed sites. High RRs in certain locations and the unusual site distribution of cases found in occupations such as chimney sweeps, or rolling-mill or non-specified glass/pottery/tile workers could hint towards agents with a local effect. A generalised increase of risk in electrical fitters and high RR in telephone/telegraph installers, as well as the not previously reported risk excess of melanoma in glass/tile/pottery workers and in fur related processes requires further investigation. Regarding site correlation, different patterns might be due to an uneven distribution of exposures or to site differences in susceptibility.
Dr Perez-Gomez is indebted to Dr Dario Garcia de Viedma for his invaluable help and support.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.