The leaf economic traits such as leaf area, maximum carbon assimilation rate, and venation are all correlated and related to water availability. Furthermore, leaves are often broad and large in humid areas and narrower in arid/semiarid and hot and cold areas. We use optimization theory to explain these patterns. We have created a constrained optimization leaf model linking leaf shape to vein structure that is integrated into coupled transpiration and carbon assimilation processes. The model maximizes net leaf carbon gain (NPPleaf) over the loss of xylem water potential. Modeled relations between leaf traits are consistent with empirically observed patterns. As the results of the leaf shape-venation relation, our model further predicts that a broadleaf has overall higher NPPleaf compared to a narrowleaf. In addition, a broadleaf has a lower stomatal resistance compared to a narrowleaf under the same level of constraint. With the same leaf area, a broadleaf will have, on average, larger conduits and lower total leaf xylem resistance and thus be more efficient in water transportation but less resistant to cavitation. By linking venation structure to leaf shape and using water potential as the constraint, our model provides a physical explanation for the general pattern of the covariance of leaf traits through the safety-efficiency trade-off of leaf hydraulic design