Salary Satisfaction of Employees at Workplace on a Large Area of Planted Land

: Salary satisfaction (SS) perception by employees can be affected by psychological impacts from the workplace setting. Landscape attributes of green and blue spaces (GBS) may account for this effect, but relevant evidence is rarely veriﬁed. In this study, a total of 56 Chinese industrial parks were chosen as study sites, where employee satisfaction was assessed by rating facial expression scores (happy, sad, and neutral emotions) in photos obtained from social networks (Sina Weibo and Douyin). The structures of the GBSs were characterized remotely by largeness of size, height, and visible ratio of green view (GVI) in a 2 km radius buffer area around the workplace. Street view images from Baidu map were selected for estimating GVI using a pre-trained deep learning model and botanical experts evaluating woody plants’ diversity. The results indicated that SS can be estimated with the maximum likelihood analysis model against the happy score, which ranged within 8.37–18.38 (average: 13.30 ± 2.32) thousand RMB. A regression model indicated SS was lowered by a larger green space area in agreement with a reduced happy score. Further, sad scores in highland areas with tall plants and a strong depression on the happy score was associated with a greater plant diversity. Interesting from this study, the designed apparent size of green space should be considered in green space construction near a workplace to prevent perceptional decline towards SS, while blue space is irrelevant in this relationship. Similarly, the diversity of woody plants should be planned to control its negative impact on the perception of positive emotions, with plant diversity beyond a comfortable level perhaps further decreasing SS.


Introduction
Industrial parks account for an inevitable proportion of built-up lands in municipal regions [1,2].An industrial park is not only a collection of corporations that shoulder regional household incomes, but also an activated source of regional heat and pollutant emissions [2].Urban landscapes with industrial parks are usually surrounded with a planned level of surrounding green-blue spaces (GBS) to mitigate environmental contamination, which can vary in density [3,4].Employees are major users of GBSs at workplaces [5], and nearly ~40% of them may spend time outdoors with GBSs [6].As a part of urban nature, GBSs at workplace can be experienced as an infrastructure that provides ecosystem services and natural elements to reduce the mental stress of workers [7][8][9].For example, office employees reported that they can perceive stress reduction at workplaces when accessing a GBS in an open outdoor space [7,8,10].The reduced noisiness of GBSs near a workplace may influence the overall perception of employee well-being [7,11].Thus, it is worthy to determine the mechanisms for promoting mental health and well-being through the landscape planning of GBS near workplaces.People perceive the natural and planned world in a variety of ways.An overall configuration of GBSs is more perceptible as a mental health promotion than any single landscape metric [12,13].For community-level landscapes, the aggregation of neighborhood GBSs has a greater impact on social cohesion and perceived happiness than any specific metrics such as building size, complexity, and fragmentation [14].Agglomerated public-space landscapes and associated GBSs (e.g., ≤9 ha) that overlap a dense road network can be perceived as benefits by people for cooling comfort and promoted mental well-being [12,15,16].Urban park visitors perceive mental well-being due to the awareness of experiences with landscape metrics of largeness [17,18], elevation [19,20], and location [21].The humanbiodiversity interaction may account for these benefits through delivering an experience with various species in nature that provide gains in perceived well-being [22].Plant biodiversity in green spaces likewise affect subjective well-being [23,24].Species diversity in blue spaces was also reported by perceived stress reduction [25,26].The current methodology for assessing plant diversity in urban GBSs includes in situ investigation [27], evaluation from a tree inventory [28], and a literature review [24].The geographical range of diversity in GBSs has been investigated from a local stand to national [24] and global scales [28].However, the range of human population for subjects' sampling has not matched these investigative scales.This is mainly because of limits in conventional methodology for data collection to assess mental perception.
The green view index (GVI) was created to evaluate the visible proportion to account for the dose of green nature in a street view [29,30].It is a flexible parameter that can be estimated using images such as street view images [31,32].Deep learning as a method has been well established for rating GVI from street view photos taken at locations around urban ecological infrastructures [30,33].These, in conjunction with place-based coordinates, can be used to locate points of interest (POI) in urban surface features.The ability to integrate programed data collection (e.g., crawling) is helpful for efficiently pooling a set of digital photos, which can be further combined with plant diversity data.The urban nature that residents can sense (e.g., through touch, sound, smell, sight, and maybe taste) is a way of perceiving GBS.These perceptions can vary by location and level of urban greening.A street view photo that can be used for rating GVI can also be a source of plant diversity data that accounts for the green color.To our knowledge, these novel techniques have rarely been used for detecting the possible driving forces of GBSs that affect the mental well-being of employees in industrial parks.
The Stress Reduction Theory (SRT) indicates that an experience with nature can reduce mental stress and evoke a restoration [34].This suggests that highly stressed people, such as employees, are a typical population who can gain perceived restoration from stress reduction.A continuous expenditure of direct attention through work affairs results in mental fatigue and a perceivable decline in well-being [35,36].The Attention Recovery Theory (ART) asserts that exposure to natural environments encourages more effortless brain function, thereby allowing it to recover and replenish its directed attention capacity [37,38].According to these two theories, people who are frequently exposed to GBSs can obtain promoted health and well-being [39,40].Exposure to GBSs means a closeness with nature, which has several benefits for visitors, including but not limited to better social interaction [41], solid social cohesion [42], improvements in psychological state [18,21,41], and improved physical health [43].These benefits result from ecosystem services provided by GBSs through mitigating air pollution [44][45][46], cooling the surrounding temperature [47], and adjusting the regional climate [19].Studies using facial expression scores have fully identified these positive effects with more details about the involvement of driving forces, such as landscape metrics [17,18,20,21], plant diversity [24], regional climate [19,48,49], and air quality [44][45][46].These parameters are also recommended for testing the SRT and ART effects on the emotional perceptions of employees exposed to GBSs at their workplace.
Salary satisfaction results from perceptions against depression or disparity about salary, which is a critical issue that attracts concerns in the evaluation of employees' work quality [50,51].Salary dissatisfaction is evoked from the accumulation of negative emotions in association with job insecurity [52,53].Dissatisfaction about salary may even lead to serious psychological disorders and mental diseases [51,54,55].Conventional evaluations of salary satisfaction depend mainly on data collected from questionnaire-based surveys that are highly labor-and time-dependent [56,57].This methodology has several limits that hinder the estimation accuracy from a desired level.Firstly, human bias cannot be avoided when a subject is reporting with a direct attention either by their inability to accurately rank and sometimes though paid respondents that just complete the work with little thought for answering [34].Occupation-related mental fatigue is a common response of employees which highly impacts their subjective well-being and satisfaction about income [35,36].Secondly, workers may inhibit their true feelings and concerns by subjective self-reserves according to the norms and values required in a group [34].In a cohort of workers, everyone can clearly recognize the individual role in a project and attach subjective well-being to the cohesion with leadership [36,58].In this social relationship frame, self-reported scores by employees are of high suspicion for their accuracy.Finally, many questionnaires have been designed and used without a formal process of validation [59].Because of these drawbacks, questionnaire-based surveys are a challenge for collecting useful for data in a large-scale investigation with a big cohort of subjects recruited.This is a reason why most of the current studies using questionnaires can draw conclusions from one or a small number of industrial parks [7][8][9].It is necessary to employ new methods and techniques for data collection, such as the evaluation of salary satisfaction.
Salary satisfaction results from emotions towards occupational experiences.Emotion is variable and fluctuates following subjective perceptions towards work affairs that may impact employees' perceptions and reduce their salary satisfaction [60].Efforts have been made to quantify emotional responses to salary satisfaction [60][61][62].Compared to self-reported emotional state, workers expose their emotions more directly through facial expressions [63,64].Facial photos or selfies are the major sources of facial expression data which can be used for assessing spontaneous emotions (facial expressions exposed without awareness of being photographed) [65][66][67] or posed emotions (facial expressions exposed with a conscious intention) [19,44,45].Self-reported perceptions can unlikely be spontaneous with regard to self-reported perceptions of salary, because the respondent is fully conscious about purpose and consequence when giving a reply [56,68,69].When emotions are exposed without any clear guidance, they can be characterized to be spontaneous.These together suggest that salary satisfaction can be estimated by modeling against facial expression scores.Spontaneous expressions are recommended firstly as independent variables for estimation, but posed emotions can also be used if the intention to pose a face is irrelevant from the awareness of its use in an academic study.
Current lives are closely twisted with time spent on social networks (SNs), where users express their emotions by posting facial photos towards life episodes [70,71].An SN user can choose to post photos in accordance with geographical information at the check-in location.This activates a probability for further analysis by matching data derived from the photo and the real-time location [63,64].As a type of SN user, employees also have a habit to use twitter micro-blogs on SNs to expose their sentiments towards working affairs [72].Hence, a big number of workers' facial photos have been pooled on SNs, which can be used to quantify sentiments therein, estimate salary satisfaction, and perform regression against the landscape metrics of GBSs at workplaces.To the best of our knowledge, it is very rare to document a study that analyzes all these responses and detects their relationships.
In this study, east China was selected as the study area, where a total of 56 industrial parks were chosen as study plots.Salary satisfaction was evaluated using model regression against facial expression scores.The regional GBS landscape was analyzed remotely using the digital data of satellite imageries, and GVI and plant diversity were analyzed using data crawling, deep learning, and artificial recognition.Our objective was to model salary satisfaction using the rated facial expression scores of employees.We also aimed to detect regressed levels of facial expression scores and salary satisfaction against the Land 2023, 12, 2075 4 of 24 driving forces of combined GBS landscape metrics, GVI, and visual plant diversity using regressing models.

