Low presence of potentially toxic elements in Singapore urban garden soils

Urban agriculture is potentially an important piece of the food security puzzle for a rapidly growing urban world population. Community gardening is also promoted as a safe and viable form of exercise for aging populations in crowded settings where opportunities to participate in other action activities may be limited. Knowledge of potential site-specific health risks to environmental contaminants is important in dialogues promoting urban farming. We assess the pseudo-total concentrations of selected potentially toxic elements (PTEs) in the soils of community gardens, public parks, and woodlands in the tropical urban island nation of Singapore. We compare concentrations of cadmium, copper, lead, and zinc with amalgamated risk guidelines to form a baseline understanding of the level of contamination in these spaces. We also perform providence tracking with lead isotopes to identify potential sources of contaminants. All pseudo-total concentrations of Cd, Cu, Pb, and Zn in the soil were below threshold concentrations considered to represent substantial risk. Further, PTE concentrations in gardens were largely equivalent to those found in community parks and woodlands, but the geographical distribution varied. Provenance tracking with Pb isotopes indicated Pb in gardens was both anthropogenic and natural, but spatially variable. The lack of strong spatial clustering of areas with the highest PTE concentrations was inconsistent with a common point source of contamination. However, the correlation between Cu and Zn suggest a common source for these elements, such as road/trafficking or atmospheric deposition. We find limited risk of urban gardeners to exposure to Cd, Cu, Pb, and Zn—elements that are commonly abundant in urban settings with dense transportation networks and substantial industrial activities. The low levels of PTEs are encouraging for the promotion of urban farming for food production and leisure in this dense urban setting. However, as concentrations were low, we did not assess bioavailability and bioaccessibility of the PTEs. These assessments would need to be determined in cases of with higher levels of contamination to provide a more thorough consideration of actual human risk.

Page 2 of 29 Goh et al. CABI Agriculture and Bioscience (2022) 3:60 (Mansfield and Mendes 2013;McClintock et al. 2018): "countries should aim to work to enhance inclusive and sustainable urbanization for participatory, integrated, and sustainable human settlement planning and management". Nearly a decade ago, Orsini et al. (2013) recognized that the food security benefit of urban agriculture was evidenced by 100-200 million people providing fresh horticultural goods to city markets. Arguably urban agriculture in developing countries plays a critical role in food availability and affordability, but there is need for improved understanding potential risks and limitations (De Bon et al. 2010;Zezza and Tasciotti 2010;Hamilton et al. 2014). While urban agriculture may not achieve food self-sufficiency per se, it is an alternative way of "feeding citizens differently" (Martin-Moreau and Ménascé 2020). Urban gardening is also a practical means of exercise, leisure, therapy, and perhaps supplementary income, in crowded urban spaces for all citizens (Orsini et al. 2013;Kalantari et al. 2018;Rogge et al. 2018;Winkler et al. 2019). Regardless of economic setting, it is increasingly important to preform foundational work addressing the potential realm of risks that await the growing wave of urban farmers, many of which may be advanced in age, and therefore in some cases, susceptible to the adverse effects of contaminated urban soils (Beckie and Bogdan 2010;Hamilton et al. 2014;Yang and Na 2017). Pollution research worldwide has revealed the potential for people engaged in outdoor urban activities to be exposed to various PTEs that originate from past and present anthropogenic sources such as transportation networks, energy production systems, mining and industrial activities, and improper handling of waste (McLaughlin et al. 2000;Mitchell et al. 2014;Jean-Soro et al. 2015;Spliethoff et al. 2016;Bechet et al. 2018;López et al. 2019). Elevated PTE concentrations in urban soils are often legacies of past activities conducted before a garden was established (Kessler 2013;Wortman and Lovell 2013). While some elements such as copper and zinc are essential in very small amounts for normal health, they may be toxic when high bioavailable fractions are encountered (De Vries et al. 2007).
Exposure to high concentrations of PTEs can cause a variety of health problems including dizziness, nausea, cramping, coughing, diarrhea, and vomiting (Jaishankar et al. 2014;Gautam et al. 2016). Long-term and repeated ingestion may induce element-specific maladies such as hypertension and weakening of bone density (Cd), renal dysfunction (Pb, Cd), anemia (Zn), and other afflictions (ATSDR 2015; Satarug and Moore 2004;Swaddiwudhipong et al. 2011;Cabral et al. 2015). The potential for adverse health risks stemming from a seemingly healthy and sustainable activity is one aspect of the "dark side of urban gardening" (Kessler 2011;Mancebo 2016)-however, others conclude the risk is sometimes more perceived than real (Brown et al. 2016;Lupolt et al. 2021). Nevertheless, assessing the potential risk is beneficial to promoting healthy outdoor activities.

