Yukawa-unified natural supersymmetry

18 pages, 7 figures, pdflatexPrevious work on t-b-\tau Yukawa-unified supersymmetry, as expected from SUSY GUT theories based on the gauge group SO(10), tended to have exceedingly large electroweak fine-tuning (EWFT). Here, we examine supersymmetric models where we simultaneously require low EWFT ("natural SUSY") and a high degree of Yukawa coupling unification, along with a light Higgs scalar with m_h\sim125 GeV. As Yukawa unification requires large tan\beta\sim50, while EWFT requires rather light third generation squarks and low \mu\sim100-250 GeV, B-physics constraints from BR(B\to X_s\gamma) and BR(B_s\to \mu+\mu-) can be severe. We are able to find models with EWFT \Delta\lesssim 50-100 (better than 1-2% EWFT) and with Yukawa unification as low as R_yuk\sim1.3 (30% unification) if B-physics constraints are imposed. This may be improved to R_yuk\sim1.2 if additional small flavor violating terms conspire to improve accord with B-constraints. We present several Yukawa-unified natural SUSY (YUNS) benchmark points. LHC searches will be able to access gluinos in the lower 1-2 TeV portion of their predicted mass range although much of YUNS parameter space may lie beyond LHC14 reach. If heavy Higgs bosons can be accessed at a high rate, then the rare H, A\to \mu+\mu- decay might allow a determination of tan\beta\sim50 as predicted by YUNS models. Finally, the predicted light higgsinos should be accessible to a linear e+e- collider with \sqrt{s}\sim0.5 TeV