Study Design and Layout
The study layout and explanatory and dependent variables are all shown in Figure 1.Industrial parks with business corporations were selected as the study plots across an area of east China.Landscape metrics for GBSs in industrial parks were chosen as the explanatory variables, which were evaluated remotely using satellite imageries.The spatial range of every industrial park was framed to a buffer area, wherein landscape metrics were evaluated for green space, blue space, and the general buffer zone containing all surface features.Plant diversity was recognized artificially by a professional group of botanical experts for evaluating GVI.Dependent variables were chosen as exposed emotions assessed through facial expression scores and recruitment wage.Facial photos were obtained from a social network platform at the locations of the objective industrial parks on working days.Facial expressions were recognized and rated with emotional scores.Hence, the recruitment wage was modelled against the facial expression scores, which resulted in evaluations of salary satisfaction.The estimated values of salary satisfaction were further regressed against all the landscape metrics as explanatory variables.
Land 2023, 12, x FOR PEER REVIEW 4 of 24 using data crawling, deep learning, and artificial recognition.Our objective was to model salary satisfaction using the rated facial expression scores of employees.We also aimed to detect regressed levels of facial expression scores and salary satisfaction against the driving forces of combined GBS landscape metrics, GVI, and visual plant diversity using regressing models.

Study Design and Layout
The study layout and explanatory and dependent variables are all shown in Figure 1.Industrial parks with business corporations were selected as the study plots across an area of east China.Landscape metrics for GBSs in industrial parks were chosen as the explanatory variables, which were evaluated remotely using satellite imageries.The spatial range of every industrial park was framed to a buffer area, wherein landscape metrics were evaluated for green space, blue space, and the general buffer zone containing all surface features.Plant diversity was recognized artificially by a professional group of botanical experts for evaluating GVI.Dependent variables were chosen as exposed emotions assessed through facial expression scores and recruitment wage.Facial photos were obtained from a social network platform at the locations of the objective industrial parks on working days.Facial expressions were recognized and rated with emotional scores.Hence, the recruitment wage was modelled against the facial expression scores, which resulted in evaluations of salary satisfaction.The estimated values of salary satisfaction were further regressed against all the landscape metrics as explanatory variables.

Study Area and Sampling Plots
East China was chosen as the study area, where a total of 56 industrial parks were selected from 16 provincial regions (3 municipalities of Beijing, Tianjin, and Shanghai and 13 provinces of Ningxia, Shanxi, Hebei, Heilongjiang, Jilin, Liaoning, Shandong, Jiangsu, Zhejiang, Fujian, Guangdong, Henan, and Shaanxi).Figure 2 shows the spatial distributions of the locations for the 56 industrial parks.The specific names and coordinates are shown in Table 1.Cities in east China have experienced fast urbanization in recent