Human urban gardening exposure chain
Community risk to exposure to potentially toxic elements (PTEs) while working in urban gardens has been a topic of concern for more than half a century (Chaney et al. 1984). Lead in the USA during the 1960s and 1970s, for example, was found in very high concentrations owing to contamination from a variety of sources including automotive exhaust, lead-based paint degradation, smelter emissions, sewage sludge, pesticides, and element-enriched manure fertilizers (Preer and Rosen 1977;Spittler and Feder 1979). Similarly, Purves (1966) found urban garden soils in Edinburgh and Dundee to be elevated in copper and boron, probably as a result of atmospheric deposition of soot, which was also applied as a garden soil amendment. Warren et al. (1971) later reported elevated concentrations of Cd, Cu, Pb, and Zn to be a growing concern for vegetable gardening in urban and industrial areas of Canada. Over the years, dozens of studies have examined the levels of these and other potentially harmful elements in urban and peri-urban soils worldwide (Table 1). The high element concentrations in the soils of many of the latter studies was often linked to legacy activities (e.g., Pb from paint; Cu, Pb, Zn from smelters and mines), trafficking phenomena (e.g., vehicle emissions, vehicle and tire wear), deposition following the release from a variety of urban and industrial sources, irrigation with wastewater, inter alia (Table 1; Additional file 1: Table S1).
These prior works have demonstrated that the "human urban exposure chain" to PTEs in the soil is a complex issue involving behavioral, biochemical, developmental, geophysical, physiological, and societal variables (cf. Clark et al. 2006;Schram-Bijkerk et al. 2018). To frame our research, which is a baseline assessment of Cd, Cu, Pb, and Zn in the soils of urban gardens in Singapore, we portray this exposure chain as a complex dynamic system. In the following narrative, the italicized terms are key variables in a simplified causal loop diagram shown in Fig. 1.
Element concentrations in the soil are influenced by several factors: mineralogy and geochemical composition of the parent rock, organic element content, particle size distribution, soil horizonation, drainage, vegetation, management, and anthropogenic inputs (Kabata-Pendias 1993; Wuana and Okieiman 2011). The elevated presence of a PTE in an urban soil could be naturally part of the background geogenic signature or associated with     enhancements caused by natural processes/sources such as volcanic eruptions, groundwater seepage, and fluvial transport (Briffa et al. 2020;Müller et al. 2020). Often, however, enrichment occurs when contaminants are transported to and deposited at a site from anthropogenic sources by atmospheric (related to wind and dry deposition) and hydrologic (e.g., rainfall, surface runoff, waste water irrigation) dispersal pathways, as well as the degradation of discarded waste or applied chemicals, including compost, fertilizers and disease/pest control agents that are applied in garden management (Kelly et al. 1996;Holt 2000;Wuana and Okiemen 2011;Su et al. 2022). While contamination may be acute, being associated with only a few pollution events (e.g., accidental spill/ release; flooding), it is often chronic, stemming from an accumulation over longer periods time (He et al. 2013;Ashraf et al. 2014). Thus, time or age are critical factors. For example, the soils associated with old homes/ neighborhoods often have elevated concentrations of some elements (inner city legacy), for example Pb from paint used to protect wooden and metal construction materials (Culbard et al. 1988;Li et al. 2001;Mitchell et al. 2014;Clarke et al. 2015;Rouillon et al. 2017;Kandic et al. 2019;Montaño-López and Biswas 2021). Gardens built on the former sites of long-running smelters, mines, factories, and waste dumps often are enriched in selected elements, representing a deposition legacy long after disuse (Kachenko and Singh 2006;Nabulo et al. 2012;Pelfrêne et al. 2012;Douay et al. 2013). Garden soils may become enriched over time from chronic atmospheric deposition (Kargar et al. 2013;Laidlaw et al. 2018). Thus, a new fill soil may not reflect the long-term pattern of contamination. Long-term application of metal-based amendments can also elevate PTEs in garden soils.
Gardens that are close in distance to mines, power generation plants, incineration facilities, factories, and manufacturing centers tend to receive higher pollution loads, but long-range dispersal is possible (Chaney et al. 1988;Marx and McGowan 2010;Hall et al. 2021). Gardens located in low-income neighborhoods are often in close distance to polluting activities (Mielke et al. 1983;Kaiser et al. 2015;Oguntade et al. 2020;Montaño-López and Biswas 2021). Roads and various phenomena related to trafficking (density/usage; road, vehicle and tire wear; emissions, etc.) are documented contributors to the PTEs found in roadside gardens (Tam et al. 1987;Hursthouse et al. 2004;Morton-Bermea et al. 2009;Folens et al. 2017). Airports, trains, and shipping are also potential sources of PTEs in urban gardens (Vazhacharickal et al. 2013;Cruz et al. 2014;Laidlaw et al. 2018;Leitão et al. 2018). As shown in the Table 1, other common sources of elevated PTE presence in garden soils are the application of contaminated sludge or irrigation water (Jamali et al. 2007;Abdu et al. 2011;Xue et al. 2012;Vazhacharickal et al. 2013) and element-based garden agrochemicals (Waterlot et al. 2011;Tran et al. 2020;Zhong et al. 2022).
Other types of interesting sources of PTEs reported in Table 1 include plastics ), E-waste processes (Luo et al. 2011), and war activities (Safi and Buerkert 2011).
Children are generally considered of greatest risk from exposure to particular PTEs because of two issues: (a) their state of health, which is generally characterized as having an undeveloped immune system; and (b) their risky behavior, which involves hand-to-mouth and object-to-mouth actions (Zeng et al. 2016;Różański et al. 2021). Elderly gardeners may also be at risk if their state of health has been compromised or they have experienced long-term exposure to particular elements (Sun and Gu 2008;Alghamdi et al. 2022). Nevertheless, the risk can apply to anyone working at contaminated sites who experiences prolonged dermal exposure (increased by lack of personal protection equipment), ingests high concentrations of elements through inhalation of dust during intense work or through hand-to-mouth transfer, or by consuming contaminated plants (Zahran et al. 2013;Antisari et al. 2015;Hough 2007;Tong et al. 2020).
Important is that the degree of toxicity experienced from exposure to PTEs in the contaminated media is not simply related to the total concentration of the element involved, but importantly, the bioavailability and bioaccessability of the element (Chapman et al. 2003;Janssen et al. 2003). Bioavailability refers to the fraction of PTEs that is released from a mineral matrix into a liquid form that can be taken up by plants or released into the gastrointestinal tract upon digestion Petruzzelli et al. 2020); and bioaccessability is the proportion that is available to induce a potentially toxic effect within the body. Bioavailability is determined by the characteristics of the element, soil, and biological organisms (Smith and Huyck 1999;Kim et al. 2015a, b;Bidar et al. 2020). Important soil properties include pH, organic matter, texture, cation exchange capacity, redox potential, oxide/hydroxide presence, and time (Petruzzelli et al. 2020). In performing risk assessments of element toxicity, bioavailability and bioaccessibility must be determined to avoid over-estimating the safety by using psuedo-total element concentrations alone (Izquierdo et al. 2015).
Psuedo-total concentrations are however useful for initial risk assessments that compare concentrations with national guidelines-although these guidelines vary greatly worldwide (Additional file 1: Table S2). If a garden soil, and/or the produce grown in it, contain elevated levels of a particular element, remediation actions may reduce concentrations to safe levels, for example by replacing the soil or constructing raised beds with safe media (Clark et al. 2008;Rai et al. 2019;Paltseva et al. 2020). Further, risk can be reduced through protection measures that reduce the longevity/intensity of exposure, limits extensive disruption of the soil, and/or restrict access to those at high risk (Briffa et al. 2020;Lupolt et al. 2021). Once concern advances to the stage where remediation is considered, for example if the concentration exceeds an intervention threshold, bioavailability and bioaccessibility studies would be highly appropriate to assess risk in greater detail. However, a viable risk may be present to some individuals in cases of repeated exposure (for example consumption of particular crops) for concentrations lower than established thresholds.
The indicated interactions in the simplified model shown in Fig. 1 affect the final outcome of the complex system, which is toxicity or illness associated with exposure to particular PTEs in gardens. Again, potential illness largely stems from ingestion/inhalation/contact with high levels of various elements in the soil and the consumption of plants with high levels of PTEs, but in consideration with their bioavailability/bioaccessibility (Rai et al. 2019;Briffa et al. 2020). One relevant feedback in the system is associated with the state of health, in the sense that a variety of types of element toxicity reduces one's health, yet various aspects of an individual's health may provide some degree of resilience (Monachese et al. 2012;Mitra et al. 2022). Nevertheless, most avenues for reducing the negative health effect are behavior risk reduction, site remediation, and appropriate garden management and mitigation actions .
Fundamentally, understanding the underlying soil and elemental properties is an important process in governing many of the processes of the human exposure chain depicted in Fig. 1 (Selim and Amacher 1996;Violante et al. 2010;Li et al. 2022). Thus, the baseline information required in investigating this issue is knowledge of the PTE concentrations in garden soils in context with the behavior of the at-risk group-this type of assessment aligns with the scope of our work and preliminary investigation (Goh 2018). Even in areas where health risks have not been reported, new and site-specific knowledge is needed to make informed decisions regarding potential risks to exposure (Chaney et al. 1984;Kim et al. 2014;Cooper et al. 2020). Once obtained, this information helps guide additional follow-on work to fully assess Page 10 of 29 Goh et al. CABI Agriculture and Bioscience (2022) 3:60 health risks to PTE exposure in contaminated urban garden soils (cf. EPA 2007).

