Research Article | Open Access
Konstantinos Koufos, Riku Jäntti, "Proportional Fair Power Allocation for Secondary Transmitters in the TV White Space", Journal of Electrical and Computer Engineering, vol. 2013, Article ID 272341, 8 pages, 2013. https://doi.org/10.1155/2013/272341
Proportional Fair Power Allocation for Secondary Transmitters in the TV White Space
The key bottleneck for secondary spectrum usage is the aggregate interference to the primary system receivers due to simultaneous secondary transmissions. Existing power allocation algorithms for multiple secondary transmitters in the TV white space either fail to protect the TV service in all cases or they allocate extremely low power levels to some of the transmitters. In this paper, we propose a power allocation algorithm that favors equally the secondary transmitters and it is able to protect the TV service in all cases. When the number of secondary transmitters is high, the computational complexity of the proposed algorithm becomes high too. We show how the algorithm could be modified to reduce its computational complexity at the cost of negligible performance loss. The modified algorithm could permit a spectrum allocation database to allocate near optimal transmit power levels to tens of thousands of secondary transmitters in real time. In addition, we describe how the modified algorithm could be applied to allow decentralized power allocation for mobile secondary transmitters. In that case, the proposed algorithm outperforms the existing algorithms because it allows reducing the communication signalling overhead between mobile secondary transmitters and the spectrum allocation database.
The main requirement for secondary spectrum usage is interference control. For single secondary transmitter, the power level maximizing the secondary transmission rate under primary system protection constraints has been derived in . However, the bottleneck for secondary spectrum usage is the aggregate interference to the primary receivers due to simultaneous secondary transmissions. Allocating the power levels to multiple secondary transmitters has already drawn considerable interest for standardization bodies all over the world [2–4] and also for the academic research community.
The SE43 working group in the Electronic Communications Committee (ECC) in Europe announced a power allocation rule for secondary operation in the TV white space (TVWS), hereafter referred to as the SE43 rule . Initially, the SE43 rule allocates the maximum permitted transmit power level to each secondary transmitter. For 2–4 simultaneous transmissions the maximum transmit power is divided with the number of active transmitters. In that case, each transmitter generates an equal amount of interference at the TV system. For a larger number of simultaneous transmissions a safety protection margin should be utilized. It has been shown that the existing SE43 rule fails to protect the TV service in all cases .
Power allocation algorithms that satisfy the TV protection constraints have been proposed by the academic research community [6–17]. In [6–10] all secondary transmitters allocated equal transmit power levels. In that case, there is a single parameter value to determine and the computational complexity is low. The drawback is that the power level for secondary transmitters located far from the TV coverage area would be limited due to secondary transmissions that are originated close to the TV coverage area. The assumption for common allocated power level is dropped in  where the sumpower is maximized. In that case, the algorithm associates transmitters located far from the TV coverage area with a high power level at the cost of transmitters in the proximity of the TV coverage area which remain practically silent.
A proportional fair (PF) power allocation rule would favor equally the secondary transmitters by equalizing their resource consumption. The available resource is the maximum permitted generated interference at the TV system. Note that the SE43 rule resembles a PF power allocation scheme in the limiting scenario with 2–4 transmitters. Unlike the SE43 rule, the scheme proposed in this paper is able to protect the TV service in all cases. When the number of secondary transmitters is high, the complexity of the proposed scheme may prohibit a real-time implementation in a spectrum allocation database. Because of that, a suboptimal method with low complexity is also introduced. The suboptimal method not only allows protecting the TV service but also enables power allocation to mobile users in a decentralized manner.
The ECC report  is not the first work targeting fairness in secondary type of networks. Fairness in secondary spectrum allocation and scheduling has been studied in [12–14, 18–21] respectively. Also, in  the minimum secondary SINR is maximized while in  the sum rate utility is maximized. In  PF power allocation is proposed but the complexity in large problem instances is not addressed. Unfortunately, the power allocation algorithms in [15–17] assume secondary transmitters with fixed and known locations. The interference that is originated from mobile users is not considered.
The remainder of this paper is organized as follows. In Section 2 we present the TV system model. In Section 3 we formulate the PF power allocation scheme under TV protection constraints. In Section 4 we propose a suboptimal alternative to the PF scheme that is extended to allow decentralized power allocation for mobile devices in Section 5. In Section 6 we demonstrate the flexibility of the proposed scheme in comparison with the current SE43 power allocation rule. In Section 7 we conclude the paper.
2. TV System Model
In TV network planning the location probability describes the percentage of locations within a square area of , also known as pixel, where the TV reception is satisfactory . For satisfactory TV operation a target SINR must be maintained with specific outage probability due to the slow fading. The outage probability is complementary to the location probability widely used in the definition of TV coverage contour . For evaluating the outage probability the TV coverage area is discretized into pixels. Each pixel is considered to be a test point where the following condition should be satisfied: where is the SINR at the th TV test point.
The probability constraint (1) must be satisfied in the presence of secondary transmissions too. A closed-form expression for (1) has been derived in  assuming that the TV signal level and the interference level due to each secondary transmitter follow the log-normal distribution within a TV test pixel. The aggregate secondary interference becomes a sum of log-normal random variables whose distribution is not available and approximating methods are commonly employed. In  the distribution of the aggregate secondary interference incorporating also the noise level is approximated by the log-normal distribution. The Fenton-Wilkinson method is used to compute the parameters of the approximating distribution.
Expressing (1) in closedform allows computing the maximum permitted mean secondary interference level which has been defined in  as the interference margin. The interference margin is a function of the TV system parameters, that is, the protection criteria of TV receivers () and the useful TV signal level at the test point. Also, it depends on the locations of secondary transmitters and their transmission power levels. In order to simplify the interference control process a lower bound for the interference margin has been derived in  that is independent of the secondary network parameters: where is the available interference margin at the th TV test point, is the TV signal level by using distance-based pathloss, in dB is the standard deviation of TV signal due to slow fading, is the inverse of the Gaussian function, and stands for the noise power level at the TV receiver.
By using the concept of interference margin one essentially turns the chance type of constraint (1) into a constraint of the form where is the number of secondary transmitters and is the mean interference level from the th transmitter to the th TV test point. The mean generated interference is a function of the secondary transmission power level , the distance-based propagation pathloss , and the mean of the slow fading. For log-normal fading the mean interference level is where in dB is the slow fading standard deviation due to the secondary transmission and is a scaling constant. The standard deviation is in general different for different secondary transmitters and TV test points due to the irregular terrain morphology.
3. Proportional Fair Power Allocation Scheme
By using the concept of interference margin the probability constraint (1) is turned into a linear constraint where is the link gain. Therefore the interference margin can be treated as an available resource. Each transmitter is allowed to take a bite out of the available resource; that is, it is allowed to generate some amount of interference at the TV test point provided that the mean interference of secondary transmissions does not exceed the interference margin.
Allocating the transmit power levels can be viewed as a resource sharing problem. In order to equalize the resource consumption we propose a PF power allocation scheme. For PF power allocation the sum of secondary transmission power levels in the log-domain should be maximized . As a result, the power allocation algorithm can be formulated as a constrained optimization problem:
The optimization function and the constraints are continuously differentiable functions of . The Lagrangian function is where is the Lagrangian multiplier. The Karush-Kuhn-Tucker conditions are
The solution has the form where the Lagrangian multiplier is positive if the th constraint is binding. Otherwise, . The PF scheme allocates the transmit power levels so that each transmitter takes an equal fraction out of the aggregate interference margin at the binding points. Unfortunately, for complex secondary network geometries and noncontiguous TV coverage areas, it is difficult to identify the binding points in advance.
One can show that the Hessian matrix of the objective function is negative definite and thus the optimization function is concave. Since the constraints are linear functions of transmit power levels, the optimization problem (4) is convex. The globally optimal solution can be found by using standard convex optimization tools. When the number of constraints and optimization parameters are high, the computational complexity may prohibit a real-time computation in the spectrum allocation database. Note that the database may need to reallocate the transmission power levels in frequent time intervals. The secondary transmitters can be mobile or new transmitters can send spectrum access requests.
4. Simplified Proportional Fair Power Allocation Scheme
We propose associateing each transmitter with the TV test point where its generated interference is maximized. In this way, we end up with a subset of TV test points that experience the highest interference levels. We call it the set of dominating TV test points. Obviously, it is sufficient to control the generated interference only at the TV test points belonging to the dominating set.
We denote by the set of transmitters dominated by the TV test point. Also, let us assume that each transmitter generates mean interference level equal to at its dominating point, . Essentially, we propose a PF power allocation scheme for secondary transmitters belonging to the same set. Next, we show how to compute the parameter .
The generated interference at the th TV test point due to the secondary transmitter can be written as a fraction of the interference : The mean interference at the th TV test point is the sum of the interference from all secondary transmitters:
Since each transmitter generates interference equal to at the test point, where denotes the cardinality of a set. Also, by using (6) in (7) the aggregate interference can be read as where + exp. For protecting the TV receivers the mean interference should be kept under the margin, . By enforcing the inequality to be tight, the parameter can be selected as the minimum over all test points belonging to the dominating set:
With the parameter at hand, we can identify the transmit power level after replacing by into (3): where the value for the distance-based pathloss and the standard deviation must be taken from the test point that dominates the th transmitter. By using (8) and (2) the allocated power level can also be written as
In Figure 1 it is assumed that both TV test points belonging to the dominating set are also binding; that is, all their interference margin is filled in by the secondary transmissions. According to the simplified PF scheme each transmitter generates mean interference at its dominating test point. The remaining interference headroom is filled in by the secondary transmitters dominated by the other test point. According to the PF scheme the aggregate interference margin is divided equally among the secondary transmitters.
(a) Simplified proportional fair power allocation
(b) Proportional fair power allocation
The two schemes are characterized by different implementation complexity. The complexity to identify the set of dominating test points in the simplified scheme is while the complexity to identify the minimum over the computed values for is . Since the overall complexity can be approximated by . In order to get insight on the complexity for solving (4) by using interior-point methods we use the complexity bound . Besides the lower implementation complexity, the simplified scheme allows power allocation to mobile UE with lower signalling overhead in comparison with the PF scheme. This is the topic of the next section.
5. Power Allocation for Mobile Users
The power allocation scheme proposed in Section 4 is applicable to secondary transmitters with fixed and known locations, for example base stations (BSs). The secondary transmitters can also be mobile user equipment (UE). If the number of UE is high, the communication signaling overhead for updating their locations in the database at frequent time intervals would be high too. Fortunately, the simplified PF algorithm can be modified to allow decentralized power allocation with reduced signaling overhead. Next, we show how the simplified scheme can group secondary pixels and control the interference from each group instead of doing it on a per-pixel basis. UE can move within pixels belonging to the same group without the need to inform the database about their locations.
Following the ECC approach, we discretize the secondary deployment area into pixels. We assume that the UE inside each pixel is distributed according to a Poisson point process (PPP). The density of the process in the th pixel is equal to the number of UE times the activity factor . The mean interference due to the transmissions originating from the th pixel to the th TV test point is  where the index is used to integrate over the area of the th pixel, is the transmit power level for UE located inside the th pixel, and where is the integration element. The describes the average number of active UE inside the integration element of the th pixel. The stands for the distance-based pathloss between the th integration element of the th pixel and the th test point.
Similar to the power allocation rule described in Section 4, we group together pixels dominated by the same test point. The mean interference level from the group to the th TV test point is . (Note that the parameter value and the set of dominating TV test points are in general different for the secondary transmitters with fixed locations described in Section 4 and for the secondary pixels described in this section. In order not to introduce unnecessary complexity we keep the same notation and provide the related explanation when needed.) The interference headroom can be divided among the secondary pixels belonging to the group by using any rule. For PF allocation each pixel takes an equal amount of resource, . In that case, pixels with low user density allocate higher power levels to the UE compared to pixels with high user density. One alternative way is to divide the interference headroom proportionally to the user density of pixels belonging to the set . In that case, the transmission power level for UE inside the th pixel becomes
The parameter can be calculated based on (8) after computing the parameter as
The transmission power level calculated based on (11) does not depend on the particular population density in the th pixel. Only the total amount of UE inside the group has to be known. Because of that, the UE can freely move between pixels belonging to the same group without the need to recalculate their transmission power levels at the database. Each UE just has to know the pixel where it belongs and set its transmission power level based on (11). The database has to recompute the parameter value and the transmission power levels only if the total number of UE in a group of pixels changes.
Allowing UE mobility within a group may violate the protection criteria at the th TV test point. One way to protect the th test point is to recompute the parameter in (12) based on the worst case interference scenario. That refers to the case where all UE inside the group is concentrated in the pixel generating the highest interference at the th TV test point. In that case the parameter becomes
Thus far, the transmission power levels for BS and UE have been calculated assuming that the full interference margin is available to BS and UE, respectively. When the BS and the UE are simultaneously active, the power allocation algorithm may have the following steps: (1)allocate some fraction of the interference margin to the BS, ,(2)compute the parameter for the BS transmissions after replacing by in (8), (3)compute the BS transmission power levels by using (9),(4)by using the identified BS transmission power levels compute the remaining interference margin at the TV test points where the is calculated from (3),(5)group the secondary pixels and identify the transmission power level for UE in each pixel by using (11). The parameter is computed after replacing by in (8). For supporting UE mobility within groups of pixels (13) must be used instead of (12).
6. Numerical Illustrations
6.1. Parameter Settings
For testing the proposed power allocation rule we select a study area in Jyväskylä in Finland where the Vihtavuori TV transmitter operates at MHz. The minimum required median field strength for fixed TV reception at QAM and code rate is dbuV/m at MHz [26, pp.185]. For other frequencies in MHz the correction factor should be added resulting in dBuV/m at MHz.
In order to identify the TV coverage area we cover the full study area with square pixels with side equal to m. A pixel belongs to the TV coverage area if the median field strength is higher than dBuV/m. The TV field strength at the center of each pixel has been calculated by using the Longley-Rice propagation model with terrain data implemented in Splat! . The slow fading standard deviation of the TV signal is dB . The TV self-interference is not considered and the noise power level is dBm. The protection criteria of TV receivers are dB and motivated from .
For testing power allocation for mobile UE, the population density map of the area is utilized; see Figure 2. The antenna height for the UE is m. Their activity factor is %. The locations of secondary BS are fixed and known. They are placed at a square lattice with km side and their antenna height is m. For modeling the distance-based propagation pathloss of secondary transmissions we utilize the modified Hata model for suburban areas . The slow fading standard deviation for BS and UE is dB .
First, we compare the allocated transmit power levels of secondary BS by using PF and simplified PF power allocation rules. For PF rule the power levels are obtained by solving (4) while for simplified PF rule (9) is utilized. For illustration purposes the dominating TV test points and the locations of secondary BS are depicted in Figure 3 for protection distance equal to km. The distributions of BS power levels for different protection distances are plotted in Figure 4. One can see that the simplified scheme almost achieves the optimal solution. In general, the PF scheme allocates higher transmission power levels for the BS located close to the TV protection area. Also, we notice that the set of dominating test points has different cardinality for different protection distances and only a single test point is binding. Hereafter, the simplified PF scheme is utilized and it is referred to as the proposed power allocation method.
Next, we compare the proposed power allocation method with other existing methods. If we are to use equal transmit power for all secondary transmitters [6–10], the maximum permitted levels are dBm for protection distances km, respectively. As a result, the common power allocation rule is conservative for secondary transmitters located far from the TV coverage area especially when the protection distance is small. Also, if the sumpower is maximized , a single transmitter is essentially favored: the transmitter located furthest away from the TV coverage area has the smallest link gain and can dominate the sum-power utility by using high transmit power level. The power level for the rest of the transmitters is practically zero. As a result, the sum-power utility is extremely unfair in our system setup.
Next, we compare our method with the existing SE43 proposal. The transmission power level for each BS by using the SE43 rule is [2, page 24] where and are the safety and the multiple interference protection margin, respectively.
Instead of using protection margins one can scale the transmission power of each BS with the total number of active BS. The modified SE43 rule allocates transmission power level equal to 
The TV SINR distribution at the test point experiencing the highest interference level is depicted in Figure 5 for protection distance km. The protection margins for the SE43 rule are taken equal to dB and dB as proposed in the ECC report . One can see that the existing and the modified SE43 rules fail to protect the TV service. The same behaviour is observed for other protection distances too. One could maintain the TV SINR target by using, for instance, a higher value for the safety margin in (14). This approach requires an iterative type of power allocation algorithm implemented in the database. On the other hand, the proposed rule has low complexity and protects the TV system in all cases.
The proposed rule can also be used to allocate the transmission power levels to the UE (11). UE mobility between pixels belonging to the same group is allowed if (13) is used instead of (12). This flexibility comes at the cost of reduced transmission power level at the UE. However, the loss is negligible especially for a large protection distance; see Figure 6.
Finally, we consider simultaneous BS and UE transmissions and the power allocation algorithm described in Section 5. In Figure 7 one can see the distribution of UE transmit power levels for different fractions of allocated interference margin to the BS. When of the margin is allocated to the BS its transmission power levels are between dBm and dBm while the UE power levels lie between dBm and dBm. At the same time, the TV SINR target is satisfied at all the dominating test points.
In this paper we proposed a power allocation algorithm for multiple secondary transmitters in the TVWS that protects the TV service in the cochannel. The proposed rule can incorporate secondary transmitters with different antenna heights, for instance, base stations and mobile UE. The impact of terrain morphology can be considered by using different slow fading standard deviations for transmitters located at different pixels. The proposed rule allows reducing the communication signaling overhead between the UE and the spectrum allocation database with the reason being that the database groups secondary pixels and controls the aggregate interference from each group instead of doing it on a per pixel basis. The proposed rule outperforms the existing SE43 rule because it protects the TV service in all cases and results in smaller communication signalling overhead. However, it does not consider the scenario where some of the secondary transmitters fall inside the coverage area of an adjacent TV channel. The SE43 rule takes into account the protection of adjacent channel TV service by using reference coexisting geometries. The reference geometry rule is conservative particularly at low user densities . Extending the proposed rule to incorporate adjacent channel TV protection can be an area of future work.
This work was partially supported by the TEKES-funded project IMANET+ and by the Academy-of-Finland-funded project SMACIW under Grant no. 265040.
- X. Kang, R. Zhang, and Y. C. Liang, “Optimal power allocation for cognitive radio under primary user’s outage loss constraint,” in Proceedings of the IEEE International Conference on Communications (ICC '09), Dresden, Germany, June 2009.
- ECC, “Report (under public consultation): technical and operational requirements for the possible operation of cognitive radio systems in the white spaces of the frequency band 470–790 MHz,” Tech. Rep. 159, 2011.
- Ofcom, “Digital dividend: geolocation for cognitive access, discussion on using geolocation to enable license-exempt access to the interleaved spectrum,” November 2009.
- “In the Matter of Unlicensed Operation in the TV Broadcast Bands: Third Memorandum Opinion and Order, Federal Communications Commission, FCC 12-36A1,” August 2012, http://transition.fcc.gov/Daily_Releases/Daily_Business/2012/db0405/FCC-12-36A1.pdf.
- R. Jäntti, J. Kerttula, K. Koufos, and K. Ruttik, “Aggregate interference with FCC and ECC white space usage rules: case study in Finland,” in Proceedings of the IEEE International Symposium on Dynamic Spectrum Access Networks (DySPAN '11), pp. 599–602, May 2011.
- M. Vu, N. Devroye, and V. Tarokh, “On the primary exclusive region of cognitive networks,” IEEE Transactions on Wireless Communications, vol. 8, no. 7, pp. 3380–3385, 2009.
- N. Hoven and A. Sahai, “Power scaling for cognitive radio,” in Proceedings of the Wireless Networks, Communications and Mobile Computing, pp. 250–255, June 2005.
- S. Zaidi, D. Mclernon, M. Ghogho, and B. Zafar, “On outage and interference in 802.22 cognitive radio networks under deterministic network geometry,” in Proceedings of the ACM Mobicom, Beijing, Japan, September 2009.
- S. Shankar and C. Cordeiro, “Analysis of aggregated interference at DTV receivers in TV bands,” in Proceedings of the 3rd International Conference on Cognitive Radio Oriented Wireless Networks and Communications (CrownCom '08), Singapore, May 2008.
- E. Larsson and M. Skoglund, “Cognitive radio in a frequency-planned environment: some basic limits,” IEEE Transactions on Wireless Communications, vol. 7, no. 12, pp. 4800–4806, 2008.
- K. Koufos, K. Ruttik, and R. Jäntti, “Controlling the interference from multiple secondary systems at the TV cell border,” in Proceedings of the IEEE Personal Indoor and Mobile Radio Communications (PIMRC '11), pp. 645–649, Toronto, Canada, September 2011.
- J. Tang, S. Misra, and G. Xue, “Joint spectrum allocation and scheduling for fair spectrum sharing in cognitive radio wireless networks,” Computer Networks, vol. 52, no. 11, pp. 2148–2158, 2008.
- B. Wang and D. Zhao, “Scheduling for long term proportional fairness in a cognitive wireless network with spectrum underlay,” IEEE Transactions on Wireless Communications, vol. 9, no. 3, pp. 1150–1158, 2010.
- L. Akter and B. Natarajan, “Modeling fairness in resource allocation for secondary users in a competitive cognitive radio network,” in Proceedings of the Wireless Telecommunications Symposium (WTS '10), pp. 1–6, April 2010.
- L. Tang, H. Wang, and Q. Chen, “Power allocation with max-min fairness for cognitive radio networks,” in Proceedings of the Global Mobile Congress (GMC '10), October 2010.
- B. Cho, K. Koufos, K. Ruttik, and R. Jäntti, “Power allocation in the TV white space under constraint on secondary system selfinterference,” Journal of Electrical and Computer Engineering, vol. 2012, Article ID 245895, 12 pages, 2012.
- D. Li and J. Gross, “Distributed TV spectrum allocation for cognitive cellular network under game theoretical framework,” in Proceedings of the Dynamic Spectrum Access Networks (IEEE DySPAN '12), October 2012.
- H. Zheng and C. Peng, “Collaboration and fairness in opportunistic spectrum access,” in Proceedings of the IEEE International Conference on Communications (ICC '05), Seoul, Korea, May 2005.
- Y. Xing, R. Chandramouli, S. Mangold, and S. Nandagopalan, “Dynamic spectrum access in open spectrum in wireless networks,” IEEE Journal on Selected Areas in Communications, vol. 24, no. 3, pp. 627–637, 2006.
- J. Bater, H. P. Tan, K. N. Brown, and L. Doyle, “Modelling interference temperature constraints for spectrum access in cognitive radio networks,” in Proceedings of the IEEE International Conference on Communications (ICC '07), pp. 6493–6498, June 2007.
- Z. Ma, Z. Cao, and W. Chen, “A fair opportunistic spectrum access (FOSA) scheme in distributed cognitive radio networks,” in Proceedings of the IEEE International Conference on Communications (ICC '08), pp. 4054–4058, May 2008.
- K. Ruttik, K. Koufos, and R. Jäntti, “Modeling of the secondary system's generated interference and studying of its impact on the secondary system design,” Radioengineering, vol. 19, no. 4, pp. 488–493, 2010.
- F. P. Kelly, A. K. Maulloo, and D. Tan, “Rate control for communication networks: shadow prices, proportional fairness and stability,” Journal of the Operational Research Society, vol. 49, no. 3, pp. 237–252, 1998.
- A. Nemirovski, Interior Point Polynomial Time Methods in Convex Programming, Lecture Notes, Georgia Institute of Technology, 2004.
- K. Koufos, K. Ruttik, and R. Jantti, “Aggregate interference from WLAN in the TV white space by using terrain-based channel model,” in Proceedings of the IEEE Cognitive Radio Oriented Wireless Networks and Communications (CrownCom '12), Stockholm, Sweden, June 2012.
- Final Acts of the Regional Radiocommunication Conference for Planning of the Digital Terrestrial Broadcasting Service in Parts of Regions 1 and 3, in the Frequency Bands 174–230 MHz and 470–862 MHz (RRC-06), ITU-R, 2006.
- “RF Signal Propagation, Loss, and Terrain analysis tool for the spectrum between 20 MHz and 20 GHz,” http://www.qsl.net/kd2bd/splat.html.
- “Monte carlo radio simulation methodology,” ERC Report 68, Naples, Italy, 2000.
- “SE43(11)63 draft report. Method of IM value setting in ECC Report 159,” http://www.cept.org/ecc/groups/ecc/wg-se/.
- L. Shi, K. W. Sung, and J. Zander, “Controlling aggregate interference under adjacent channel interference constraint in TV white space,” in Proceedings of the IEEE Cognitive Radio Oriented Wireless Networks & Communications (CrownCom '12), Stockholm, Sweden, June 2012.
Copyright © 2013 Konstantinos Koufos and Riku Jäntti. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.