Study Area and Sampling Plots
East China was chosen as the study area, where a total of 56 industrial parks were selected from 16 provincial regions (3 municipalities of Beijing, Tianjin, and Shanghai and 13 provinces of Ningxia, Shanxi, Hebei, Heilongjiang, Jilin, Liaoning, Shandong, Jiangsu, Zhejiang, Fujian, Guangdong, Henan, and Shaanxi).Figure 2 shows the spatial distributions of the locations for the 56 industrial parks.The specific names and coordinates are shown in Table 1.Cities in east China have experienced fast urbanization in recent decades, and they account for the major body of gross domestic product (GDP) in China [44,45].According to Meng et al. (2022) [2], most of the chosen industrial parks made a big contribution to regional heat island effects.Local employees may perceive remarkable positive emotions due to the cooling effect [15].The chosen industrial parks stationed businesses with ecommercial retails, mainly in high-tech and big-scale manufacture industrials.Businesses in these industrial parks also recruited positions with high levels of wage in their host regions.
Land 2023, 12, x FOR PEER REVIEW 5 of 24 decades, and they account for the major body of gross domestic product (GDP) in China [44,45].According to Meng et al. (2022) [2], most of the chosen industrial parks made a big contribution to regional heat island effects.Local employees may perceive remarkable positive emotions due to the cooling effect [15].The chosen industrial parks stationed businesses with e-commercial retails, mainly in high-tech and big-scale manufacture industrials.Businesses in these industrial parks also recruited positions with high levels of wage in their host regions.Jinan Int.Biol. 6Med. 7Park 117.129 36.68093

Facial Images and Emotional Scores
The ethical norm for employing facial photos of human subjects was reviewed and approved by the Ethic Review Committee of Ice and Snow Tourism Resorts Equipment and Intelligent Service Technology, Ministry of Culture and Tourism of The People's Republic of China.All photos were obtained from open twitters on social networks, which were posted to the public with full awareness of their users.Human-related data were obtained and analyzed only as numeric values that were used in this study without any consent or permission to be released to any third parties, to avoid the potential release of biological information attached to the subjects' faces and their personal privacy.
A rule was set for selecting industrial parks that the number of employees' facial photos had to be at least 100 during working times in the days from 1 January 2020 to 31 July 2023.Thus, the period for photo searching was as long as 3.58 years, which is much longer than most of the current studies collecting photos from urban parks [19,21,24,45].This procedure was conducted due to three reasons.Firstly, photos taken by people at workplaces during working times were much fewer than those taken in open urban parks, and the collection period was long enough that the number of photos could be met.Secondly, we aimed to conduct a national investigation that covered a geographical range of most provinces or municipalities in east China, and the number of photos for a single industrial park should be more than 100 according to current evidence [18,45,48].Any shorter time of photo collection failed to fulfill the total number of photos.Finally, most times during 2020-2022 fell in episodes of the COVID-19 pandemic, which reduced the number of workers in industrial parks due to lockdowns.
Sina Weibo (or namely Chinese twitter) [73] and Douyin (also being aware of Chinese Tik Tok) [74] were used as social networks for photo collection.Sina Weibo is a twittering platform where photos could be obtained from the micro-blogs of users with the time and place exposed when the photos were posted [63,64].Douyin is an online platform where users share short videos on their private pages but open updates where date and place are also available (if no privacy ban was set by the user).Micro-videos downloaded from Douyin were evenly cut to 10-20 sliced photos according to the total time used.In this study, all targeted photos that were documented had to follow all of the following rules: (1) All photos had to be taken at the places in any IP listed in Table 1 and the time had to be from 08:00 a.m. to 21:00 p.m. from Monday to Friday every week, excluding holidays and sudden events.(2) A photo had to contain at least one intact human face with all sensing organs fully exposed.
(3) Any photos with extremely high and low ages of people, who looked too young (e.g., toddlers and infants) or too old (e.g., senior citizens with visible disabilities), were excluded, because these people were unlikely workers.
The photos were firstly cropped to make sure only a face of one person left therein, which accounted for about 75% of the whole projected area.Thereafter, the photos were rotated to make the nose axis vertical to the horizontal line to ensure a high recognition rate by the software [45,66,67].All the prepared photos (initial n = 5600, 100 photos per IP) were sent to be analyzed for facial expression scores.The FireFACE ver 1.0 software (Hainan Guanzhong A&F Inc., Sanya, China) was used to rate happy, sad, and neutral emotional scores with the positive response index (i.e., happy score minus sad score) calculated.This instrument has been validated to identify that its recognition accuracies for happy and sad emotional scores can be fully accepted in comparison to all similar instruments [75].The FireFACE ver 1.0 software has been successfully used in studies testing the facial expression scores of visitors experiencing GBSs in urban parks at multiple geographical scales [21,24,44,45].
Finally, a total of 973 photos failed to be recognized by the software, leaving another 4627 photos successfully rated for emotional scores.The recovery rate of available recognitions was 82.63%.As it is shown in Table 1, the generated dataset contained a data pool of photo numbers to range within 40-98 individuals per industrial park with an average of 82.63 ± 8.68 (±standard error) individuals and a coefficient of variance (CV) of 0.11.Even the minimum number of photos per industrial park was similar or higher than that in other studies (>40) at regional scales [20,48,49,76].Overall, our facial expression scores can be accepted for further analysis.