Objectives
In this study we were concerned with assessing the presence of PTEs in the soils of community gardens in highly developed Singapore where urban gardening in housing estates is being promoted to encourage a growing, aging community to be active and selfsufficient (Tan and Neo 2009). Community gardens are also anticipated to supplement and diversify local food sources in the city (Chandran 2019). To that end, the nation aims to double the number of community gardens and allotments from 1500 to 3000 plots by 2030 (Begum 2020). The growing prevalence of gardens within public housing estates in Singapore that are visited frequently by onsite residents potentially represent an uncertain health risk because the gardens are inherently located in a dense matrix of industrial and residential land-uses with high road densities (Ng et al. 2006;Joshi and Balasubramanian 2010;Yuen et al. 2012). Further, other green landscapes such as public parks and accessible woodlands are also frequently visited locations of potential exposure to PTEs (Han et al. 2017). We focus on the presence of four PTEs (Cd, Cu, Pb, Zn) that are often associated with contamination in urban garden soils worldwide (Xia et al. 2011;Kaiser et al. 2015; Table 1). These elements are also among the most studied in past assessments of urban contamination set in Singapore (discussed below). Lead has a legacy presence in Singapore soils despite most leaded petrol being phased out in Southeast Asia in the 1990s (Lee et al. 2014), as well as the regulation of Pb in other important sources such as paint ). Our focus is on gardens constructed on soils; soil-less, raised bed, rooftop, and hydroponic systems are not considered. We also examine key soil properties (pH, Total Organic Carbon (TOC), texture) that may influence the mobility of elements in the soil (Zwolak et al. 2019). Absent in our assessment with respect to the exposure model (Fig. 1, "Human urban gardening exposure chain" section) are the determination of bioavailability, bioaccessibility, and PTE presence in vegetables. In foreshadowing our results, we did not perform these advanced experiments because pseudo-total element concentrations in soil were not alarmingly high. Finally, we use Pb to track the potential anthropogenic pollution sources of element contaminants in the Singapore urban soils (Lee et al. 2021).