Remote Evaluation of Green and Blue Space Landscape Metrics
All the industrial parks were outlined in ArcGIS software (version 10.2, Esri-Branch, Shanghai, China).Shanghai International Tourism Park had the largest radius of an approximate square (~2 km) and Zhejiang Xiangyuan Cultural Tourism Park had the lowest (~30 m).A buffer area was set for every industrial park area with a circling buffer in a radius of 2 km, as this is used as a common geographical criterion for outlining GBS areas (Figure 3).approximate square (~2 km) and Zhejiang Xiangyuan Cultural Tourism Park had the lowest (~30 m).A buffer area was set for every industrial park area with a circling buffer in a radius of 2 km, as this is used as a common geographical criterion for outlining GBS areas (Figure 3).Landsat 8 OLI satellite imageries (30 m × 30 m) were used for the GBS landscape analysis.Green space area was evaluated as the largeness of the normalized difference vegetation index (NDVI) that was calculated as [77,78]: where Bandni and Bandred are surface reflections at the near-infrared and red bands, respectively.Blue space area was evaluated as the largeness of the normalized difference water index (NDWI) that was calculated as [79]: where Bandgreen is the surface reflection at the green band.Elevation was evaluated using the data of the digital elevational model (DEM) from the Aster GDEM 30M map [80].The height of the tallest surface feature was evaluated using the data of the digital surface model (DSM) using the AW3D30 DSM map [81].The heights of the mapped surface Landsat 8 OLI satellite imageries (30 m × 30 m) were used for the GBS landscape analysis.Green space area was evaluated as the largeness of the normalized difference vegetation index (NDVI) that was calculated as [77,78]: where Band ni and Band red are surface reflections at the near-infrared and red bands, respectively.Blue space area was evaluated as the largeness of the normalized difference water index (NDWI) that was calculated as [79]: where Band green is the surface reflection at the green band.Elevation was evaluated using the data of the digital elevational model (DEM) from the Aster GDEM 30M map [80].The height of the tallest surface feature was evaluated using the data of the digital surface model (DSM) using the AW3D30 DSM map [81].The heights of the mapped surface features were evaluated from the averages of a layer superimposing the DSM data onto the DEM data [82].

Green View Index and Plant Diversity Evaluation
A total of 100 street view photos were obtained from the buffer area of an industrial park.Street view images were the source for evaluating GVI, which were obtained from the Baidu application programming interface (API) [83].As described above, the buffer areas had a radius of 2 km, resulting in an area of 12.56 km 2 .There was a probability of ~25% for successfully matching a coordinate and a photographed image in the Baidu map street-view [84].At least 400 grids were set in a buffer to expect that a hundred out of them could be crawled through API to match a host location with existing street view images.Therefore, a grid had an area of 3.14 ha, wherein a street view image was expected to be obtained from coordinates.
Deep learning was developed using Python 3.6 (Python Software Foundation, Python, Beaverton, ON, USA).We employed deep learning to train the machine recognizing pixels occupied by green plants (leaves, needles, twigs, and crowns, etc.) based on an initial dataset of cityscapes [85].A total of 2000 street view images were downloaded randomly by coordinates in grids with a radius of 3.5 km 2 arranged across the municipal area of the host city where an objective IP was constructed.A convolutional neural network was used to train the machine, recognizing the semantic tags of objectives in a street image.The DeepLabV3+ model was employed for dividing objectives through recognizing and categorizing every pixel by the major network from the residual network-101 (ResNet-101).The trained program was coded to an application (Pycharm, Praha, Czech Republic), which was validated to recognize objective pixels as a part of green organs in a plant with an accuracy over 85% [75].All the GVI values estimated from a buffer area were averaged for the host industrial park.Thereafter, the application was run for recognizing GVI in every street view image, which ranged from 0.00% to 87.94% with a median of about 11.00% (Figure 4).
Land 2023, 12, x FOR PEER REVIEW 9 of 24 features were evaluated from the averages of a layer superimposing the DSM data onto the DEM data [82].

Green View Index and Plant Diversity Evaluation
A total of 100 street view photos were obtained from the buffer area of an industrial park.Street view images were the source for evaluating GVI, which were obtained from the Baidu application programming interface (API) [83].As described above, the buffer areas had a radius of 2 km, resulting in an area of 12.56 km 2 .There was a probability of ~25% for successfully matching a coordinate and a photographed image in the Baidu map street-view [84].At least 400 grids were set in a buffer to expect that a hundred out of them could be crawled through API to match a host location with existing street view images.Therefore, a grid had an area of 3.14 ha, wherein a street view image was expected to be obtained from coordinates.
Deep learning was developed using Python 3.6 (Python Software Foundation, Python, Beaverton, ON, USA).We employed deep learning to train the machine recognizing pixels occupied by green plants (leaves, needles, twigs, and crowns, etc.) based on an initial dataset of cityscapes [85].A total of 2000 street view images were downloaded randomly by coordinates in grids with a radius of 3.5 km 2 arranged across the municipal area of the host city where an objective IP was constructed.A convolutional neural network was used to train the machine, recognizing the semantic tags of objectives in a street image.The DeepLabV3+ model was employed for dividing objectives through recognizing and categorizing every pixel by the major network from the residual network-101 (Res-Net-101).The trained program was coded to an application (Pycharm, Praha, Czech Republic), which was validated to recognize objective pixels as a part of green organs in a plant with an accuracy over 85% [75].All the GVI values estimated from a buffer area were averaged for the host industrial park.Thereafter, the application was run for recognizing GVI in every street view image, which ranged from 0.00% to 87.94% with a median of about 11.00% (Figure 4).All the images that were successfully analyzed for GVI were sent to a professional botanical group with one teacher and four undergraduates.All the plant species were recognized visually by two of the undergraduates, and after the two agreed with each other for recognizing all species, the teacher made the final determination.Due to the limit of resolution, only trees and shrubs were recognized, and these two plant types are also major botanical drivers that can cause visitors' attention and result in emotional perceptions [24].The general level of GVI varied to a large extent (Figure 4), leading to some images containing no plants while others contained a rich amount.Finally, a total of 5250 photos were successfully recognized to contain at least one plant species (recovery: 93.75%).Plant species diversity was quantified through two equations, namely the Simpson's diversity index (Simpson index) and the Shannon-Wiener diversity index (Shannon index) [24]: where p i is the proportion of species i relative to the total number of all species up to the final species S and lnp i is the natural logarithm of the proportion p i of species i.Again, all variables assessing plant diversity were averaged as a mean for the host IP.

Recruitment Wage and Salary Satisfaction
In this study, we put forth a simple but efficient equation for quantifying the salary satisfaction (SS) of employees in industrial parks: where ∂ is the estimate of the maximum likelihood parameter that needs to be estimated; Var(y i ) is the variance model of emotional response y i ; expression scores for happy, sad, and neutral emotions can be assigned with i = 1, 2, and 3, respectively, assuming a visible face harboring salary satisfaction can be assessed based on a mixture of scores for three basic component emotions; and I is the intercept value in the model for analyzing the maximum likelihood of parameters.A generalized linear model (GLM) was defined as classes from an extension of traditional linear models that allow for the mean of a population to depend on a linear predictor through a nonlinear link function and allows the response probability distribution to be any member of an exponential family of distributions [86].According to histogram analyses (Figure S1), happy scores showed a normal distribution pattern, and sad scores and recruitment wage (R W ) both showed approximate patterns with positive values in skewness by about ~35%.Hence, a normal regression model was employed from the GLM family to estimate salary expectation against perceived happy, sad, and neutral emotional scores [87]: where the variance of the response y i [Var(y i )] depends on the mean µ through a variance function V; Ø is a dispersion parameter which needs to be estimated; w i is a known weight for each observation; and µ i is the expected value of y i , which are formulized as: where the quantity x i is a column vector of the covariates for the observation i that is the number of explanatory variables; the vector of unknown coefficients β is estimated by at least squares fit to the data y; and ε i is assumed to be an independent and normal random variable with zero mean and constant variance.

Data Analysis and Statistics
SAS software (ver.9.4, SAS Statistics Institute, Cary, NC, USA) was used for the data analysis and statistics.Data from explanatory variables (landscape metrics, GVI, and plant diversity indexes) were shown as spatial distribution patterns.R W and salary satisfaction were grouped to show their distribution patterns.Emotional scores (happy, sad, neutral, PRI) were ranked to make a distribution-free dataset [88], which can be regressed without statistical limits using general linear modeling (GLM) models [19,21,44,67].Spearman correlation was used to detect relationships between pairs of facial expression scores.The satisfied wage was estimated with the maximum likelihood parameter method (using Genmod procedure in SAS) using the normal-likely algorithm converged in the GLM model.Estimates could be accepted for further correlation only when the Wald Chi-square (ChiSq) was large enough to make sure the Pr (>ChiSq) was lower than 0.05.Linear correlation was used to bridge the relationship between emotional scores and R W for scattered significant estimates, with all negative estimates excluded to meet the objective logistic recognition (no one would expect a wage with a negative amount).The multivariate linear regression (MLR) model was employed using a backward regressing manner, with all variables entered and eliminated until significant parameter estimates were left (p < 0.05).Dependent variables were chosen as the four emotional scores and salary satisfaction.

Spatial Distributions of Landscape Metrics
As it is shown in Figure 5, the industrial parks constructed in Hebei, Jiangsu, Shanghai, and Guangdong were generally located at lower elevations than those in other regions (Figure 5A).The surface features were generally taller in buffer areas of industrial parks located in northeast China (Heilongjiang, Jilin, and Liaoning), north China (Beijing, Henan, Shanxi, and Shaanxi), and parts of Zhejiang and Guangdong (Figure 5B) than those in the other industrial parks.The feature height was lower in buffer areas of industrial parks in Jiangsu, Shanghai, and parts of Guangdong than that in other regions.Green space area mainly ranged within 350-800 m 2 per industrial park, with obvious smaller ones located in Heilongjiang, Beijing, Shanxi, Shandong, Jiangsu, and Guangdong (Figure 5C).Green space height was generally lower in south China such as in regions of Jiangsu, Shanghai, Fujian, and Guangdong (Figure 5D).Blue space area was generally larger in southern parts of China than in the north, with occasional large areas of about ~200 m 2 per industrial park in Ningxia, Shanxi, and Tianjin (Figure 5E).

Spatial Distributions of GVI and Plant Diversity
GVI ranged within 0.16-0.20 with occasionally lower values in buffer areas of industrial parks concentrated in parts of Jiangsu, Shanghai, and Guangdong (Figure 6A).Both the Simpson and Shannon indexes showed alternative high and low levels along the northsouth geographical gradient (Figure 6B,C).The buffer areas of industrial parks had lower diversity indexes in southern Liaoning, western Shandong, Fujian, and parts of Guangdong.These regions also contained industrial parks whose buffer areas were found to have fewer species amounts (Figure 6D).

Spatial Distributions of Facial Express Scores
Employees showed faces with high happy scores (up to ~50%) in industrial parks in three general agglomerations in north China (Hebei, Beijing, and Tianjin), east China (Jiangsu, Shanghai, and Zhejiang), and south China (Guangdong) (Figure 7A).In contrast, these regions were also more common with low sad scores (Figure 7B).Extremely high sad scores (up to ~20%) were found in industrial parks located in Shaanxi, Henan, Liaoning, and eastern Zhejiang.Neutral scores showed a similar distribution pattern with sad scores (Figure 7C).Again, high PRI scores were found to agglomerate in industrial parks around Beijing, Shanghai, and Guangdong (Figure 7D).

Spatial Distributions of GVI and Plant Diversity
GVI ranged within 0.16-0.20 with occasionally lower values in buffer areas of industrial parks concentrated in parts of Jiangsu, Shanghai, and Guangdong (Figure 6A).Both the Simpson and Shannon indexes showed alternative high and low levels along the

Spatial Distributions of Facial Express Scores
Employees showed faces with high happy scores (up to ~50%) in industrial parks in three general agglomerations in north China (Hebei, Beijing, and Tianjin), east China (Jiangsu, Shanghai, and Zhejiang), and south China (Guangdong) (Figure 7A).In contrast, these regions were also more common with low sad scores (Figure 7B).Extremely high sad scores (up to ~20%) were found in industrial parks located in Shaanxi, Henan, Liaoning, and eastern Zhejiang.Neutral scores showed a similar distribution pattern with sad scores (Figure 7C).Again, high PRI scores were found to agglomerate in industrial parks around Beijing, Shanghai, and Guangdong (Figure 7D).
Spearman correlation indicated that happy score was negatively related to sad and neutral scores, but had a positive relationship with PRI score (Table 2).Sad score had a positive relationship with neutral score, but a negative relationship with PRI score.Neutral score also had a negative relationship with PRI score.Spearman correlation indicated that happy score was negatively related to sad and neutral scores, but had a positive relationship with PRI score (Table 2).Sad score had a positive relationship with neutral score, but a negative relationship with PRI score.Neutral score also had a negative relationship with PRI score.