Singapore geography
Singapore has a tropical equatorial climate (mean annual temperature = 27 °C; total annual rainfall = 2300 mm) with two monsoon seasons and short dry periods. The geology of Singapore is somewhat diverse, with major formations including igneous rocks such as granite and norite, various sedimentary rocks (sandstone, mudstone), isolated metamorphic rocks, Quaternary sediments, as well as marine clay deposits (Sharma et al. 1999). Soils are generally acidic (pH 3.7-6.2), have low cation exchange capacities (maximum mean value: 100 mmolc/kg), and have low P concentrations (mean values ≤ 0.28 g/ kg) irrespective of geology (Leitgeb et al. 2019). Soils developed on sedimentary rocks and sediments tend to have higher concentrations of soil organic carbon, total P, base forming elements, exchangeable base cations and moderate levels of some PTEs, for example Cr, Cu, Ni, Pb, and Zn (Leitgeb et al. 2019). Nguyen et al. (2019Nguyen et al. ( , 2020 determined that relatively high concentrations of PTEs may be associated with particular types of geological material (e.g., norite) and native soils.
Singapore's land-use is characterized by a dense matrix of urban areas with some lands dedicated to green spaces (Fig. 2). These areas differ somewhat in the degree of human disturbance. Woodland areas are made of (replanted) secondary forests located in either nature reserves or urban parks (Gaw and Richards 2021;Ramchunder and Ziegler 2021). Some have been disturbed by old village settlements, military activity, road building, and other activities requiring forest clearance followed by subsequent reforestation. However, these woodlands are not currently being disturbed due to strict land-use regulation. Open parks, in contrast, are grassland fields with few trees, but may be surrounded by substantial forested areas in some locations. Although they are considered open green spaces, soils tend to be compacted, both from construction and subsequent recreational activities where trampling occurs (Henderson 2013).
While some of the woodland soils may be in situ, the surfaces on all parks and urban gardens contain substantial fill material that was likely either imported from several quarries that were once operational in Singapore, or from other reclamation and sand mining projects (cf. Kog 2015). Thus, soils in gardens likely do not correspond with the underlying geology. Gardens are often built on marginal pockets of land sandwiched between housing blocks. Most of the construction of modern Singapore, including the spaces we investigated, took place in the last 50 years (Guo 2016). One garden we studied was built by top-filling over an old asphalt parking lot; another was built over a concrete basketball court. These fill soils tend to be deficient in important plant nutrients, which necessitates the addition of fertilizers or construction of raised beds using garden soils purchased from local manufacturers. Store-bought soils and inorganic fertilizers may originate from Singapore or be imported from overseas. Gardeners also apply mulches and other organic materials (eggshells, coffee grounds, compost) to enrich soil nutrients. These additions may modify the amount of PTEs present in the garden. Garden soil disturbances are mainly related to reworking of the soil with hand tools.

Singapore urban pollution
Several studies have examined selected PTEs in a variety of Singaporean soils/sediments, including forests (Nguyen et al. 2018(Nguyen et al. , 2020Leitgeb et al. 2019;Wang et al. 2022), landfills (Patra et al. 2017), mangroves (Cuong et al. 2005), water bodies (Sin et al. 1991;Chen et al. 1996Chen et al. , 2016Nguyen et al. 2018), marine sediments (Rahman et al. 1979;Goh and Chou 1997;Wood et al. 1997;Cuong et al. 2008), military lands (Nguyen et al. 2019), road dust (Joshi and Balasubramanian 2010;Joshi et al. 2009;Yuen et al. 2012); fly ash (Tan et al. 1997;Wu and Ting 2006); and various types of urban lands (Zhou et al. 1997;Ng et al. 2006). The studies of Chen (1999), Zhou et al. (1997), and Wang et al. (2022) Fig . 2 The map shows major land-cover/land-use types in Singapore, with sample locations shown as diamonds (community gardens), circles (public parks) and triangles (woodlands). (Adapted from a base map by Gaw et al. 2019). Below the maps are representative images for gardens (left), parks (center), and wooded areas (right)  Most values are ranges, means ± standard deviation, or medians ± MAD; a indicates a range of means/medians; ~ signifies the values were estimated from a graph. ND is not determined; bdl is below detection limit References: (a) This study; see "Methods and materials" section for details on the analysis methods  are particularly informative to our study because they show loading of Cu, Pb, and Zn in soils are higher in the proximity to industrial areas compared with natural or residential areas (Table 2). Chen (1999) found very high maximum values of Cu (485 mg/kg) and Zn (3594 mg/ kg) in industrial soils. In the most recent study, Wang et al. (2022) found maximum concentrations of Cu and Zn were an order of magnitude higher in soils adjacent to industrial activities, compared with parks. The main sources of high PTE concentrations in Singapore soils are likely vehicular phenomena and industrial activities (Zhou et al. 1997;Wood et al. 1997). For example, Yuen et al. (2012) reported roads in or near industrial areas had higher levels of PTEs in surficial dust compared with residential areas (Table 2). Joshi et al. (2009) reported very high concentrations of Cu (9069 ± 3742 mg/kg) and Zn (1696 ± 446 mg/kg) for road sediment in industrial areas. Wood et al. (1997) reported that elevated Pb and Zn levels in sediments of the Johor Straits along the Malaysia-Singapore causeway likely resulted from gasoline and tire wear, respectively. Fly ash from incinerators in Singapore have very high concentrations of PTSs (Table 2): Cu (570-087 mg/kg), Pb (2000-2150 mg/kg), and Zn (6288-8072 mg/kg). Potential contamination sources in Singapore are indicated in the map in Additional file 3: Fig. S2.
A prior study of airborne element-bearing aerosols in Singapore showed an enrichment of a variety of PTEs (Pb, Zn, Cd, V, Ni, Cr, and Cu) that originate from several sources including oil-fired power plants, element processing industry, land reclamation, construction activities, municipal solid waste incinerators, and traffic emissions (Balasubramanian and Qian 2004). Ng et al. (2006) determined that element accumulations in lichens across Singapore were the result of atmospheric deposition, with peak concentrations of Cu (45 mg/kg), Pb (17 mg/kg), and Zn (84 mg/kg) found in locations associated with heavy petroleum refinement, shipping industries, and road traffic, respectively. Rivellini et al. (2020) found evidence that trace elements in aerosols may originate from shipping emissions and the oil refinery industry. Once airborne, these contaminants can be deposited in the environment during rain events and via dry dust deposition (López et al. 2019). Hu and Balasubramanian (2003) determined that the removal of the trace elements from the atmosphere by precipitation was influenced by the rainfall volume as well as pH. Further, the magnitude of mean annual wet deposition fluxes of combustion-generated elements (V, Ni, and Cu) were high relative to other locations because of abundant rainfall. The southwest monsoonal winds may cause slight increases in PTE concentrations of airborne particles due to the presence of industrial estates and petrochemical plants in the southwestern region of Singapore Rivellini et al. 2020). Another potential source for the occasionally elevated concentrations of Cu and Zn found in residential areas of Singapore is the degradation of buildings materials.

Sampling and analyses
We adopt an exploratory survey of PTEs in soils across the breadth of the island of Singapore to eliminate spatial bias that may be related to proximity to probable sources of element pollution. In May-June 2018 we collected 300-g samples from 10 sites in the following three land uses (Fig. 2): woodland; parks; and community gardens (raw data reported in Additional file 1: Table S2). Soil pH was measured in a soil suspension with 1:5 soil:water ratio using a calibrated 9157BNUMD RossUltra pH/ATC Triode with AgCl-Hg electrolytes (Thermo Fisher Scientific, Waltham, USA). Soil particle size was determined using the Mastersizer 2000 (Malvern Instruments, Malvern, UK) via the laser obstruction method. Total organic carbon was determined by combusting 0.03 g of ground soil that was first pre-wetted, fumigated for 6 h to remove inorganic carbonate, and then oven-dried. The samples were combusted at 950 °C in the Vario TOC Cube (Elementar Analysensysteme GmbH, Langenselbold, Germany).
Soil element psuedo-total concentrations were determine using the PerkinElmer Optima 8300 Inductively-Coupled Plasma Optical Emission Spectrometer (ICP-OES). Prior to analysis we used a microwave-assisted acid approach for digestion (US EPA Method 3051A17). First, 0.5 g of the finely ground material was placed in digestion vessels before being reacted with 9 ml of Nitric acid (65%) and 3 ml of Hydrochloric acid (37%) in a fume hood for 15 min. The samples were then heated in a microwave digester, ramping up to 175 °C within 10 min, then maintaining that temperature for another 15 min. Multi-elemental standard solutions of 0.5 mg/kg, 1 mg/kg, 5 mg/kg and 10 mg/kg in 2% HNO 3 were used to calibrate the ICP-OES instrument.
We used the USGS 2711a Certified Reference Material to perform quality control and quality assurance of both the accuracy and repeatability of PTE measurements (Table 3). The percentage recovery of PTEs in this research (% Lab recovery) is generally higher than the recovery rates from the USGS Contract Laboratory Study Programme (CLSP). We found higher recoveries compared to those obtained in the USGS CLSP, which mainly used the weaker EPA Method 200.7 (US EPA 2001) for soil digestion via hotplate reflux (Mackey et al. 2010). Microwave-assisted digestion results in higher recoveries because the high temperatures and pressures achieved in microwave vessels allows for a more complete digestion of the soil's PTEs (Chen andMa 1998, 2001). Closed microwave vessels also reduce the potential loss of trace elements through volatilization, which is a feature in open hotplate digestion (Link et al. 1998).
The accuracy and bias of Cd, Pb, and Zn are acceptable; however, the over estimation of Cu (34%) is noted and considered in our interpretation of the results. We deemed that ICP-OES and di-acid digestion approach was acceptable for the expected range of total concentrations that would indicate potential risk (Olesik 2020). Over the years, a variety of methods have been used to determine pseudo concentrations of PTEs in garden soils (summarized in Additional file 1: Table S1).
Finally, we determined lead isotope ratios in five samples to investigate the origin of Pb at selected sites (Komárek et al. 2008). In doing so we assume that Pb is a representative indicator of general pollutant dispersal pathways (Wiederhold 2015;Harvey et al. 2017;Carrasco et al. 2018;Crispo et al. 2021). We analyzed the soil from five sites, one woodland (#4) and two each for parks (#1 and 7) and gardens (#2 and 10), to track likely sources of Pb in the three community spaces (see locations on Fig. 2). The Pb isotope values were determined by the Environmental Chemistry and Microbiology Unit at the Charles Darwin University in Darwin, Australia. The samples were selected because they were of relatively high concentrations. The finely-ground samples underwent acid digestion with HNO 3 and HClO 4 . Isotope values were determined using an Inductively-Coupled Plasma Mass Spectrometer.

Calculations and statistical handling
As the number of soil samples collected at each site were few (n = 10), we used nonparametric statistics for data inference, including the Spearman rank correlation test, the Kruskill Wallis (KW) for multiple treatment comparisons, and the Mann-Whitney U (MW-U) as a post-hoc test. Differences were considered significant at α = 0.05 (P-values < 0.05). We used an Enrichment Table 3 Quality assurance/quality control information for Al, Cd, Cu, Pb, and Zn determinations for one standard reference material (USGS 2711a) Certified value is the concentration reported by the SRM manufacturer; analyzed is the mean ± Stdev determined by ICP-OES (this study); recovery = analyzed value/certified value x 100%; bias = ((analyzed value − certified value)/certified value) × 100; precision = Stdev/mean (ie., coefficient of variation of analyzed values)

Cd Cu Pb Zn
Certified value (mg/kg) 54.1 ± 0.5 140 ± 2 1400 ± 10 414 ± 11 Analyzed value (mg/kg) 58.6 ± 0.8 190. Factor (EF) to assess the anthropogenic enrichment of PTEs in the soils of different land covers (Chester and Stoner 1973). Although the validity of EF values is debated (Reimann and Garrett 2005), we find them useful for indicating sites of potentially high contamination (Szolnoki et al. 2013

Soil element concentration risk guidelines
Finally, as Singapore does not have published guidelines pertaining to health risks associated with concentrations of PTEs in soils, we determined three thresholds based on a review of guidelines for nearly 50 countries/ regions worldwide (Table 4) (Table 4) are determined as the 20th, 50th (median), and 80th percentile of the values summarized globally (Additional file 1: Table S3). The ERC values are of approximately the same magnitude as those employed by Wang et al. (2022) as "intervention" values. We caution that while these guidelines are useful for examining potential levels of risk in our assessment, they may not be appropriate for policy/ management applications without a more critical review.
In Tables 1 and 2 we list studies based on a firstorder ranking of level of contamination. For each study we calculate a composite score that is determined by awarding 1 point if the maximum value of an element exceeds the TRC; 2 points if the PRC is exceeded; and 3 points if the ERC is exceeded. Four points are assigned if the maximum is tenfold the ERC (twofold for Cd). For studies, not reporting maximum values we use the mean (median) + stdev (MAD) as the upper value. If an element was not determined (ND) we adjust the ranking to align with level of contamination relative to other studies. These rankings are made to facilitate discussion in this paper, not to compare cities/locations directly (see below).

Pseudo-total concentrations
The main difference between PTE presence in soils in gardens, parks, and woodlands was in concentrations of Table 4 Number of samples exceeding threshold risk concentrations (TRC), potential risk concentration (PRC), and elevated risk concentrations (ERC) for four metals and three treatments a The pseudo-total concentration indicators of level of risk are determined as the 20th (TRC), 50th (PRC), and 80th (ERC) percentile of the ranked concentrations from values reported for numerous countries worldwide (summarized in Additional file 1: Table S2). In the bottom panel, n = 10 for each land cover  Cu and Zn ( Unlike Cu and Zn, Cd was rarely detectable in any sample (detection limit = 1.6 mg/kg Cd): three garden soils had Cd concentrations ranging from 1.7 to 3.2 mg/ kg Cd; one woodland soil sample had a relatively high value of 6.4 mg/kg Cd (Table 5). No detectable Cd was found in park soils. Similarly, Pb was also rarely present in concentrations above our (high) detection limit of 12.6 mg/kg (n = 5, 5, and 4 in gardens, parks, and woodlands, respectively). Detectable Pb concentrations ranged from 13 to 61 mg/kg Pb across the three treatments (Table 5). The maximum concentrations found in gardens, parks, and woodlands were about 36, 61, and 33 mg/kg Pb, respectively. The generally low PTE concentrations we found are generally comparable to those reported by Chen (1999), Zhou et al. (1997), Ng et al. (2006, Nguyen et al. (2018Nguyen et al. ( , 2020, and Wang et al. (2022) for a variety of urban soils in Singapore.

PTE enrichment
In comparison with elemental background concentrations from soils in the forested Nee Soon Catchment, the Cu and Zn values found in soils of the three land uses investigated are considered "enriched" to some degree. For example, the median enrichment factor (EF M ; Eq. 1) for Cu in gardens was 7; median EF M values were 3 and 4 in parks and woodlands, respectively. Simply, these values indicate the medians of the samples measured in gardens, parks, and woodlands are all threeto sevenfold higher than the background concentration of 5 mg/kg Cu. Enrichment factors for Zn ranged from 5 in garden soils to 2 for both park and woodland soils (background = 22 mg/kg Zn). Here, enriched does not necessarily indicate high risk. It may simply indicate an anthropogenic contribution; alternatively, it may simply the particular mineralogical composition of the soil-which again may be imported or natural for that particular location.

Correlations
The garden soils were slightly sandier (median of 36%) than those in parks or woodlands (24 to 30%; Additional file 3: Table S4). Copper and Zn concentrations were negatively associated with clay content in garden soils (ρ = − 0.663 and − 0.725, respectively; α < 0.01; Additional file 3: Table S4). In the woodlands, Cu concentrations were positively correlated with silt percentage (0.674), and negatively correlated with the sand fraction (− 0.609; α < 0.01). This same textural association was present to a lesser degree for Zn in woodland soils: silt (ρ = 0.371; α < 0.05); and sand (− 0.255; α < 0.05). The pH of garden soils ranged from 5.8 to 7.8, with a neutral median value (6.9 ± 0.6; Additional file 3: Table S4). In comparison, the soils of the other two land uses were more acidic (median pH was 5 in the woodland soils; 6 in park soils) and displayed larger ranges: 4.8-8.0 (parks); 4.0-8.6 Table 5 Summary statistics for PTE pseudo-total concentrations and soil properties for the three land covers Values reported are medians ± standard absolute deviations (n = 10 samples). Ranges are reported in brackets. bdl refers to below detection limit, which are 1.62 (Cd), 5.82 (Cu), 12.6 (Pb), and 3.54 (Zn). Values in rows with different letters in "()" are significantly different, based on post-hoc testing (KW, followed by MW-U; α = 0.05) *The differences among landcovers is significant at α = 0.5 (KW Test), with garden values being significantly higher than the other treatments (KW; α = 0.05). Values in rows with different letters in "( )" are significantly different, based on post-hoc testing (MW-U; α = 0.05)

Element
Units Gardens Parks Woodlands (woodlands). The only significant correlation involving pH was a weak correlation with Cu concentration in garden soils (0.387; α < 0.01; Additional file 3: Table S4). None of the soils had particularly high TOC, with median percentages ranging from 3.4 to 4.5 (Table 5). The highest value was found at a woodland site (9%); low values < 2% were found among all three treatments. TOC concentrations were positively correlated with Cu and Zn in garden soils (ρ = 0.761 and 0.667; α < 0.01; respectively), as well as park soils (0.447 and 0.667, but with different levels of significance). Collectively, these significant correlations between element concentrations and soil properties provide some indication of elements being bound to the organic matter or particular texture fractions in soils, but without a consistent pattern across land uses. Unexpected is the negative association with clay content (Additional file 3: Table S4).
With regard to the PTEs, the correlation between Cu and Zn was significant (ρ = 0.78; α < 0.001; Fig. 3). Correlation analysis involving Cd and Pb was hindered by so many samples having concentrations below detection limits.

Fig. 3
Relationship between copper (Cu) and zinc (Zn) in garden, park, and woodland soils in Singapore. The Spearman rank correlation coefficient is ρ = 0.78 (p-value < 0.001) Fig. 4 Triple Pb isotope plot of Pb sources in Singapore, which includes roadside sediments in blue circles (Yuen et al. 2012), aerosols in orange crosses (Lee et al. 2014;Chen et al. 2016), incinerator ash in red diamonds (Chen et al. 2016), gas station soil residuals in brown asterisks, Singapore coastal water in purple triangles and urban runoff in green circles (Carrasco et al. 2018), MacRitchie reservoir sediments (Chen et al. 2016) and Nee Soon Catchment soil samples (unpublished data) Page 18 of 29 Goh et al. CABI Agriculture and Bioscience (2022) 3:60 Lead isotope ratios Chen et al. (2016) suggested that due to the close proximity of industrial estates with incinerator plants, Pb pollution in Singapore usually comes from multiple sources including industrial activities, incineration, the past use of leaded petrol, and the weathering of parent rock material. The Pb isotopes for garden site G2 were similar to those found in urban runoff from petrol stations and shipyards (Fig. 4). This result suggests that even though Pb was banned from petrol in 1998 (Chen et al. 2016) the legacy of Pb contamination has persisted, probably because of the low solubility (Sutherland and Tolosa 2001;Yuen et al. 2012). Abrasive blasting of element surfaces in shipyards during repair and construction works are another potential source of Pb aerosols. Considering the down-wind location of G2 from the Tuas shipyards, the aerial transport of Pb particles to the site is plausible. Wang et al. (2022) reported high concentrations of selected elements in the plant foliage in this areaalthough lead was not elevated. Again, potential sources of pollution are shown in the map in Additional file 3: Fig. S2. Sample W4 is similar to the roadside sediments collected from Tuas and Senoko Industrial Estates and the incinerator from Keppel Seghers (Fig. 3). Again, incinerator ash from Singapore contains high concentrations of Pb (≥ 2000 mg/kg; Table 2). The presence of Pb at W4 could therefore be due to the upwind (S and SW winds) location of the Keppel Seghers incinerator and the airborne transport of Pb in fly ash from industrial activity in Tuas and West Coast, followed by deposition (Hu and Balasubramanian 2003). Industrial emissions of Pb also include plastics production and element processing plants (Ragaini et al. 1977;Siddiqui et al. 2009), which are found at the Tuas (southwest) and Senoko industrial estate in the north.
The Pb isotope ratio at park site P1 is similar to that of coastal waters of the Central Business District, Kranji and Sembawang (Fig. 3). Rahman et al. (1979) state that element accumulation in Singapore waters and marine sediments must come from natural sources or be introduced from land-and marine-based activities, including leakage from landfills, oil/fuel spills, and/or runoff from shipyards, wood mills, motor works shops, and small industries. Lead concentrations in marine sediments at seven locations ranged from 20 to 28 mg/kg Pb (Rahman et al. 1979). The Pb ratios for P1 likely reflect the mixing of natural Pb and that associated with landbased runoff; however, this inference has uncertainty.
Two of the soil samples tested did not contain an anthropogenic Pb signature (Fig. 3). Isotopic ratios for sites P7 and G10 are close to those found at the base of the MacRitchie Reservoir sediment core, which is believed to reflect naturally occurring Pb in Singapore that was transported with eroding material from the forested headwaters. The G10 signature is most similar to those from soils in the relatively undisturbed Nee Soon Catchment on the island.

Spatial patterns
The observed median values of Cu (18-30 mg/kg) and Zn (50-52) for the three treatments are at the high end of the range of central tendency values (reported means, median, midpoints of ranges) associated with soils in forests and natural spaces determined by others for Singapore (Table 2; Fig. 5). The median concentrations for garden soils were slightly higher (43 and 125 mg/kg Cu and Zn, respectively), and therefore, resembled the lower values associated with disturbed soils in Singapore (Table 2). Comparison with reservoir sediment values supports the idea of some degree of anthropogenic enrichment above the natural geochemical signal on the island. With respect to geography, the highest Cu concentrations had no obvious pattern, while Zn concentrations tended to be higher in the western and south-eastern region of Singapore (Fig. 6). Higher Pb concentrations were found in the northern part of Singapore, although a more obvious relation was the occurrence of detectable concentrations in park soils (n = 3) versus garden (1) and woodland soils (1). Again, most Cd concentrations were low and had no strong geographic pattern. Elevated concentrations of elements in samples from the western and northern area of Singapore were somewhat expected because of the close proximity to industrial element processing, wood combustion, and waste incineration facilities, which are shown in gray shading in Fig. 6b (Nriagu 1979;Wang et al. 2022). Prior analysis of emissions from two incineration plants in the southwest of Singapore showed that high levels of Cu and Zn attached onto particulate matter (Tan et al. 1997), which may lead to elevated concentrations via deposition. We also expected some enrichment in this area because of the role of the southwest monsoonal winds blowing over industrial estates and petrochemical plants . Of note, five of the western-most sites had concentrations of Cu and Zn that were higher than their respective TRCs ( Fig. 6d; Table 3).
The great variability in Pb isotopic ratios (Fig. 4) indicates different potential sources, both natural and anthropogenic (Fig. 4). The absence of common sources  Table 4). The background in each map varies to allow comparison of the study sites with the road network (A), major parks, industrial areas, water bodies and shipyards (B) and large geological formations (C). The base maps in A-C do not necessarily reflect an impact on the element that is mapped Page 20 of 29 Goh et al. CABI Agriculture and Bioscience (2022) 3:60 for lead aligns with the general lack of spatial pattering in the concentrations (Fig. 6). Relevant is that many of the samples enriched in Cu and Zn found in prior studies were often associated with industrial activities in similar locations (Zhou et al. 1997;Chen 1999;Wang et al. 2022). Both our results and those of Wang et al. (2022) suggest that the footprint of many potential pollution sources is relatively small-and therefore contribute to the "spotty" nature of enrichment across the island.

Singapore urban garden soil contamination context
The apparent combination of several activities/processes that contribute to the variable enrichment of particular PTEs makes Singapore different from several studies conducted in very contaminated urban environments (Table 1). Briefly, the sites in Table 1 (Table 1). Again, one should be careful in making literal distinctions among sites/studies owing to differences in methods, timing, and rationales for each study. For example, many early studies were conducted in areas where contamination from known activities was. Generally, soil PTE concentrations in gardens worldwide tend to be affected by legacy contamination from the following (Table 1): leaded paint and various construction materials used in old neighborhoods (e.g., Baltimore, Christchurch, London, New York, New Orleans, San Francisco, Washington DC), smelting/ mining (Baia Mare, Dalnegorsk, Gijon, Palmerton, Paris area) and unsafe waste disposal/usage in gardens as irrigation or fertilizer (Kampala, Paris, Marseille, New York, New Castle). The proximity to heavily used roads may also be an important source of Cd, Cu, Pb, and Zn in gardens, for example, as is the case in Bristol, Cardiff, Hong Kong and New Orleans (Table 1). Older roads may have high levels of Pb stemming from past leaded petrol use.
Noticeable in Singapore is the generally low PTE concentrations in garden, park, and woodlands soils, compared with road-associated sediments and fly-ash ( Table 2). As concentrations of Cu and Zn decrease dramatically with increasing distance from roads (Chen et al. 2010;Kim et al. 2017), the moderate PTE concentrations we found are plausible because many of these spaces not often located directly adjacent to roads. Gardens are often ensconced within housing complexes and occasionally sheltered by buildings and trees.
Enriched values in garden soils may also be result from degradation of building materials associated with the residential units (Chaney et al. 1984;Lim et al. 2021).
Many authors in the reviewed studies refer to general urban atmospheric deposition of elements in dense cities with much industrial activity and traffic congestion. This type of "new urban deposition" phenomena is likely an apropos description explaining the variations we found in Singapore. In support, Balasubramanian and Qian (2004) reported that the enrichment of aerosols in Singapore was more than 1000 times for Cd, Pb, and Zn; 100 times for Cu. Deposition is almost certainly superimposed on the signatures of the garden soil and may have a legacy signature associated with the age of the building complex. Management activities might be partly responsible for our finding of enriched Cu and Zn concentrations despite the coarse nature of the garden soils, as reported worldwide (Table 1).
For example, the high TOC percentages we found in gardens, which are maintained by application of composted dead organic matter containing elementfixing humic acid, potentially contribute to the enrichment of Cu and Zn in surface soils (Veeken et al. 2000;Amir et al. 2006;Sánchez-Monedero et al. 2002). In addition, NPK fertilizers may facilitate the leaching of Cu and Zn from the topsoil (Young 2013;Bolan and Duraisamy 2003;Loneragen and Webb 1993;Carbonell et al. 2011;Rajneesh et al. 2017;Warman and Cooper 2000). Further, Cu-and Zn-based fertilizers often contain insoluble compounds (e.g., CuSO 4 , ZnO, and ZnSO 4 ) that potentially elevate Cu and Zn concentrations (Tsang et al. 2014;Jones and Belling 1967). Garden patrons often remedy acidity by applying Ca-rich fertilizers including eggshells to raise the soil pH (Ok et al. 2011), which may in turn facilitate the precipitation of Cu and Zn compounds in the topsoil (Degryse et al. 2006;Kabata-Pendias 2011;McBride 1994). Fertilization may also elevate concentrations in park soils. Focused research would be needed to understand these processes and relationships more thoroughly.

Environmental risk of gardening
In general, comparison of the pseudo-total concentrations of PTEs with risk thresholds do not indicate substantial risk. In the Singaporean context where gardening is an activity that is promoted as beneficial for the growing, aging community to be active and self-sufficient, our analysis indicates limited risk of exposure to contaminated soil. Most Singapore gardeners only spend a few hours per day in their gardens (Goh 2018); and therefore, inhalation and exposure rates are very low (Additional file 3: Table S5). For completeness, however, those concerned about health risks would Page 21 of 29 Goh et al. CABI Agriculture and Bioscience (2022) 3:60 also need to consider all exposure routes, including the consumption of produce grown in gardens (visualized in Fig. 1). Others report that consuming the products of community gardens comprise about 60% of overall PTE exposure because some garden produce can bioaccumulate elements via root intake or foliar uptake (Spliethoff et al. 2016;Finster et al. 2004). Plants can potentially accumulate up to 10 times the concentration of Cd and a third of the concentrations of Cu, Zn and Pb in the soil (Kabata-Pendias 2011; Byrne et al. 1976). Dalenberg and van Driel (1990) state that atmospheric deposition is a significant source of Pb for field crops, accounting for 73-95% of total Pb found in plants.
Gardener activities such as applying fertilizers, pesticides, mulching, or soil amendment may influence the presence of elements in the soil (Mortvedt 1996;Atafar et al. 2010), but our analyses did not reveal a consistent association between element concentrations and soil variables in Singapore urban gardens. Nor did we identify sources of substantial contamination. Even for such cases of low levels of contamination, the presence of pollution sources (e.g., industry, road networks) necessitates that urban garden activities should be monitored from time to time to safeguard and inform the public of the level of risk, even if none (cf. Chaney et al. 1984;Kim et al. 2014). Patrons of parks and gardens should be aware that the presence of elements in the soil are a function of background residual signatures and potential enrichment via atmospheric deposition, surface runoff deposition, and application of soil amendments.
Finally, with respect to adverse health risk from consuming elements, bioavailability/inaccessibility should be considered to assess the potential for a toxicity effect if concentrations in soils and in farmed produce are high (Intawongse and Dean 2008). Again, we did not perform this assessment because of low pseudototal PTE concentrations. We assume that the Singapore context is analogous to that of reported by Izquierdo et al. (2015) in Madrid with slightly higher levels of PTE contamination. They concluded that only in the worstcase scenario of children playing in urban gardens and eating produce grown from them, would risk exceed limits of acceptability. This assessment is in agreement with those of Sialelli et al. (2011) for urban soils in Torino (Italy), and Sipter et al. (2008) studying contaminated vegetables in Gyongyosoroszi (Hungary), but not of Bielicka-Giełdoń et al. (2013) in Koszalin (Poland) for conservative guidelines.

Limitations
We recognize that our elemental recovery methods were not sensitive enough to determine many low-level concentrations associated with Cd and Pb samples. While this issue did not affect our ability to assess environmental risk to these elements, it did hinder our attempt to identify associations between all elements and soil physical factors. Further the recovery of Cu from the certified reference material was high, indicating our Cu determinations are potentially elevated. However, we did not determine the observed pseudo-total concentrations of Cu to be unrealistically high. With regard to the methods employed, others working in Singapore have had better luck deriving total concentrations using the techniques employing strong or mixed acids and analysis with an ICP-MS, particularly for low concentration elements such as cadmium (Yuen et al. 2012;Nguyen et al. 2018).
We are also somewhat uncertain of how the observed associations in soil pH, TOC, and texture play a role in element retention in garden soils of Singapore. For example, the negative correlation between elements (Cu and Zn) and clay content is inconsistent with most reports of higher element concentrations in clayey soils due to the greater presence of negatively-charged sites for element cation adsorption (Scokart et al. 1983;Quenea et al. 2009;Huang et al. 2014;Thomas and Lavkulich 2015). Nevertheless, some studies have reported greater PTE adsorption on coarse sediments in roadside dust and soils (Varrica et al. 2003;Wang et al. 2006;Acosta et al. 2011).
Finally, we were also limited by the inability to source the location of the soil material used in the construction of the gardens and parks. Doing so would have allowed us to fully understand element enrichment over the background concentrations that have taken place since construction of the gardens. Related, we also do not know the age of the gardens, with respect to element accumulation time.

Conclusions
While provenance tracking using lead isotopes indicated a range of potential sources for PTEs entering urban gardens, measured pseudo-total concentrations of Cd, Cu, Pb, or Zn did not indicate substantial health risk to those practicing urban gardening in Singapore. Further, estimated intake rates determined in a preliminary study (Goh 2018; Additional file 2: Table S5) were much below guideline specifications of maximum allowable doses, although the method to estimate these rates was coarse. Given the low concentrations of PTEs in the garden soils, we did not expand the analysis to determine element concentrations in produce grown in the gardens or the bioavailability/bioaccessability of the elements if ingested or inhaled. However, when combined with findings from similar studies performing these assessments worldwide, the observed low pseudo-total concentrations observed leads us to conclude that urban gardening activities represent a viable and safe form of leisure, with the potential to produce food for consumption in Singapore. Risk would increase if vegetables in very contaminated soils were consumed, particularly by small children. These positive results add to the discussions on the viability and sustainability of urban agriculture in most urban settings-in agreement with Brown et al. (2016) and Lupolt et al. (2021) that the benefits often outweigh the risks. With respect to urban contamination in general, the growing work in Singapore, summarized in Table 5, supports the notion of a highly heterogenic signal of PTEs, whereby high concentrations associated with particular point sources, especially roads and industrial activities, are mapped onto a variable natural signal that has probably been also been enriched diffusely by the deposition of PTEs over time.