Salary Satisfaction Estimate and Spatial Distribution
According to GLM regression, the I' for salary satisfaction was estimated to be RMB 8000 with Wald 95% confidence limits within 7998.04-8001.96.Facial expression scores for happy, sad, and neutral emotions were all input to GLM, regressing recruitment wage, and only happy scores were retained to obtain a series of significant estimates of the maximum likelihood parameters.A total of 53 parameters were estimated to be significant, but 18 of them were estimated to be negative values, which had to be excluded from further analysis under the rule following a general logistic recognition.As a result, salary satisfaction can be estimated as: where x Happy is the input variable of the observed happy scores that resulted in positive estimates (n = 35).The coefficient determination was 46.76% with the probability of significance of 0.0046, both indicating a highly linear distribution of data.
As it is shown in Figure 8, the published recruitment wage was higher in industrial parks located in agglomerated regions around Beijing, Shanghai, and Guangdong (Figure 8A).However, estimated salary satisfaction showed a distribution pattern with a much higher homogeneity (Figure 8B).Most salary satisfaction values ranged around RMB 19,000 and some lower values around RMB 9000 distributed in industrial parks in northern regions and eastern Zhejiang and Fujian.
National 863 Central Software Park was the only industrial park where corporates paid salaries as high as the satisfactory levels (RMB 742.12) in central China (in Zhengzhou Henan) (Figure 8C).In Shanghai, businesses in three industrial parks paid satisfactory salaries, i.e., Jiading Industry Development Park (RMB 38.43), Shanghai Chemical Industry Park (RMB 7875.70), and Shanghai Fine Chemical Engineering Industry Park (RMB 5359.31).Only one industrial park's corporates provided satisfactory salaries (RMB 1645.31) in Jiangsu, which was Jiangsu Tianmu Lake Tourism Joint Stock Company Park.Two industrial parks were found with satisfactory salaries in Zhejiang, which were Zhejiang Xiangyuan Cultural Tourism Inc. Park (RMB 5045.93) and Zhejiang Haizheng Medine Industry Inc. Park (RMB 3133.19).Sihui Fine Chemical Industry Park was the only industrial park with a satisfactory salary (RMB 2237.69) in Guangdong.
National 863 Central Software Park was the only industrial park where corporates paid salaries as high as the satisfactory levels (RMB 742.12) in central China (in Zhengzhou Henan) (Figure 8C).In Shanghai, businesses in three industrial parks paid satisfactory salaries, i.

Estimation of Salary Satisfaction
Employee satisfaction is important for organizational well-being and responsible for decreased turnover [54].Expectations for a higher salary are associated with career satisfaction, which was reported by international enterprise professionals [51], subway workers [56], humanitarian workers [52], and clinical social workers [54,89].Satisfaction about salary was also reported to be associated with job satisfaction by public hospital doctors [90], treated school teachers [69], and market access professionals [68].These types of stud-

Estimation of Salary Satisfaction
Employee satisfaction is important for organizational well-being and responsible for decreased turnover [54].Expectations for a higher salary are associated with career satisfaction, which was reported by international enterprise professionals [51], subway workers [56], humanitarian workers [52], and clinical social workers [54,89].Satisfaction about salary was also reported to be associated with job satisfaction by public hospital doctors [90], treated school teachers [69], and market access professionals [68].These types of studies mainly consider that demographical effects are the major driving force that impose impacts employees' satisfactions on salaries [68,69].Other explanatory factors have occasionally been reported from aspects in skillful fields, e.g., years in practice, education background, and specialty for hospital physicians [90].The current estimation of salary satisfaction is limited to workers from a narrow range of occupational types.To break this restriction, Sanchez-Gomez, Breso, and Giorgi [57] tried to put forward a model to predict salary gain against quantified emotional intelligence (EI).They found that combined increases in both EI and emotional repair capacity together attributed to a higher level of salary.This study initiated the probability to predict salary satisfaction by modeling the emotional performance of an employee.However, technical limits in emotional data collection existed in this study and in those without novel instruments.
Our study employed a novel methodology using inter-disciplinary pieces of technique and knowledge for estimating salary satisfaction.In this study, emotional scores were employed as explanatory variables that were used for modeling salary satisfaction.About this, one may argue that salary satisfaction is perceived in response to recruitment wage, hence, the estimation against emotional scores did not have any applicative meaning.This argument fails to consider that salary satisfaction cannot be directly reflective following awareness of recruitment wage, unless through perceiving various emotions.That is why our salary satisfaction was estimated against all types of emotional scores.We also considered the condition that significant estimates may fluctuate to be either positive or negative values, which had to be screened to remove negative values to meet the common sense before being modeled by another linear correlation.This process was achieved using statistics that other studies took into account for estimating salary satisfaction value [91,92].Therefore, our estimation of salary satisfaction was accurate at least to a greater extent than others in the current methodology.

Facial Expression Scores of Employees in Industrial Parks
In our study, the emotional states of employees were quantified through facial photos that were posted by social network users.This type of emotional data may suffer queries about unavailability to reflect real-time subtle emotions, because facial emotions are mostly posed with prepared or intended sentiments [71].We argue that this is not the issue that determines the accuracy of our results.All the photos were collected from twitters posted by people exposed to GBSs at workplaces, which has been identified to be a moderate psychological stimulation that can evoke sentiment response through perceiving changes in landscape metrics [18], regional meteorology [19,48], and atmospheric pollution [44,45].The magnitude to which an employee posed emotions was elicited by the same source of psychological hints, even though they knew what emotions they intended to pose.Secondly, the intended emotions were posed to a camera, but subjects can rarely be aware that their faces will be taken as a source in an academic study.That is, the intention to pose sentiments had no relevance to any existing dataset used for testing facial emotions.Hence, the results will not be affected by the subjects' intentions, no matter whether they meant to do so or not.Finally, we controlled our subjects to be a group of employees who were targeted at working times in workplaces, which can mostly control the systematic errors in big social network data caused by individuals' intentions.Overall, it is acceptable to employ facial photos as a source of emotional data from social networks.
Our facial expression scores showed a spatial pattern that positive sentiments were exposed in agglomerated regions around Beijing, Shanghai, and Guangzhou.These three core regions are also three biggest metropolises of China, with huge residential populations (>15 million) and large areas of built-up lands.According to findings unraveled in urban parks [21,63], places with crowded populations on highly urbanized lands tend to evoke more negative sentiments of visitors, which disagree with our results.However, other national-scale studies have revealed that people would be likely to pose more smiles or positive emotions in municipal regions where the total product was high in retailing sales [44,45].Activated retail records precondition high incomes for local households, hence, people tended to look happier than those from regions with smaller retailing markets.
Because of the negative relationship between happy and sad scores, their spatial distributions also showed contrasting patterns.Similar spatial patterns between happy and PRI scores suggest that employees generally showed more positive emotions (20-60%) than negative ones (10-25%).Although neutral scores also had a negative relationship with happy scores, their spatial distribution patterns were not different to the extent as those between happy and sad scores.We found that, in three agglomerated regions with high positive sentiments, neutral scores were exposed in industrial parks alternatively with high happy scores.These suggest that employees tended to post neutral emotions when they felt not so happy to make a smile.Regarding that neutral scores were positively correlated with sad scores, parts of exposed neutralities may be an exterior performance of sadness.

Driving Forces When Being Exposed to Green Spaces
It is surprising to find that exposure to green space caused totally negative emotional consequences perceived by employees.Plant diversity quantified by the Simpson index was perceived as a trigger of neutral emotions or parts of drivers that depressed positive emotions.These are totally against previous findings in the urban parks of China [24] and private gardens of Germany and New Zealand [93].Our findings also disagree with the nature of ART and SRT, and their former theory of Biophilia hypothesis.We surmise that this contradiction results from the different ways to experience exposure to green space.People visit parks and gardens without necessary stressful burdens, hence, their attention are easily attached to the focus on plant biodiversity.In contrast, employees were concerned with affairs at working time when they perceived heavy mental stresses deeply in their work.Deeply stressed cohorts of people may not find it easy to recover with an experience with nature.The diversity of plant species is far away from workers' concerned interests, hence they showed more neutral faces.A view with highly diverse species may make workers consider a hint of attention transfer, which was perceived as a negative contribution to exposed happiness.
The elevation of the location of the industrial park was found to have a contribution to exposing sadness, which disagrees with findings for visitors of tropical forest parks with degraded ecological states [19].Again, a visit to an urban park is unlikely or rarely attached to the necessity of preconditioned mental burdens; hence, a high elevation in a subtropical mountain would be accompanied with higher wind velocities that bring thermal comfortable feelings.However, the high elevation of an industrial park can make workers notice more harsh conditions for everyday accessibility to the workplace.This surmise was reasonable, because industrial parks located at high elevations (200-1000 m a.s.l.) were mostly located in montane regions whose host cities were subjected to a terrain with moderate to low elevations, such as industrial parks located in Heilongjiang, Shanxi, Shaanxi, and Ningxia.
Both GVI and green space height both contributed to the effects that caused an upregulation of exposed sadness, suggesting that employees disliked landscapes with a high ratio of seen greenness provided by groups of tall plants at their workplaces.It has been found that the height of plants accounts for a large proportion of GVI in urban green space [94].Although tall plants are responsible for a high GVI, they may be perceived as a trigger of negative emotions due to their intervention of expected openness but upregulated enclosure in workers' views.According to He et al. [95], an urban view with attributes of low openness but a high enclosure can be a trigger of unstable senses.Green space with tall plants also means a larger understory space where more shrub diversity can be easily seen [24], which is also has a depressing effect on exposed smiles.
The area of green space was perceived as a depression against posted smiles, hence performed as a sole facet that reduced salary satisfaction.This totally fails to concur with previous findings that have demonstrated a positive relationship between green space largeness and positive emotions [16][17][18]20].Continuous to the above discussed contents, a large area of green space is associated with a higher chance of experiencing more plant species, which did not benefit the perceptions of positive emotions as joint effects with plant diversity.In industrial parks where the areas for workplaces were fixed to a geographical range in a radius of ≤2 km, larger areas of green spaces meant a longer distance of path between constructions which impacted accessibility to the workplace.

Applicative Meaning of This Study
The whole process and way that we estimated salary satisfaction can be referred to by further studies that aim to assess salary satisfaction against employees' emotional performances.The synthesis of data collection through remote sensing, online crawling, and deep learning can be strengthened by more skillful manipulations to harvest a more accurate dataset.
Among the 56 industrial parks, 8 of them were found to have corporates that could provide salaries meeting employees' satisfactions.Among them, three were in chemical industry and two were in tourism; hence, these two types of industries were gaining decent profits and had capacities to recruit with high wages.The highest difference between recruitment wage and salary satisfaction (RMB 7875.70) was also provided by an industrial park engaging in the chemical industry, namely the Shanghai Chemical Industry Park.The other chemical industry park in Shanghai (Shanghai Fine Chemical Engineering Industry Park) also provided salaries with a high recovery from the difference of recruitment wage and salary satisfaction (RMB 5359.31).Hence, the chemical industry was identified as paying higher salaries for employees than other industries.On the other hand, among the eight industrial parks, six were located in regions around Shanghai, including those in Zhejiang and Jiangsu, which suggest corporations in industrial parks of this area were certified to provide high salaries.These findings can be referred to by coming investments in industrial park construction in China.
Our results about the driving forces of landscape metrics and plant diversity corrected a common recognition that more green space may be better at the workplace.According to our findings, an industrial park should be constructed with a small green space decorated in dwarf plants with a low number of species.This would avoid the impairment of positive sentiments in workers and reduce their perceptions of low satisfaction about monthly salaries.This manner of industrial park landscape construction will be referred to by labor management for controlling wage cost through controlling high-income desires.

Limits of This Study
Our study has three limitations: Firstly, the time of the study fell in the COVID-19 pandemic, which modified the real-time situation at workplaces and changed workers' perceptions with uncertain consequences.Further work should be conducted in ordinary times to avoid unknown interventions on results by sudden social hotspot events.
Secondly, our methodology to estimate salary satisfaction can be validated in cooperation with self-reported expectations.It is a novelty for our study to publish the new methodology and it is obligated for following studies to validate and optimize the methods.
Finally, it will be better if all industrial parks were controlled to have large walled areas to fully make sure exposure to GBSs occurred at workplaces.

Figure 1 .
Figure 1.Layout and design the study.Abbreviations: DSM, digital surface model; DEM, digital elevation model; NDVI, normalized difference vegetation index; and NDWI, normalized difference water index.

Figure 1 .
Figure 1.Layout and design the study.Abbreviations: DSM, digital surface model; DEM, digital elevation model; NDVI, normalized difference vegetation index; and NDWI, normalized difference water index.

Figure 2 .
Figure 2. Spatial distributions of industrial parks across a geographical range of East China.

Figure 2 .
Figure 2. Spatial distributions of industrial parks across a geographical range of East China.

Figure 3 .
Figure 3.The illustration of buffer areas in a circle with a radius of 2 km for industrial parks with the largest (Shanghai Int.Tour.Park), medium (Taiyuan Qinghua S&T Park), and minimum (Zhejiang Xiangyuan Cul.Tour.Park) areas.

Figure 3 .
Figure 3.The illustration of buffer areas in a circle with a radius of 2 km for industrial parks with the largest (Shanghai Int.Tour.Park), medium (Taiyuan Qinghua S&T Park), and minimum (Zhejiang Xiangyuan Cul.Tour.Park) areas.

Figure 4 .
Figure 4. Typical views in Baidu street view images with green view indexes (GVIs) variation from 0.00% (A) through a median of 11.09% (B) to 87.94% (C).All images are randomly crawled at coordinates of grids (3.14 ha) from a 2 km radius buffer area surrounding an industrial park in China.

Figure 5 .
Figure 5. Spatial distributions of elevation (DEM) (A) and surface feature height (B), area (C) and height of green spaces (D), and blue space area (E) in buffer areas of industrial parks of East China.

Figure 5 .
Figure 5. Spatial distributions of elevation (DEM) (A) and surface feature height (B), area (C) and height of green spaces (D), and blue space area (E) in buffer areas of industrial parks of East China.

Figure 6 .
Figure 6.Spatial distributions of green view index (GVI) (A), Simpson diversity index (B), Shannon Wiener diversity index (C), and plant species amount (D) in buffer areas of industrial parks of East China.

Figure 6 . 24 Figure 7 .
Figure 6.Spatial distributions of green view index (GVI) (A), Simpson diversity index (B), Shannon Wiener diversity index (C), and plant species amount (D) in buffer areas of industrial parks of East China.Land 2023, 12, x FOR PEER REVIEW 14 of 24

Figure 7 .
Figure 7. Spatial distributions of facial expression scores for happy (A), sad (B), neutral (C), and PRI scores (D) of employees in industrial parks of east China.
e., Jiading Industry Development Park (RMB 38.43), Shanghai Chemical Industry Park (RMB 7875.70), and Shanghai Fine Chemical Engineering Industry Park (RMB 5359.31).Only one industrial park's corporates provided satisfactory salaries (RMB 1645.31) in Jiangsu, which was Jiangsu Tianmu Lake Tourism Joint Stock Company Park.Two industrial parks were found with satisfactory salaries in Zhejiang, which were Zhejiang Xiangyuan Cultural Tourism Inc. Park (RMB 5045.93) and Zhejiang Haizheng Medine Industry Inc. Park (RMB 3133.19).Sihui Fine Chemical Industry Park was the only industrial park with a satisfactory salary (RMB 2237.69) in Guangdong.

Figure 8 .
Figure 8. Spatial distributions of values for recruitment wage (Rw) (A), salary satisfaction (SS) (B), and their difference (Rw-SS) (C) in industrial parks of East China.

Figure 8 .
Figure 8. Spatial distributions of values for recruitment wage (Rw) (A), salary satisfaction (SS) (B), and their difference (Rw-SS) (C) in industrial parks of East China.

Figure 9 .
Figure 9. Multivariate linear regression (MLR) models of ranked scores for happy (A), sad (B), and neutral (C) emotions as well as PRI (D) for employees against combined explanatory factors of elevation (Buffer DEM), surface feature height in green space (Green space height), green space area, GVI, and Simpson index in buffer areas of industrial parks.Dots present significant parameter estimates (p < 0.05) for independent factors with error bars marking standard errors.

Figure 9 .
Figure 9. Multivariate linear regression (MLR) models of ranked scores for happy (A), sad (B), and neutral (C) emotions as well as PRI (D) for employees against combined explanatory factors of elevation (Buffer DEM), surface feature height in green space (Green space height), green space area, GVI, and Simpson index in buffer areas of industrial parks.Dots present significant parameter estimates (p < 0.05) for independent factors with error bars marking standard errors.

Land 2023 , 24 Figure 10 .
Figure 10.The MLR model of salary satisfaction for employees against combined explanatory factors of elevation (Buffer DEM), surface feature height in green space (Green space height), green space area, GVI, and Simpson index in buffer areas of industrial parks.Dots present significant parameter estimates (p < 0.05) for independent factors with error bars marking standard errors.

Figure 10 .
Figure 10.The MLR model of salary satisfaction for employees against combined explanatory factors of elevation (Buffer DEM), surface feature height in green space (Green space height), green space area, GVI, and Simpson index in buffer areas of industrial parks.Dots present significant parameter estimates (p < 0.05) for independent factors with error bars marking standard errors.

Table 1 .
List of industrial parks with names, coordinates, and numbers of posed employees' photos.

Table 1 .
List of industrial parks with names, coordinates, and numbers of posed employees' photos.Order Industrial Park Name Longitude ( • ) Latitude ( • ) Photo No.

Table 2 .
Correlations between pairs of facial expression scores.

Table 2 .
Correlations between pairs of facial expression scores.